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J. Hadamard's ideas of correct formulation of problems of mathemati- 
cal physics as well as related Banach's theorem on the inverse operator are 
analyzed. Modern techniques of numerical simulations are shown to be in 
drastic contradiction to the concepts of J. Hadamard, S. Banach and a num- 
ber of other outstanding scientists in the sense that the priority is given 
to the realization of inefficient algorithms, based on a belief that ill-posed 
problems are adequate to real phenomena. 

A new method of the solution of problems, traditionally associated with 
Fredholm integral equations of the first kind, is developed. Its key aspect 
is a constructive use of possibilities of the functional space I2 to ensure 
the conditions of correctness. A well-known phenomenon of smoothing of 
information is taken into account by means of a special composition that 
explicitly involves the sought function and is infinitesimal in the space L2. 
A finite error in the determination of the data by measurements can also 
be represented in the same manner. By relatively simple transformations, 
the outlined class of problems is reduced to the solution of Fredholm integral 
equations of the second kind with properties most favorable for the numerical 
realization. 

We demonstrate a reduction to Fredholm integral equations of the first 
kind and, correspondingly, a possibility to extend the suggested approach 
to wide classes of linear boundary-value and initial-boundary-value prob- 
lems characterized by variable coefficients, non-canonical domains, as well 
as by other factors complicating their solution. Basically similar techniques 
arc developed for the solution of initial-boundary-valuc problems for non- 
linear differential equations of the evolutionary type. Also considered are 
boundary-value problems for substantially nonlinear equations and equa- 
tions of a mixed type, the cases of nonlinearities in the boundary conditions, 
the presence of a small factor by the highest-order derivative, the inverse 
problem of the restoration of the variable coefficient of the differential oper- 
ator, and the problem of the Stefan type. 

We put forward arguments that the determination of causal relation- 
ships, based on the formulation restricted to a primitive renaming of known 
and unknown functions of the corresponding direct problem, is essentially 
illegitimate. The suggested approach may be considered as a constructive 
realization of J. Hadamard's ideas of the existence of correct formulations of 
physically meaningful problems. 
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Chapter 1 
Introduction 



At the beginning, we should explain the title of the work and, in the first 
place, the meaning of the employed notions. In this regard, we assume 
the availability of information allowing us to formulate a mathematical 
model of a certain phenomenon in a traditional way. Correspondingly, 
the determination of unknown functions using the data of the problem 
is implied. If the dependence of the solution to the problem on these 
data with respect to the norm of the chosen space is continuous, such a 
problem, as a rule, belongs to the domain of analysis or, in other words, 
its formulation is direct. 

However, the investigation of a concrete phenomenon in a variety 
of the determining factors with the aim of obtaining, as a final re- 
sult, of qualitatively new information (the synthesis of knowledge) also 
envisages the realization of problems in their inverse formulation, i.e., 
the restoration of data using the hypothetically known solution: In 
other words, the restoration of the cause using its consequence, which 
is usually identified with the necessity of solving incorrectly formulated 
problems. 

The purpose of the present investigation consists in the justification 
of the illegitimacy of this statement and, on the contrary, in a construc- 
tive development of J. Hadamard's ideas of the existence of correctly 
formulated problems, adequately describing real processes and phenom- 
ena. Note that the difference between these two notions in the context 
of the book is unessential. However, the term "process" accentuates a 
time factor. 
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CHAPTER 1. INTRODUCTION 



In the focus of the attention is a natural, to our mind, issue that, 
as an example, can be explained by the evaluation of the integral 

1 

(AV) {x) = Jk {x, (0 dC^fix), xe [0, 1] , (1.1) 



which amounts to the determination of the function / (x) using given 
k (a;, ^) and ip (x). This procedure can be easily associated with a lot 
of physical, as well as other, interpretations. Its realization, at least in 
the case of the bounded integrand, does not pose any problem. 

On the other hand, if the kernel k {x, ^) and the right-hand side 
/ (x), evaluated according to Eq. (1.1), are given, it is obvious that the 
corresponding function ijj (x) is objectively existent. Thus, the problem 
is whether it is legitimate to restore this function by means of the 
solution of the Frcdholm integral equation of the first kind (1.1), just 
renaming the known and unknown components, i.e., by assuming that 
the function / (x) is given and ip (x) is to be determined. And, generally 
speaking, what is the basis to argue that mathematical formulations of 
the direct the inverse problems can be absolutely identical? 

By the way, it is very difficult to imagine the way of the formulation 
of the inverse problem on physical grounds explicitly, without any rela- 
tion to the direct one. The fact of primitive renaming of the known and 
the unknown function, without any additional corrections, itself causes 
a protest. Thus, wc put forward a thesis that a constructive approach 
to the formulation of inverse problems should differ from the common 
one. This orientation caused the presence in the title of the work the 
notion of methodology. 

As a matter of fact, we hope to find reserves of the synthesis of the 
whole complex of knowledge about the phenomenon by investigating it 
from different sides, using formulations whose mathematical definitions 
are not identical. Although the use of Eq. (1.1) for the determination 
of / (x) is fully justified, the problem of the restoration of the function 
■0 (x) does not need to coincide with it with regard to its representation 
and, correspondingly, does not need to be reducible to the Fredholm 
integral equation of the first kind, characterized by instability of algo- 
rithms of numerical realization. 
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However, if there exists an alternative to the formulation of the in- 
verse problems by means of the above-mentioned renaming of known 
and unknown components, it is logical to assume that corresponding 
formulations may possess much more attractive properties in a compu- 
tational sense. Prom this point of view, the arguments of J. Hadamard 
acquire a rather concrete meaning, stimulating a search of adequate, 
in the given situation, and, at the same time, correct formulations of 
problems of mathematical physics. A constructive realization of the 
outlined orientation turned out to be possible in the context of the 
following considerations. 

The reasons for the difficulties related to the solution of incorrectly 
formulated and essentially mathematically senseless problems are, in 
principle, well understood. As regards the Fredholm integral equation 
of the first kind (1.1), the key role is played by the mismatch between 
the function / (x) and a hypothetically exact solution of the corre- 
sponding direct problem (the result of integration) in the conditions of 
approximations and different rounding of significant digits. 

Correspondingly, considerable attention is paid to the so-called phe- 
nomenon of smoothing of information in procedures of integration. At 
the same time, the data of the problem, i.e., the free term f (x) and 
also the kernel k{x,$,), are usually determined experimentally, which 
inevitably incurs a considerable error in Eq. (1.1). In this regard, we 
should point out the dominance of the methodology of A. N. Tikhonov 
that is based on objective incorrectness of the formulation of most prob- 
lems of mathematical simulation. 

There appears a rather obvious, as it seems, question: Why not take 
into account in practice the above-mentioned errors in the formulation 
of problems, instead of merely bearing them in mind when identifying 
the reasons for computational discrepancies? One can assume that an 
adequate simulation of the error may contribute to a correct formulation 
of problems that are considered to be inverse. 

The adequacy implies, in the first place, the structure of the rep- 
resentation of the error. In this regard, let us turn to the procedure 
of integration (1.1). On the basis of general considerations, it is logi- 
cal to represent the loss of information about the function ip {x) in the 
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evaluation of / {x) in the form 



{Sf){x)^^{x)-X lh{x,0^{Od^, xe[Q,l]. 



(1.2) 



-1 



Here, the function ip (x), x E [—1,0), the kernel h {x, C,) and the param- 
eter A should satisfy the requirement of the realization of the condition 
of the type {5f) (x) = in L2 (0, 1) for -0 (x), x e [0, 1] from a rather 
representative class. As will be shown, this requirement does not pose 
a considerable problem. 

Thus, the outlined structure embodies the difference between the 
function 4'{x), subject to integration, and a close expression that, in 
turn, appears as a result of the execution of an analogous procedure. 
One should note the absence of any a priori premises of self-sufficiency 
of (1.2) in achieving the goal, namely, a correct formulation of the 
problem of the restoration of the function ip {x) on the basis of the 
data of (1.1). In other words, we rather put forward a hypothesis about 
priority importance of a qualitative side of the model of smoothing of 
information. 

In light of the above, instead of the traditional incorrectly formu- 
lated problem (1.1) for the determination of the function ijj{x), the 
following system of equations is employed: 



where /x is a parameter, analogous to A. 

As a result of comparatively simple transformations, we managed 
to reduce the considered problem to the solution of the Predholm in- 
tegral equation of the second kind with more favorable properties for 
numerical realization. We assume the absence of singularities as well 
as oscillations of the kernel, not caused by k (x, ^), i.e., those that are 
enforced by the employed algorithm. 

Aside from explicit presence of the sought function and, in gen- 
eral, the structure of the representation of the error (1.2), a necessary 
condition turned out to be the completeness of the kernel 



f,{A^){x)^f,f{x) + {Sf){x); {Sf){x)^0, x e [0, 1] , 



h{x,0 



1 — 2r cos [27r {x — ^) + r^] ' 



< |r| < 1, 
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whereas the simphfication of the evaluation proved to be possible due 
to its symmetry h{x) = h{x ±1) and dependence on the difference of 
the arguments. 

However, a broad class of linear boundary- value and initial-boundary- 
value problems of mathematical physics can be rather elementarily re- 
duced to the Frcdholm integral equations of the first kind. As a re- 
sult, the above-mentioned procedure can be considered as an efficient 
method of their numerical reahzation. This statement is illustrated by 
the following example: 



d^u + dyU 



-1, ti = 0, X = y = Q\ X = y = 1. 
Form the notation d^u — ip, in view of (1.3), it follows: 



:i.3) 



u {x, y) 



'X X 

Jix-O-xJil-O 



(1.4) 



u{x,y) = -y{l-y) 



y y 

j {y-v)-y j i^-v) 



ijj{x,r])dr]. (1.5) 



The elimination of the function u {x, y) leads to a two-dimensional Fred- 
holm integral equation of the ffrst kind. In the case of conditions on the 
contour of a more general type, one can use, for example the following 
relation: d^u + u — t/j. The outhned scheme is practically indifferent 
to the type and order of differential operators, the presence of variable 
coefficients, the configuration of the boundary of the domain of the 
function and some other factors that usually complicate the realization 
of numerical algorithms. 

It should be noted that the dimensionality of the boundary- value 
problem, regarding its further reduction to the Fredholm integral equa- 
tion of the first kind, is of no principal importance, because this proce- 
dure is, in fact, performed with respect to one variable while the rest 
play the role of particular parameters. By the way, in this context, 
the conventional division of problems of mathematical simulation into 
direct and inverse becomes rather conditional. 

Nonlinear boundary- value and initial-boundary-value problems are, 
in general, analogously reduced to nonlinear integral equations of the 
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second kind that are specially adopted for an effective realization of 
iterative algorithms. 

Perhaps, the motivation of the proposed investigation could be of 
certain interest. At the beginning, the author was confused by the ab- 
sence in the specialized literature of a clear statement of the universality 
of the presented method of the reduction of problems of mathematical 
physics to the Fredholm integral equations of the first kind. Given the 
complexity of the solution of incorrect problems, the factor of the elim- 
ination of boundary and initial conditions as well as a placement of 
the whole lot of initial data in the kernels of derived equations seems 
to be rather attractive. Indeed, a conventional classification of prob- 
lems according to the complexity of their numerical realization, in fact, 
vanishes, and the construction of an effective algorithm of the solu- 
tion of the Fredholm integral equation of the first kind appears in the 
foreground. 

At the same time, a specific degeneracy of the relations for the 
new unknown function ip, generating the representation (1.4), (1.5), is 
worth noting. There appears a natural question: Is there a possibility 
to reduce boundary- value problems to equations with more acceptable 
properties (in the sense of a numerical realization) by means of relations 
between the functions u and -0 of a more general type. For example, 
setting 



can one complete the definition of the kernel K (,x, y, ^) to reduce the 
problem (1.3) to the Fredholm integral equation of the second kind with 
respect to the function ip7 However, an attempt to do something like 
this did not succeed. 

Analogously, no results were yielded by attempts to reduce the con- 
sidered problems to integral equations of the first kind Btp = Ai/j + /, 
with the operator B having a suitable inversion form [by analogy with 
(1.1), A and / characterize the data of a concrete problem]. At the 



X 



X 




Lo 







where 



K{x,y,x) = 0, x,ye[Q,l], 
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same time, this did not mean that one could not reahze, in princi- 
ple, the algorithm of a derivation by means of transformations of the 
Predholm integral equation of the second kind. 

The use of the model of smoothing of information by means of 
integration of the type (1.2) proved to be constructive. In this regard, 
a profound elaboration of methodological aspects of correct formulation 
of problems of mathematical physics proved to be necessary. 

In section 2, we analyze J. Hadamard's arguments concerning the is- 
sue of correct formulation of problems for partial differential equations. 
Both related and alternative positions on this issue of known special- 
ists are also illuminated. Different interpretations of contcxtually close 
Banach's inverse operator theorem, one of the three main principles of 
functional analysis, are given. 

Section 3 contains a critical review of methods of the solution of 
incorrect problems, based, in the first place, on the concepts of A. 
N. Tikhonov and V. M. Fridman. We express considerations on the 
reasons for a crisis in the field of mathematical simulation, related to 
an orientation at the use of incorrect formulations that reduce to a 
numerical realization of ill-defined algebraic systems. 

The material of section 4, in a sense refracts principle difficulties, 
accompanying the solution of incorrectly formulated problems, by the 
prism of fundamental concepts of J. Hadamard and S. Banach. In 
this context, we formulate general premises for correct formulation of 
problems of mathematical simulations, traditionally associated with the 
Fredholm integral equations of the first kind. 

Section 5 is devoted to the presentation of an original approach to 
the solution of a canonical equation of the above-mentioned type and is 
basic constructively. Here we investigate issues related to the simulation 
of the error and to the transformation of the initial formulation. We 
give an algorithm of the reduction of the problem to the Fredholm 
integral equation of the second kind, and also consider corresponding 
aspects of methodological character. 

In section 6, we concretely demonstrate the universality of the tech- 
nique of the reduction of linear boundary-value and initial-boundary- 
value problems to Fredholm integral equations of the first kind with 
respect to one of high-order derivatives. We also demonstrate an ex- 
tension to them of the algorithm of the previous section. 
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Section 7 develops the outlined orientation involving into the sphere 
of effective realization sufficiently nontrivial applications (including fac- 
tors of nonlinearity, singular perturbations and some other). The pre- 
sentation of the material has the form of sketches. 

In section Conclusions, we summarize the main points of the pro- 
posed investigation from the position of priority importance of correct 
formulation of problems of mathematical physics for their numerical 
realization. 

Mathematical techniques employed in the presentation of the ma- 
terial is comparatively simple: basics of the classical theory of integral 
equations; elements of functional analysis; general principles of for- 
mulations of problems of mathematical physics and methods of their 
solution. 

The hterature to each section is given in reference order. We refer 
mainly to original publications and those ones that treat the subject in 
most detail. Sections and subsections (chapters and paragraphs of the 
literature sources) are denoted, respectively, as section 1, section 1.1, 
sections 1.1, 1.2. 

The author has pubhshed about 30 papers on the subject of this 
investigation in leading scientific journals and conference proceedings. 
All these publications are written by the author on his own. They 
clearly state the general orientation of the investigation and the means 
of its practical realization. However, only by turning to the method- 
ology of correct formulation of problems of mathematical physics, the 
author managed to provide the justification of the developed approach. 

The concepts, described below, were formed at different stages un- 
der the influence of well-known specialists, in particular: Profs. S. M. 
Belonosov, P. M. Varvak, E. S. Wentzel, N. A. Virchenko, Yu. V. Gan- 
del, Yu. A. Danilov, V. A. Dobrushkin, M. D. Dolberg, N. T. Zhadrasi- 
nov, P. A. Zhilin, B. Y. Kantor, V. V. Kopasenko, K. A. Luric, S. G. 
Mikhlin, N. F. Morozov, M. V. Paukshto, A. S. Sakharov, E. A. Simson, 
V. C. Ryabenkii, L. A. Filshtinskii, N. P. Fleishman, D. G. Khlcbnikov, 
A. V. Cherkaev, I. Y. Ghudinovich, I. D. Chueshov, V. A. Shterbina, 
as well as Corresponding Member of The Academy of Sciences of the 
USSR A. I. Luric, Academician of the National Academy of Sciences of 
Ukraine V. L. Rvachev, Prof. E. B. Koreneva, and Senior Researcher 
I. Z. Reutfarb. 
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Chapter 2 



The issue of the correct 
formulation of problems of 
mathematical physics 

2.1 Hadamard's definition of correctness 

J. Hadamard has defined two conditions that should be satisfied by a 
correctly formulated boundary-value (initial-boundary-value) problem 
for partial differential equations: existence and uniqueness of the solu- 
tion ( [1], p. 12).^ At the same time, the third condition of Hadamard's 
definition of correctness is well-known that concerns continuous depen- 
dence on the data of the problem. Indeed, he paid serious attention 
to the investigation of this issue with regard to Cauchy-Kovalevskaya's 
theorem concerned with the solution of the differential equation 

d^u = f (t, xi, X2,..., Xn, dtu, d^,u, dro^u, d^^uj (2.1) 

(a system of analogous equations), where / is an analytical function of 
its arguments in the vicinity of the origin of coordinates, with initial 
conditions 

U (0, Xi, X2,..., Xn) = (fio (xi, X2,..., Xn) ] 

^For the first time, tiie concept of correct formulation was put forward by J. 
Hadamard in his article of 1902. 
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(2.2) 



(0, Xi, X2,..., Xn) = (fik-l {Xl, X2,..., Xn) ■ 



As is pointed out by J. Hadamard, the consideration of the problem 
(2.1), (2.2), named after Cauchy, raises three questions ([1], p. 17): 

1) Does it admit a solution? 

2) Is the solution unique? (In general, is the problem correctly 
formulated?) 

3) How the solution can be actually derived? 

Cauchy- Kovalevskaya's theorem (in its authors' interpretation) states 
that, except for some special cases, the above-mentioned problem ad- 
mits a unique solution that is analytical at the origin of coordinates. 
Moreover, the functions (po,..., (pk-i can be not only analytical but reg- 
ular, i.e., continuous together with their derivatives up to a certain 
order. A possibility of a uniform approximation of <fo,---, fk-i by Tay- 
lor series expansions in powers of Xi, Xn , retaining all operations 
on analytical functions (including differentiation up to a corresponding 
order), is implied. 

However, such an approach was strongly criticized by J. Hadamard. 
In his opinion, the question is not how such an approximation affects 
the initial data, but rather what is an effect on the solution? He em- 
phasized the non-equivalence of the notion of small perturbation for 
given Cauchy 's problem and of the solution to this problem ([1], p. 
39). In this regard, J. Hadamard presented his prominent example of 
a solution of the differential equation 



where is a rapidly decreasing function of n. 

The expression on the right-hand side of (2.4) can be arbitrarily 
small. Nevertheless, the problem admits the solution 



d^u + dlu = 0, 



(2.3) 



subject to the conditions 



u (x, 0) = 0; dtu (0, x) — an sin {nx) , 



(2.4) 



u (x, t) 



a. 



■n 



sin {nx) sinh {nt) . 



(2.5) 



n 



2.1. HADAMARD'S DEFINITION OF CORRECTNESS 
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For a.„ = 1/n or or e~^, this solution is rather large for any 

nonzero t, because of the prevailing growth of e"* and, correspondingly, 
of sinh {nt). Thus, the function (2.5) does not depend continuously on 
the initial data and, as a result, the problem (2.3), (2.4) is incorrect. 

Concerning the regularity of the right-hand side of (2.2), J. Hadamard 
remarked: " ...actually, one of the most curious facts of the theory is that 
equations, seemingly very close to each other, behave in a completely 
different way" ([1], p. 29). 

A large number of investigations devoted to the issue of the correct 
formulation of Cauchy's problems. The authors of these investigations 
concerned themselves with specification of corresponding classes of dif- 
ferential equations and with minimization of requirements imposed on 
the initial data (see [2]). However, we are mostly interested in the actual 
character of the dependence of the solution on the data of the problem 
and, in this regard, the classic J. Hadamard's statement that "an ana- 
lytical problem is always correctly formulated in the above-mentioned 
sense, when there exists a mechanical or physical interpretation of the 
question" ([1], p. 38). 

As was pointed out by V. Y. Arsenin and A. N. Tikhonov [3], the 
latter questioned the legitimacy of studies of incorrect problems, speci- 
fied by the authors as the following: the solution of integral equations of 
the first kind; differentiation of approximately known functions; numer- 
ical summation of Fourier series whose coefficients are approximately 
known in the metric I2; analytical continuation of functions; the so- 
lution of inverse problems of gravimetry and of ill-defined systems of 
linear algebraic equations; minimization of functionals for divergent se- 
quences of coordinate elements; some problems of linear programming 
and of optimal control; the design of optimal systems and, in partic- 
ular, the synthesis of aerials. It is emphasized that this list is by no 
means complete, because incorrect problems appear in investigations 
of a broad spectrum of problems of physics and engineering. 

In his talk at the meeting of the Moscow Mathematical Society de- 
voted to J. Hadamard's memory, G. E. Shilov said the following [4]: 
"Our time has brought about corrections in Hadamard's instructions, 
because it turned out that incorrect, according to Hadamard, problems 
could have meaning (as, e.g., the problem of restoration of a poten- 
tial from scattering data). However, the studies of correct problems. 
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proclaimed by Hadamard, was a cementing means for the formation of 
the whole theory" (functional analysis is implied). This quotation is 
borrowed from a biographical sketch by E. M. Polishtuk and T. O. Sha- 
poshnikova [5] , where it is also pointed out that in the course of time 
J. Hadamard's opinion about the importance for practice of exclusively 
correct problems was understood in a less absolute sense. 
At the same time, rather sharp statements were made: 
" And what is more, Hadamard put forward a statement that incor- 
rect problems had no sense at all. Since (as can be seen from a modern 
point of view) most applied problems, represented by equations of the 
first kind, are incorrect, this statement of the outstanding scientist, ap- 
parently, strongly slowed down in 1920-1950's the development of the 
theory, methods and practice of the solution of problems of this class" 
([6], p. 12). 

" Until quite recently, it was thought that incorrect problems had no 
physical sense and that it was unreasonable to solve them. However, 
there are many important applied problems of physics, engineering, 
geology, astronomy, mechanics, etc., whose mathematical description is 
adequate although they are incorrect, which poses an actual problem 
of the development of efficient methods of their solution" ([6], p. 225).^ 

"From the results of this work [of A. N. Tikhonov] followed a limi- 
tation of the well-known notion of J. Hadamard [1] of a correctly for- 
mulated problem of mathematical physics, which was of indisputable 
methodological interest, and inconsistency of Hadamard's thesis, wide- 
spread among investigators, that any incorrectly formulated problem 
of mathematical physics was unphysical." ([7], p. 3). 

"For a long time, activities related to the analysis and solution of 
problems called incorrect used to be relegated (by famous mathemati- 
cians too) to the domain of metaphysics" ([8], p. 126). "A prevailing 
number of mathematicians (including Hadamard) expressed.their atti- 
tude towards this problem in the following way: If a certain problem 
docs not meet the requirements of correctness, it is of no practical 
interest and, hence, does not need to be solved" ([8], p. 127) (I. G. Pre- 
obrazhcnskii, the author of the section "Incorrect problems of mathe- 
matical physics"). 



^In the context of what follows, we draw attention to the "adequate description" . 
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Note that the latter paper most distinctively reveals the style that 
causes a principal objection. Thus, A. Poincare is accused of incon- 
sistency of methodological views on the nature of causal relationship 
([9]) ("The Last Thoughts"). Indeed, the text does not contain any 
evidence that he makes a fetish of the problem of restoration of the 
cause from the effect. On this basis, a conclusion is made about the 
great scientist's misunderstanding of the essence of instability of com- 
putation procedures inherent to incorrect problems and, in particular, 
to integral equations of the first kind.^ 

The adequacy of employed models to considered concrete processes 
is not even touched on. Thus, a quite legitimate question arises: How 
does one know that A. Poincare, if necessary, could not find a way of 
a mathematically correct formulation of the same physical problems? 
Anyway, is there any contradiction in general arguments for the ex- 
istence of such a possibility, including the aspects of its constructive 
realization? 

By the way, exactly A. Poincare repeatedly mention J. Hadamard 
while establishing a relationship between the correct formulation of 
problems and a practical realization of employed models. We draw 
attention to an expressive thesis: "If a physical problem reduces to an 
analytical one, such as (2.3), (2.4), it will seem to us that it is governed 
by a pure occasion (according to Poincare, it means that determinism 
is violated) and it does not obey any law" ([1], p. 43). 

In light of the above, the arguments of I. Prigogine and I. Stengers 
[10] are of interest: "...one can speak of a 'physical law' of some phe- 
nomenon only in the case when this phenomenon is 'coarse' with respect 
to a limiting transition from a description with a finite accuracy to that 
with an infinite accuracy and thus inaccessible to any observer, whoever 

^In particular, the exact statement reads: "However, one must remember that 
vagueness of philosophical positions of some scientists in the West, even rather 
renown, results in the fact that, based on correct starting points, they draw rather 
inaccurate conclusions, repeating old mistakes of, for example, A. Poincare, who 
writes: 'If two organisms are identical, or simply similar, this similarity could not 
occur by chance, and we can assert that they lived under the same conditions...' In 
other words, the fact of possible incorrectness of the inverse problem is completely 
ignored." However, one would hardly mention Poincare's mistakes... if modern 
" spontaneous supporters of the principle of determinism did not repeat them" ( [8] , 
p. 134). 
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he may be" (p. 9). "Scientist in a hundred different ways expressed 
their astonishment that a correct formulation of the question allows 
them to solve any puzzle suggested by nature" (p. 44). 

Thus, underlined are: first, methodological importance of correct 
formulation of problems of mathematical physics; second, a leading 
role of the employed procedures and, finally, substantial influence of 
the quality of their realization on the degree of complexity of obtain- 
ing the final result. In other words, one and the same problem can be 
better or worse formulated. The above-mentioned statement, which is, 
essentially, the same hypothesis of J. Hadamard, as a matter of fact, 
implies a possibility of a "good" (correct) formulation of any meaning- 
ful problem and, consequently, can be interpreted as having a global 
orientation. 

In this regard, one can establish an obvious relationship to D. 
Hubert's comments on his 20th problems that suggested a possibility 
of correct formulation of arbitrary boundary-value problems of mathe- 
matical physics by means of special requirements on boundary values 
of corresponding functions (a type of continuity or piecewise differen- 
tiability up to a certain order) and, by necessity, by giving an extended 
interpretation to the notion of the solution ([11], pp. 54-55). 

For the first time, the three conditions of the correctness of problems 
of mathematical physics were clearly pointed out by D. Hilbert and R. 
Courant ( [12], pp. 199-200): existence, uniqueness and continuous 
dependence of the solution on the data of the problem. Concerning 
the last, they say: "...it has crucial importance and is by no means 
trivial... A mathematical problem can be considered adequate to the 
description of real phenomena only in the case when a change of given 
data in sufficiently narrow limits is matched by an alike small, i.e. 
restricted by predetermined limits, change of the solution" . 

V. A. Steklov's position is quite analogous ([13], p. 62): "...if differ- 
ential equations with the above-mentioned initial and boundary condi- 
tions are not constructed on erroneous grounds, are not in direct con- 
tradiction to the reality, they must yield for each problem a unique and 
completely definite solution...". Along the same fines, I. G. Petrovskii 
writes ([14], p. 87): 'The above-mentioned arguments for the correct 
formulation of Cauchy's problem show that other boundary- value prob- 
lems for partial differential equations are of interest for natural science 
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only in the case when there is, in a sense, continuous dependence of the 
solution on boundary conditions" . 

S. L. Sobolev is less categorical ([15], p. 38): "The solution to an 
incorrectly formulated problem in most cases has no practical value" . 

Of considerable interest is the opinion of V. S, Vladimirov ([16], p. 
69): "The issue of finding correct formulations of problems of mathe- 
matical physics and methods of their solution (exact or approximate) is 
the main content of the subject of equations of mathematical physics". 

V. V. Novozhilov, in fact, drew attention to the potential of vari- 
ation of the formulation of the considered problem with the aim of 
the simplification of the procedure of its numerical realization ([17], p. 
352): "The absence in the term "a mathematical model" of the indi- 
cation of its inevitable approximate character leaves way for a formal 
mathematical approach to models, disregarding those concrete prob- 
lems for whose solution they were intended, which is, unfortunately, 
wide-spread at present" . 



2.2 J. Hadamard's hypothesis and incor- 
rectness of "real" problems 

Thus, J. Hadamard and other outstanding scientists thought that any 
physically intcrprctable problem could be correctly formulated. How- 
ever, a quite opposite point of view dominates in modern publications. 
Indeed, a visibly larger part of practically important problems con- 
sidered therein are incorrect. However, is the actual methodology of 
mathematical formulation of these problems and, correspondingly, the 
results of its refraction with respect to realities adequate? 

Here we will not elaborate on something like general principles of 
the construction of differential equations, and, generally speaking, it is 
reasonable at the beginning to restrict the question to the following: 
What arguments allow one to conclude that an incorrectly formulated 
problem adequately describes an observable phenomenon or a poten- 
tially real process? In this regard, let us turn to the procedure of the 
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solution of the Fredholm integral equation of the first kind 
1 

jk{x,i),l^{i)di = f{x), a:e[0,l], (2.6) 



which is a classical incorrect problem: the kernel /c(a;, ^)and the free 
term / {x) are given; the function -0 [x) is to be restored. 

For definitiveness, let the kernel k be symmetric and closed, i.e. 
k{x,^) = k{^,x) and its eigenfunctions 'ipni^), being nontrivial solu- 
tions of the integral equation 

1 

^l^{x)^\jk (x, ^ (0 di = f{x), xe [0, 1] 



with characteristic numbers A = A„, n = 1,2, . . ., form a complete in 
L2 (0, l)orthogonal system of elements. Besides, the function / (x) e 
L2 (0, 1). In this case, according to Picard's theorem, the solution to 
Eq. (2.6) exists and is unique under the condition (see, e.g., [18]) 

00 J. 

alXl < 00, = / / (x) VS„ (x) dx. (2.7) 

n=l ^ 

If we assume that all the above-mentioned conditions are fulfilled, 
there is still the third condition of correctness that, as is known ([19], 
p. 69), is not satisfied by Eq. (2.6). Numerous literature references 
clearly illustrate an inadequately strong influence on the solution of 
small perturbations of the data of the problem, in the first place of / (x) . 
As a rule, this function is determined experimentally and mismatch the 
kernel k{x,^) with respect to smoothness. Thus, Eq. (2.6), strictly 
speaking, looses sense. At the same time, a possibility of an equivalent 
description of realistic problems by means of integral equations of the 
first kind is indisputably admitted at present, which is confirmed by 
their colossal list [6]. 

Let us specify Eq. (2.6): 

/ (1-^)^, 0<e<,T; . 1 . , , 
k{x,t) = < ,^ ^^^-1 t (x) — T^smimnx), 

(2.8) 
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where m is an integer. 

For this choice, ipn i^) = sin {nnx); A„ = (nvr)^, n = 1, 2, . . . . The 
use of Mercer's theorem [18], according to which 

n=l 

and a representation of ip {x) as a series expansion in terms of ■0„ {x) 
with undetermined coefficients allows one to find the solution to Eq. 
(2.6): 

ip (x) = sin {mirx) . (2-9) 

However, the procedure of calculations turned out to be so simple 
owing to a special choice of the data of the problem. If this is not the 
case or in the case of the solution of Eq. (2.6) with the kernel and the 
free term (2.8) by means of one of numerical methods, the complexity 
of the realization of a an approximation of sufficiently high order is 
practically identical to the most general situation, characterized by an 
error in the determination of / (x) and k (x,^).'^ As a matter of fact, 
even if the data are objectively compatible, the incorrectness of Eq. 
(2.6) is incurred by an error in the approximation of / and k and also 
by rounding off the digits by the computer. 

The key factor of the incorrectness of Eq. (2.6) follows from a 
comparison of the free term (2.8) with the solution (2.9). Indeed, by 
increasing m, the function / (x) can be made arbitrarily small without 
any change of the amplitude value of ip (x) . Correspondingly, any error 
in the calculations with / {x) is projected onto {x) with the factor 
m^. The mechanism of this effect related to the smoothing of informa- 
tion about the sought function by the procedure of integration will be 
repeatedly discussed in what follows. 

However, let us return to the question of the relation of an incorrect 
formulation to the reality. In this regard, we draw attention to the 
following. By considering (2.6) as the Fredholm integral equation of the 
first kind (the function ip is to be determined), we mean the solution 
of the inverse problem (I). However, equation (2.6) can be used for the 

*This complexity implies an ill definition of systems of linear algebraic equations 
obtained as a result of some sort of discretization. The methods allowing one to 
overcome this difficulty are discussed in section 3. 
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solution of the corresponding direct problem (D): the determination of 
the function / (x) from the data k (x,^) and ip (x). This procedure is 
correct and thus is radically simple than the problem I. It is sufficient 
to note the absence of any principal difference between an essentially 
numerical realization and the integration of expressions (2.8), (2.9) in 
the analytical form. 

Here we want to draw attention to an issue that seems to be of 
substantial importance. The problem D, as a rule, is transparent: in its 
categories, we quite naturally model current processes and phenomena 
by, which should be emphasized, explicit means of linear superposition. 
Correspondingly, if, for instance, k {x, ^) is a characteristic of the media 
and ip (x) is intensity of external influence, a resulting effect in this or 
that subject sphere is to be elementarily summed up. 

The situation is diametrically different for the problem I. One could 
hardly point out any realistic process (phenomenon) for which it could 
be formulated in mathematical terms directly on the basis of the subject 
sphere. In other words, without any relation to the problem D, which 
commonly implies a transformation of the latter into the problem I by 
means of mechanical renaming of known and unknown components. 

In our opinion, the methodology, which states the adequacy of the 
problem I, obtained by the above-mentioned renaming of the compo- 
nents, to the realities on the basis of a quality information about a 
concrete problem D, is profoundly erroneous. Correspondingly, the 
statements of experts who reject J. Hadamard's hypothesis of the ex- 
istence of correct formulations of problems of mathematical physics 
should be considered unjustified. 

Let us turn to the problem D that describes some realistic process 
(2. 6). For this process, the determination of ijj (x) from the data k {x, ^) 
and f{x), i.e. the formulation of the corresponding inverse problem 
that will be denoted as I , is, of course, reasonable. Suppose that in 
this case J. Hadamard's hypothesis is correct, and, hence, the problem 
I is correct. However, the problem I', the solution of the Frcdholm 
integral equation of the first kind (2.6), is incorrect by definition. 

The conclusion is obvious: Mathematical formulations of the prob- 
lems I and I are non-identical. As a result, the formulation of the 
problem I cannot be restricted to readdressing the status of the un- 
known variable between the functions / and ip in the problem D. Note 
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in this regard that constructive methodology of the correct formulation 
of the problem I is the basis of the subject and the main objective of 
the present investigation. 

The above arguments seem to be rather convincing, however, at 
this stage of our consideration, we can neither prove the correctness 
of J. Hadamard's hypothesis nor illustrate its constructive character 
with respect to separate classes of problems. One should also bear in 
mind that, using special methods, the solution of the incorrect problem 
I, as a rule, can be obtained with accuracy that is considered to be 
practically acceptable. In this regard, the question arises: Should one 
aim at the correct formulation I , if the algorithm of the calculation 
of the function ijj (x) in the formulation of the problem I in some way 
realizes its regularization? This implies a well-known deformation of 
the formulation I with the use of a small parameter that yields the 
property of correct solvability. 

Thus, can the algorithm to a full extent, including the efficiency of 
numerical realization, level off the principal difficulties inherent to the 
incorrectness of the problem I in the form (2.6)? It is clear that the 
answer is definitely negative: Otherwise, the deep-rooted differentiation 
of problems in correctly and incorrectly formulated ones would make 
no sense. 

Furthermore, the indicated difference is of exceptional importance, 
because correctness of the formulation is a criterion of a qualitative 
level, whereas the efficiency of a method of the solution of the Fred- 
holm integral equation of the first kind can be estimated only in terms 
of quantitative factors of a palliative property. The latter is caused by 
a direct relationship between a degree of regularization and the defor- 
mation (distortion) of the problem I. 

What is, however, the actual difference in the interpretation of the 
formulations I and I ? The answer to this question is the basic concept 
of the suggested approach. At this stage, we only note that a transfor- 
mation of the formulation I into the formulation I' will be mostly qual- 
itative and will be realized by means of zero in L2 (0, 1) perturbation 
of Eq. (2.6), modeling the phenomenon of smoothing of information. 
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2.3 Banach's theorem on the inverse op- 
erator 

Let us quote ([5], p. 175): "First, Hadamard defined the correctness of 
tlie problem by the conditions of solvability and uniqueness and strongly 
insisted on continuous dependence of the solution on the initial data 
only in the consideration of Cauchy's problem. In the book "The theory 
of partial differential equations", published in Peking a year after his 
death, he wrote: 'This third condition that we introduced in "Lectures 
on Cauchy's problem..." but did not consider as part of well formulated 
problems, was added, quite justified, by Hilbert and Courant [12]. Here, 
we accept their point of view." ' 

E. M. Polishuk and T. O. Shaposhnikova made the following com- 
ment on this text [5], pp. 175-176]: "From a mathematical point of 
view, the question of the necessity of the requirement of the continuity 
of the solution with respect to the data seems to be rather delicate. 
As a matter of fact, according Banach's well-know theorem on closed 
graph, unique solvability of a linear problem leads to boundedness of 
the inverse operator and, thus, continuous dependence of the solution 
on the right-hand sides." It is pointed out that variations of the coef- 
ficients of differential equations and of the boundary of the considered 
domain can also influence the solution of the problem; hence, the use 
of the three conditions of the correctness is preferable. 

At the same time, Banach's theorem on the inverse operator ([20], 
pp. 34-36), being a consequence of the above-mentioned one, is more 
closely related to the considered issue. Its formulation, given by A. 
L Kolmogorov and S. V. Fomin, is the following ([21], pp. 259-260): 
Let A be a linear bounded operator that maps a Banach space Bi in a 
one-to-one fashion onto a Banach space B2. Then the inverse operator 
A~^ is unique. 

In addition, L. A. Lyusternik and V. 1. Sobolcv ([22], pp. 159- 
161) emphasized that a one-to-one mapping of the whole Banach space 
Bi onto the whole Banach space B2 is implied. Besides, a situation is 
discussed when " ...an operator, being the inverse of a bounded operator, 
although linear, turn out to be defined not on the whole space B2 but 
only on a certain linear manifold and unbounded on this manifold" . 
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A somewhat different interpretation of tfie same theorem in ([23], 
p. 60) is reduced to the following: If a linear bounded operator A that 
maps a Banach space Bi onto a Banach space B2 has an inverse A~^, 
then A~^ is bounded. It is pointed out that this statement becomes 
invalid if one gives up the requirement of completeness of one of the 
spaces. There is also a clarification: The existence and uniqueness of 
the solution of the equation Aip = f with an arbitrary right-hand side 
from B2 leads to continuous dependence of the solution t/j — A"^/ on 
/■ 

L. V. Kantorovich and G. P. Akilov made an addition concerning a 
mapping under the specified conditions onto a closed subspace of the 
Banach space B2 ([24], p. 454). As a matter of fact, a closed subspace 
of a Banach space is itself a Banach space. 

S. G. Mikhlin gave a proof of the theorem ([25], p. 507): For the 
linear problem Aip = / to be correct in a pair of Banach spaces Bi, B2, 
it is necessary and sufficient that the operator A^^ exist, be bounded 
and map the whole space B2 onto Bi . At the same time, the author 
made a clear division between the category of the existence and unique- 
ness of the solution of the boundary- value problem and its correctness 
as a whole, which implies, as a result, continuous dependence on the 
data (the third condition according to Hadamard). In this regard, note 
a typical definition: "A boundary- value problem is called correct in a 
pair of Banach spaces Bi, B2 if its solution is unique in Bi and exists 
for any data from B2, and if an arbitrarily small change of the solution 
in the norm Bi corresponds to a sufficiently small change of the initial 
data in the norm B2" (p. 204). 

It is emphasized that the considered problem may turn out to be 
correct in one pair of spaces and incorrect in another one. Besides, 
the incorrectness of the Predholm integral equation of the first kind 
(2.6) is established by contradiction: If the problem is correct, there 
exists a bounded operator A~^ and, hence, the identical operator / = 
A~^A is compact in an infinite dimensional space, which contradicts 
the fundamentals of the general theory [23]. S. G. Mikhlin also quite 
encouragingly pointed out the approach of an approximate solution of 
incorrect problems, headed by A. N. Tikhonov. 

In an analogous, as to its content, course [26], S. G. Mikhlin reit- 
erated the above-mentioned formulations. However, A. N. Tikhonov is 
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not mentioned at all, whereas the discussion of Eq. (2.6) found a rather 
interesting continuation, (p. 171). It is shown that the problem of its 
solution becomes correct a pair of spaces Bi , B2 is replaced with such 
one that the operator A is no longer compact. The general consider- 
ations are illustrated by the following example. Let k {x, ^) and / (x) 
satisfy the conditions of section 2.2, including (2.7). Then, according 
to Picard's theorem, the solution to Eq. (2.6) has the form 

00 

1p{x)^Yl ('nXni^n (x) , [0, 1] . (2.10) 

n=l 

It turns out that if one retains L2 (0, 1) as B2 and for Bi also takes 
a Hilbert space of functions normalized according to (2.7), i.e. I2, the 
solution of Eq. (2.6) becomes a correctly formulated problem, and, cor- 
respondingly, the operator A~^ is bounded. Indeed, the value domain 
of the operator A in I2 is not closed ([27], p. 34). 

A decade later, S. G. Mikhlin, in fact, gave up the investigations 
related to the issue of correctness [28]: "The author adheres to the 
classical point of view, according to which the problem being solved by 
mathematical methods should be considered as correctly formulated. 
Of course, there are other opinions (p. 7)... Thus, we neglect the 
so-called incorrigible errors related to the formulation of the above- 
mentioned problem as a problem of natural science or of social studies 
(measurement errors, insufficient accuracy of basic hypotheses, etc.)" 
(p. 17). 

M. M. Lavrentiev and L. Y. Savehev characterized investigations of 
the issue of the solvability of Eq. (2.6) on the basis of considerations 
of the type of [26] as trivial, because it is difficult to imagine that for 
experimentally determined / (x) the corresponding error may prove to 
be small in the norm of the space I2 ([29], p. 217). At the same time, it 
is pointed out that, generally speaking, for any operator equation, one 
can choose pairs of spaces such that the problem of its solution will be 
correct. 

G. M. Vainikko and A. Y. Veretennikov draw attention to the com- 
plexity of the description of such spaces. And what is more, even the 
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Volterra integral equation of the first kind 

X 

jk{x,i)i;{i)di^f{x), xe[o,i], 



which admits the regularization 

X 



and is elementarily solvable by quadratures, for reasons of the norm 
for ( rule, has to be considered as an incorrectly formulated 

problem ([30], p. 6). 

As regards the pair of spaces that realize the conditions of the cor- 
rect formulation, an original remark of K. I. Babenko is of interest ([31], 
p. 304): "Hadamard's well-known example (2.3), (2.4) that yields the 
solution of Cauchy's problem of the type (2.5) by no means tells of the 
absence of continuous dependence on the initial data, as it is usually 
interpreted. It rather tells of the fact that small changes of the initial 
data may result in leaving the totality of the initial data for which the 
solution of Cauchy's problem exists." By the way, R. Richtmyer demon- 
strated the correctness of the procedure of a numerical realization of 
a rather complicated problem of the above-mentioned type with the 
representation of sought functions by two-dimensional power series and 
with the use of special methods of suppression of errors of arithmetical 
operations ([32], section 17. B). 

In the context of the present consideration, two theorems, given by 
V. A. Trenogin, are of interest ([33], p. 225): 

Let El and E2 be infinite dimensional normalized spaces, with E2 
being complete. If yl is a compact linear operator from Ei to E2, 
different from finite dimensional, its value domain R {A) is not a closed 
manifold in E2. 

Let A be a compact operator from an infinite dimensional normal- 
ized space El to a normalized space E2, with the inverse operator A~^ 
existing on R {A). Then A~^ is bounded on R (A). 
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2.4 The premises of the reahzation of the 
conditions of correctness 

Let us assume that f = f* (x) is the result of exact integration of the 
function ip (x) e L2 (0, 1) and of the symmetric kernel k {x, ^) by means 
of the formula (2.6). In this case, /* e ^2- Moreover, the corresponding 
set, which is important, is not compact. This leads to the boundedness 
of the operator A~^ from I2 to L2 (0, 1) defining the solution of the 
equation 

{Ai;){x) = f4x), xe[0,l], (2.11) 

as mentioned above. There is no need to explain that such a situa- 
tion would be, in principle, rather favorable, in view of the fact that 
Hadamard's third condition of correctness, directly following from the 
existence and uniqueness of the solution, is fulfilled automatically. 

The reason is the stability, by definition, of the inverse procedure of 
the restoration of ip {x) with respect to small, in a sense, perturbations 
of the data of the problem. Correspondingly, the solution to Eq. (2.11) 
can be obtained by the use of the formula (2.10), and the error of the 
evaluation of the coefficients a„, following (2.7) with / = /* (x), and of 
the values of A„ will adequately influence ip (x) in L2 (0, 1). 

However, in reality, the coefficients a„ are determined with an error 
equivalent to {5f) (x), which is inherent to the definition of the free 
term in (2.6). That is, in reality, the right-hand side of Eq. (2.11) has 
the form 

fix) = Mx) + iSf)ix). (2.12) 

As a matter of fact, a necessary condition of the correctness of Eq. 
(2.11), namely / G I2, proves to be illusive, because it is impossible to 
obtain any quantitative information about {5f) (x). 

Note that specific" diffusion" of the space I2 is caused by the repre- 
sentation of the normalizing functional (2.7) as an infinite series, whose 
terms are determined on the basis of spectral characteristics of the 
kernel k{x,^). In this sense, such a space as L2 (0, 1) is much more 
explicit for the function / {x) . However, the use of it incurs rather 
negative consequences. 

Indeed, in this case R (A) does not belong to the closed space 
L2 (0, 1), the operators A and A~^ become, respectively, compact and 
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unbounded. As a consequence, an investigation and a numerical real- 
ization of the original equation (2.6), in fact, turn out to be beyond the 
sphere of the application of Banach's fundamental theorem on the in- 
verse operator. Isn't it a too high price to pay for seemingly ephemeral 
clarity in the formulation of the problem? 

We draw attention to a known point of view that a choice of ap- 
propriate spaces for the solutions to problems of mathematical physics 
should be done on the basis of practical applications, which can hardly 
be disputed. As the same time, a wide-spread opinion that, for example, 
a sociologist should formulate a problem to be solved by mathematical 
methods with a specification of appropriate spaces for its data. This, 
as a rule, admits variety, which is a prerequisite for an increase in the 
efficiency of procedures of numerical realization. 

Are there any prospects to overcome the above-mentioned complex- 
ity in mating the free term of the Predholm integral equation of the first 
kind (2.6) with the adequate space ^2? In this regard, let us turn to 
Eq. (2.11) that by virtue of (2.12) takes the form 

(AV) (x) = / (x) - (5f) (x), xe [0, 1] . (2.13) 

It seems that there is a chance of a reduction of the given function 
/ (x) to G I2 by means of adaptive simulations of the error (Sf ) (x). 
Indeed, Sf and, essentially, the smoothing of information by the proce- 
dure of integration can be naturally reflected by the difference between 
the explicit form of the sought function ip (x) and an integral over this 
function whose kernel would not impose any additional restrictions on 
the formulation of the problem.^ As a development of analogous con- 
siderations, there appears the condition 

P/IIl.(o,i) = 0. (2.14) 

Thus, instead of traditional restoration of the function ip (x) directly 
from the solution of the Predholm equation of the first kind, we suggest 
to employ a nonessential, at the first sight, perturbation (zero in L2), 
which is adequate to the problem (2.13), (2.14). As will be shown below, 

realization of this approach is a key aspect of the constructive part of the 
present consideration (see section 5). 
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as a result of comparatively simple transformations, the determination 
of ip (x) will be reduced to the solution of the Fredholm integral equation 
of the second kind, which is rather preferable from the point of view of 
the efficiency of procedures of a numerical realization. 

In conclusion, we want to point out that known methods of the 
solution of equations of the type (2.6) are widely used when there 
is a cardinal mismatch between R{A) in L2 (0, 1), or in some other 
space, and the functions / (x) that are determined by means of mea- 
surements. Prom this point of view, the fulfillment of the condition 
(2.14) can hardly be interpreted in a direct way. Nevertheless, the 
above-indicated approach still holds by interpreting, figuratively, an 
approximation of the considered problem to the correctly formulated 
one. Thus, the convergence of the series (2.7) must be implied with 
regard to a corresponding variation of the kernel k{x,^). As far as 
algorithms of the solution of resulting Fredholm integral equations of 
the second kind are concerned, their stability is practically not affected 
by the error in the determination of / (x) e L2 (0, 1) (see, e.g., [19]). 
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Chapter 3 

The existing approaches to 
the solution of incorrect 
problems 

3.1 A. N. Tikhonov's methodology 



The consideration of this subsection is based on the material of the 
monograph by A. N. Tikhonov and V. Y. Arsenin [1] that is, hterally, 
pieced by the concept of the adequacy of incorrect formulations and, in 
particular, of integral equations of the first kind to problems of math- 
ematical physics. As an illustration, we show that the solution to the 
Predholm integral equation of the first kind 



{A^) {x) = jk (x, ^ (0 di-fix), xe [a, b] , (3.1) 



with k{x,C,) and dxk{x,C,) being continuous with respect to 
undergo arbitrarily considerable changes both in the metric C and L2 
for small in L2 (a, b) variations of the right-hand side in the form 



The situation with the perturbation of the kernel /c (x, ^) is, in fact, 
analogous. In this regard, the authors pose the question: What should 



b 



a 



b 




a 
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be understood as the solution of Eq. (3.1) when k and / are known 
approximately? In their opinion, a problem of this type should be con- 
sidered as " underspecified" , and, correspondingly, a choice of possible 
solutions should be made taking into account " usually" available addi- 
tional qualitative or quantitative information about the function i/j (x). 
In this regard, we draw attention to N. G. Preobrazhenskii's consider- 
ations concerning a system of linear algebraic equations, obtained by 
the discretization of (3.1) ([2], p. 130): 

"An analysis shows that choosing sufficiently high order of an ap- 
proximation, we transform [the above-mentioned problem] into an ar- 
bitrarily ill-defined one... Under these conditions, it is necessary to add 
to the algorithm some a priori nontrivial information, only by the use 
of which we can expect to filter out veiling false variants and single 
out the solution, closest to the sought one. any purely mathematical 
tricks that do not employ additional a priori data are equivalent to 
an attempt to construct an informational perpetuum mobile producing 
information from nothing." 

The so-called method of the selection of the solution to incorrectly 
formulated problems is based on a priori quantitative information. It 
is shown that if a compactum M of a metric space Ei is mapped in 
a one-to-one and continuous manner onto a set F of a metric space 
E2, the inverse map F onto M is also continuous. Correspondingly, 
an assumption that the solution, in particular, to Eq. (3.1) belongs 
to the compactum M allows us to consider the operator A~^ to be 
continuous on the set F = AM. A practical realization is reduced to 
an approximation of M by a series with parameters that change within 
certain limits (for M to represent a closed set of a finite dimensional 
space) and should be determined from the condition of the minimum 
of the error of closure of (3.1). Note the absence of any more or less 
general recommendation with respect to the choice of M. 

In light of the above, M. M. Lavrentiev has formulated the notion 
of correctness according to Tikhonov for an equation of the type (3.1), 
with the functions ip and / belonging to Banach spaces Bi and B2, 
respectively [3]: 

1) It is a priori known that the solution to the considered equation 

exists and belongs to a set M of the space Bi. 

2) The solution is unique on the set M. 
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3) The operator is continuous on the set AM of the space B2. 

If M is a compactum (this case is called " usual" ) the last condition 
becomes a consequence of the first two conditions. 

Those problems in which the operator A"^ is unbounded on the set 
AEi and the set of possible solutions Ei is not a compactum are called 
substantially incorrect. For such problems, A. N. Tikhonov has put 
forward an idea of a regularizing operator G, in a sense close to A~^, 
whose value domain for the map from E2 into Ei admits matching to 
the right-hand side of (3.1), known approximately. Moreover, G must 
contain a regularization parametera that depends on the accuracy of 
the initial information. 

The operator G (/, a) is called a regularizing operator for Eq. (3.1) 
if it possesses the following properties: 

1) It is defined for any a > and f E E2. 

2) For Aip^ — where ip^ and /* are corresponding exact expres- 
sions, there exists such a (S) that for any < e < pe^ {i^^ii^a) there is 
5 (e) > pE2 ii',,^a)- Here, = G{f,a). 

It is implied that there is a possibility of a choice of a (6) such that 
for 5 — > the regularized solution ■j/'^ ^ i-e., e — > 0. At the same 
time, it is pointed out that the construction of the dependence a (5), 
for which the operator G (/, a (5)) is a regularizing one, is algorithmi- 
cally complicated for classes of practically important problems. There 
are a lot of publications of A. N. Tikhonov's followers devoted to the 
resolution of this difficulties, which will be discussed below. 

In [1] , the construction of G (/, a) is carried out by the use of tech- 
niques of calculus of variations that reduce the evaluation of ijj (x) to 
the minimization of the functional 



For Eq. (3.2), its stabilizing component is recommended to be taken 
in the form 



(3.2) 



b 



^W\ = j {po {x) ijj^ (x) + Pi (x) [ijj' {x)f^ dx, 



(3.3) 



a 



where Po, Pi > are given functions. 
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In the case of a symmetric kernel k {x, the procedure of the mini- 
mization is equivalent to the solution of the integrodifferential equation 



a {po (x) {x) - [pi {x) i^' {x)^}HM) {x) ^f(x), xe [a, b] , (3.4) 



Here, v (x) is an arbitrary variation of ip (x) in the class of admissible 
functions. 

In the opinion of the authors of [4], an overwhelming majority of in- 
verse problems are incorrectly formulated, and attempts to solve them, 
in view of their great practical importance, were being undertaken for 
a long period. "But only as a result... of the appearance of fundamen- 
tal publications of academician A. N. Tikhonov, the modern theory of 
the solution of inverse problems, based on the notion of a regulariz- 
ing algorithm, was constructed" (p. 7). In what follows, the authors 
demonstrate the efficiency of the procedure of a numerical realization 
of the Predholm integral equations of the first kind, related to the inter- 
pretation of astrophysical observations, by means of the selection of the 
compactum of possible solutions in the class of monotonically bounded 
functions. 

As is pointed out by O. A. Liskovets [5], "...the correctness accord- 
ing to Tikhonov is achieved at the expense of the reduction of the 
admissible manifold of solutions to the class of correctness" (p. 13). 
The following quotation from the above-mentioned monograph is also 
of considerable interest: "In contrast to a previously prevailing opinion 
that all the problems describing physical reality are correct, according 
to the modern point of view any realistic problem can be regularized, 
i.e., it has at least one rcgularizer" (p. 14). 

Let us refer to V. A. Morozov's conclusion ([6], p. 9): "A. N. 
Tikhonov's method of rcgularization turned out to be simple in prac- 
tice, because it did not require actual knowledge of the compactum M 
that contained the sought solution to Eq. (3.1)... The main difficulty 
of the application of this method consists in the formulation of algo- 
rithmic principles of the selection of the parameter of rcgularization a" . 
According to his own monograph ([7], p. 4), "the importance of A. N. 



under the conditions 




(3.5) 
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Tikhonov's paper [8] can hardly be overestimated. It served as impetus 
for a number of publications by other investigators in different fields of 
mathematical analysis and natural science: spectroscopy, electron mi- 
croscopy, identification and automatic regulation, gravimetry, optics, 
nuclear physics, plasma physics, meteorology, automation of scientific 
research and some other spheres of science and engineering" . 

V. V. Voevodin's considerations [9] are rather typical: " The success 
of the application of the regularization method to the solutions of unsta- 
ble systems of algebraic equations is explained to a large extent by the 
fact that A. N. Tikhonov and his followers did not restrict themselves to 
an investigation of separate fragments of this complicated problem but 
considered the whole complex related issues. This, in the first place, 
concerns a clear formulation of the problem itself, the construction of 
a stable with respect to perturbation of the input data algorithm of its 
solution, the development of an efficient numerical method, estimates of 
a deviation of the actually evaluated object from the sought one taking 
into account a perturbation of the input data and errors of rounding" . 

A quotation from the preface to the collected volume by A. A. 
Samarsky and A. G. Sveshnikov [9] reads: "A clarification of Audrey 
Nikolaevich Tikhonov of the role of incorrect problems in classic math- 
ematics and its applications (inverse problems) is of fundamental im- 
portance for the who;e modern mathematics. He proposed a principally 
new approach to this class of problems and developed methods of the 
construction of their stable solutions based on the principle of regular- 
ization" . 



3.2 A brief review of the development of 
the outhned concepts 

The results of investigations devoted to the determination of the regu- 
larization parameter a are summarized in [10]. Based on the assump- 
tion that errors in the determination of the free term / (x) and the 
kernel k {x, ^) of Eq. (3.1) are known, one uses different methods of the 
minimization of the error of closure of the type 



E2 



E2 ' 



/^e (0,1). 
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The evaluation of the parameter a as a root of the corresponding equa- 
tion does not pose any problem. However, a choice of is, in fact, 
related to considerable uncertainty. The main obstacle is that a reli- 
able estimate of the error caused by incompatibility of the concretely 
considered equation Aip ~ f is rather questionable. 

Considerable efforts were undertaken to reduce the volume of infor- 
mation necessary for the evaluation of the parameter a. A noticeable 
step in this direction was made by A. N. Tikhonov and V. B. Glasko who 
suggested a criterion of the minimization of the functional ||Q;d'0a/dQ;|| 
with respect to a > [11] (sec also [1], section 2.7). However, its theo- 
retical justification proved to be possible only for rather narrow classes 
of problems. A number of methods of the determination of a is related 
to the use of solutions to Eq. (3.1) for a special form of the expressions 

In [10], the status of the studies of estimates of the accuracy of 
methods of the solution of the integral equation (3.1) is also discussed. 
If ip (x) belongs to a compactum, any serious complications, as a rule, 
do not arise, and the main interest is focused on the algorithm of reg- 
ularization. If pi = in (3.3) and the parameter a is finite, Eq. (3.4) 
becomes a Predholm integral equation of the second kind, to which, 
under the assumption that the error in the determination of k (x, ^) 
and / (x) is known, the whole general theory of approximate methods 
of L. V. Kantorovich applies ([12], section 14.1). 

At the same time, as shown by V. A.Vinokurov [13], when a priori 
information about the solution to Eq. (3.1) is missing, the estimate 
of the error of the evaluation of ip (x) by means of regularization is 
impossible in principle. Justified is only a formulation of the question 
of the convergence of the procedure of computation or of a possibility 
of the regularization of the corresponding problem. 

In this regard, we note the arguments of A. B. Bakushinskii and A. 
V. Goncharskii ([14], p. 13): "Unfortunately, in the general case, it is 
impossible to estimate the measure of closeness of G (/, a) to A^^ (/*) 
without additional information about the solution to Eq. (3.1). This 
is a characteristic feature of incorrect problems. In the general case, 
a regularization algorithm ensures only asymptotic convergence of an 
approximate solution to the exact one for 6—^0". 

The name of M. M. Lavrentiev is associated, in fact, with a particu- 
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lar case of a practical realization of A. N. Tikhonov's method consisting 
in the reduction of the problem (3.4), (3.5) to the solution of the Fred- 
holm integral equation of the second kind 

h 

aijj {x)+Jk {x, ^ (0 d^ = f{x), xe [a, h] , (3.6) 

a 

where a > is a small parameter. 

It is shown that ||V'a — — > for 5 — > 0, 7 — > and (5 + 7) /a {6, 7) 
0. Here, 7 is an error in the determination of the kernel k{x,$,), anal- 
ogous to 6 (see section 3.1). 

V. K. Ivanov's method [15] allows one to find the so-called quasiso- 
lution minimizing the error of closure of (3.1) for a class of functions 
■0 (x) e Mr, where Mr e is a compactum whose specification, gen- 
erally speaking, are analogous to that employed in Picard's theorem 
(see section 2.2). The quasisolution to (3.1) on such a compactum has 
the form 

00 

'4'{x) = Y1 (A + \n) {x) , x^ [a, 6] . (3.7) 

n=l 



Here 



/ {x) {x) dx; 



Xn and 'ipn (x) are, respectively, characteristic numbers and eigenfunc- 
tions of the kernel k{x,^); the parameter A = and represents a posi- 
tive root of the equation 

under the conditions, respectively, 

oo oo 

E«<^'; E«>^'- (3-9) 

n=l n=l 

Special methods of regularization are developed for the situation 
when there is considerable volume of information of statistical char- 
acter (spectral densities, mathematical expectations, etc.) about the 
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solution to an equation of the type (3.1). Thus, V. N. Vapnik [16] rather 
constructively employed the specifics of problems concerned with recog- 
nition of images, related to nonuniqueness and, as a result, to extreme 
behavior of the sought functions. We point out a definition in the 
above-mentioned monograph (p. 8) that, apparently, was implied by 
many authors but did not receive such a clear formulation. 

" The problem of the restoration of dependencies from empirical data 
was and, probably, will always be central in applied analysis. This 
problem is nothing but mathematical interpretation of one of the main 
problems of natural science: How to find the existent regularity from 
random facts." 



3.3 V. M. Pridman's approach 

Let k{x,^) be symmetric, positive definite kernel and Eq. (3.1) be 
solvable. Then, as shown by V. M. Pridman [17], a sequence of functions 
determined by iteration 



f{x)+jk (X, V'n (0 di 



n = 0, 1,..., 



(3.10) 

converges in L2 (a, h) to the solution of Eq. (3.1) for an arbitrary choice 
of the initial approximation ^iIjq (,t) G L2 (a, h) and < A < 2Ai, where 
Ai is the smallest characteristic number of the kernel k (x,^). 

M. A. Krasnoselskii [18] extended this result to an arbitrary solvable 
equation of the type (3.1) with a linear bounded operator A in a Hilbert 
space H. A theorem on the convergence of successive approximations 

il;n+i = {I -iyAi)iJn + i^fi (3.11) 

to the solution is proved. Here, Ai = A* A; fi (x) = {A*f ) (x); I is the 
identity operator; A* is the conjugate operator to A; < z/ < 2/ \\Ai\\; 

ipo {x) e H. 

Note that in the case of the integral operator (3.1) 
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where 

b 

a 

A number of procedures are known that improve convergence of 
iterations according to Fridman (see [10]). For example, under the 
conditions that are specified with respect to the procedure (3.10), 

V'n+1 (x) = — - ^ ipn (x) , (3.12) 

^ + n=0 

where (fo (x) e L2 (a, b); 

b 

(fin (x) = (fin-l (x) + f{x)- j k (X, ^) (pn-1 (0 71 = 1, 2, ... . 

a 

G. M. Vainikko and A. Y. Veretennikov [19] studied an iteration 
algorithm of an implicit type: 

b 

aipn+i {x)+ j k {x, (0 = oitpn (x) + f (x) , n = 0, 1, . . . , 

(3.13) 

where ■0o (x) & L2 (a, b); the parameter a> 0. 

Note that in contrast to the regularization of the type (3.6), based on 
the smallness of a, the considered approach is characterized by multiple 
iteration with, on the contrary, sufficiently large value of this param- 
eter. Moreover, one of the merits of the procedures (3.10)-(3.13) is a 
possibility of a constructive application of an a posteriori estimate of 
the error to accomplish the iteration. 

In the simplest case, one finds the number n for which for the first 
time 

llV'n+l - V'n|L,(a,6) < Ci5 + C27, 

where 6 and 7 are errors in the determination of f {x) and k{x,^), 
respectively; ci, C2 are constants meeting a number of requirements to 
ensure the stability of the procedures of computation. 

The influence of errors, small in a probabilistic sense, on the con- 
vergence of successive approximations is also investigated [19]. 
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The authors of [20] gave arguments for usefulness of the combination 
of the regularization of the equation of the type (3.1), whose parameter 
is the number of iterations, with algorithms of the saddle-point type. 
This approach has its origin in the pubhcation by V. M. Pridman [21] 
and is realized, in particular, according to the scheme 

V'n+l = - PnA* (A^n - /) , (3.14) 

where 

_ \\A*{A^l;^-f)f 

^" \\AA*{A^P^-f)t 
which is adequate to the choice of the step of the descent from the 
condition of the minimum of the error of closure 

An+i = WAtjjn+l - /II ; II. II = \\.\\L,{a,b) " 

3.4 Inverse problems for differential equa- 
tions of mathematical physics 

Modern viewpoints on this subject are characterized in the monograph 
by O. M. Alifanov, E. A. Artyukhin and S. V. Rumyantsev [20]. In the 
procedure of mathematical formulation of the problems, structural and 
parametric identification is emphasized, which implies, respectively, a 
qualitative description of the considered processes by means of differ- 
ential operators and allotting quantitative information to the model. 
Interpretation in terms of causality is also given. The cause includes 
boundary and initial conditions with their parameters, coefficients of 
the differential equations and also the domain of the problem. The ef- 
fect reflects the status of the investigated object and represents, mostly, 
fields of physical quantities of different types. 

The restoration of the cause from the information about physical 
fields is considered as an inverse problem. A key consideration is as 
follows (p. 19): "A violation of a natural causal relation that takes 
place in the formulation of the inverse problem can lead to its mathe- 
matical incorrectness, such as, in most cases, instability of the solution. 
Therefore, inverse problems constitute a typical example of incorrectly 
formulated problems" . 
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According to the sought function, the following types of inverse 
problems of the identification of physical processes for partial differen- 
tial equations are singled out: 

1) Retrospective problems: the determination of the prehistory of a 
certain state of the problem. 

2) Boundary problems: the restoration of boundary conditions or 
of the parameters contained therein. 

3) Coefficient problems: the restoration of the coefficients of the 
equations. 

4) Geometrical problems: the determination of geometrical charac- 
teristics of the contour of the domain or of the coordinates of points 
inside. 

A principal difference between inverse problems of identification and 
those of regulation is pointed out, concerning the width of classes of 
possible solutions. Whereas in the former case their increase leads to 
complications in the numerical realization, in the latter case, on the 
contrary, this is a favorable factor. By the way, the algorithmic means 
[20] are almost completely based on the methods of the solution of 
integral equations of the first kind, to which the considered problems 
of heat exchange are reduced. 

In the formulation of inverse problems of mathematical physics, the 
proof of corresponding theorems of existence and uniqueness is of cru- 
cial importance. In this regard, a general approach, outlined schemati- 
cally by A. L. Buchgeim [22], is of interest. Thus, the following formu- 
lation is considered: 

Pu^f; Qf^g, (3.15) 

where P is an operator of the direct problem; Q is an "information" 
operator describing the law of the change of the right-hand side; g 
is given, whereas u and / are the sought elements of corresponding 
functional spaces. 

The apphcation of the operator Q to the first equation (3.15) yields 
QPu — g, which is equivalent to 

PQu = [P, Q]u + g, 

where [P, Q] = PQ — QP is the commutator of the operators P and Q. 
The meaning of the commutation lies in the fact that, as a rule, there 
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is no information, except for (3.15), about the function /. Therefore, 
it is easier to study the operator on the solution of the direct problem 
u that satisfies some manifold of boundary conditions. It is important 
that in typical apphcations the operator Q does not " spoil" the part of 
boundary conditions that reflects the domain of the operator P. As a 
result, one gets a speciflc factorization of the inverse problem (3.15) as 
a product of two direct problems, induced by the operators P and Q 
under the condition that the commutator is, in a sense, "subordinate" 
to them. 

In the trivial case [P, Q] — 0, the initial problem decomposes into 
two simple ones: Pv = g; Qu = v. For the description of properties of 
the employed operators, a priori estimates are used. 

Of interest is also a quotation from the introduction to the mono- 
graph by R. Lattes and J.-L. Lions [23]: "In this book, we suggest a 
method of quasiinversion, intended for the numerical solution of some 
classes of incorrect, according to Hadamard, boundary value problems. 
Practical and theoretical importance of such problems is being more 
and more realized by investigators" . Moreover: The main idea of the 
method of quasiinversion (universal in numerical analysis!) consists in 
an appropriate change of operators entering the problem. This change 
is done by the introduction of additional difl^erential terms that are 

i) sufficiently "small" (they can be set equal to zero); 

ii) "degenerate on the boundary" (to prevent, for example, the ap- 
pearance of complicated boundary conditions and of such conditions 
that may contain unknown, sought variables)". 

In particular, the incorrect formulation of the problem of thermal 
conductivity 

dtu - dlu = 0, (3.16) 

u (0, t) = u {l,t) = 0; u {x, T) = K (x) , 

where k {x) is an unknown function, can be replaced by the following, 
with a small parameter e: 

dtU — d^u — ed^u = 0; (3-17) 

u — df^u — 0, x — 0; X = 1; u{x,T) = k{x). 

The authors point out (p. 36): "In a numerical realization, it is 
natural to choose e as the smallest possible one. However, in problems of 
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the considered type, one should expect numerical instability for e 0. 
Therefore one can expect at most that for any problem there exists a 
certain optimal value of e equal to eo" . The absence of convergence "in 
a usual sense" of the solution of the problem (3.17) to the exact one 
for e ^ was pointed out by A. N. Tikhonov and V. Y. Arsenin ([1], 
p. 52). 

3.5 Alternative viewpoints and develop- 
ments 

In Y. I. Liubich's opinion, any more or less general theory of integral 
equations of the first kind is absent, and only in some cases it is possible 
to use special methods. An example is given by known Abel's equation 
([24], p. 83). 

K. I. Babenko's remark ([25], p. 310) is rather typical: "Although 
from the point of view of the loss of information algorithms are not 
estimated, it seems to us that this is an important characteristic and it 
should be taken into account" . In what follows, the lack of optimality 
of the traditional approach to a numerical realization of incorrectly 
formulated problems is concretely demonstrated. 

A profound analysis of methodological aspects of this sphere is given 
by R. P. Fedorenko ([26], sections 40, 41). In particular, he failed to 
establish the value of the regularization parameter a by minimizing the 
functional (3.2), because for small values the sought function began to 
oscillate, whereas with it increase the value of considerably exceeded 
the admissible one. It is found that the reason is the inadequacy of 
the theory [1] to problems of control, characterized by discontinuity of 
solutions. 

In the context of investigations of the problem (3.16), R. P. Fe- 
dorenko brought up the following consideration: " All the methods of 
the solution of incorrect problems more or less consist in preventing 
the appearance in the sought solution of higher harmonics with large 
or even simply finite coefficients. But what is "high frequency"? Be- 
ginning with what number n should we consider the function sin (nirx) 
redundant, only spoiling the solution? This, of cause, depends on T" . 
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It is implied that a hypothetically known solution of the corresponding 
direct problem can be expanded into a Fourier series 



oo 



u{x) — u (x, 0) = 2^ sin (mrx) . 

n=l 

It is shown that the use in [23] of the value T = 0.1 and the errors 
in L2 (0,1) of the satisfaction of the last condition (3.16), with S of 
order 10^^, imposes the restriction n = 2. In this context, the method 
of P. Lattcs and G.-L. Lions came under criticism. This authors, while 
solving the problem (3.17) on a grid with a step of Ax = 0.02, obtained 
an absolutely unacceptable component Uq, namely, lO^sin (67ra;). This 
occurred for S at the level of 0.05, under the conditions when \k,{x)\ < 
1... . 

Note also the remark [26] that, aside the fact of the boundcdncss 
of the regularizing operator G (see section 3.1), its norm is an 
exceptionally important characteristic whose value directly influences 
a relation between the accuracy of the given function k and the solution 

uo = Gk. ^ 

Indeed, let us consider Eq. (3.6), written in the canonical form 

1 ? 1 
^l^ix)^— k{x,^)xl^{^)d^+-f{x), xe[a,h]. (3.18) 
a J a 

a 

Let a = 0, 6 = 1, and the kernel k (,x, ^) is determined by the expres- 
sion (2.8). In this case, its characteristic numbers and eigenfunctions 
are A„ = — (nvr)^, ipn {x) = sin {uttx), n = 1, 2, . . ., respectively. Thus, 
on the basis of general theory (see, e.g., [27]), for ^ A^, the solution 
of Eq. (3.18) is 

1 1 ~ c 

i^{x) 



^ -f (x) - - Yl r—, — 7 — ^2 sin (nnx) , 



where 

1 



Cfi — J f (x) sin (nnx) dx. 




^By the way, in most specialized publication this issue is not accentuated. 
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It is not difficult to notice that for small values of a the error in the 
determination of the function / (x) can considerably distort ip {x). 

In a constructive aspect, R. P. Fedorenko recommends to use tra- 
ditional formulations of inverse problems of differential or variational 
character with an application of additional conditions that rationally 
restrict classes of possible solutions. As the main factor to achieve the 
desired efficiency, a comprehensive analysis of qualitative peculiarities 
of solutions to the considered problems, involving elements of numerical 
simulations, is suggested. 

What arc the values of the regularization parameter a, typical of 
computational practice? The authors of [28] point out that for problems 
of restoration of time- dependent density of thermal flux on the surface 
from the results of temperature measurements at internal points of the 
samples the corresponding range is rather representative: 10~^ — 10~^. 
The editors of the above-mentioned book have a different point of view: 
" One can give a lot of examples of solutions to inverse problems thermal 
conductivity, when the range of acceptable values of a turn out to be 
rather narrow" (p. 141). 

The main technique of a numerical realization [28] is interpreted 
by the authors as a complement to the method of least squares by 
a procedure that smooths oscillations of the solutions in high order 
approximations. In this regard, they point out a relationship between 
Tikhonov's regularization and algorithms of singular expansions and 
ridge regression (or damping) that are widely used for the suppression 
of the instability of the method of least squares [29] . 

In a number of publications, one can see an orientation towards reg- 
ularization of Eq. (3.1) without the distortion of the original operator 
along the fines of (3.4) or (3.6). Thus, A. P. Petrov [30] suggested a 
formulation of the problem with f (x) & R (A) by means of the represen- 
tation / = A'lp + uj, where is a random process reflecting errors of the 
data and of the calculations. At the same time, the author failed to use 
his formally achieved correctness to construct an efficient algorithm of 
a numerical realization. The reason, in our opinion, is non-constructive 
structure of uj from the point of view of adaptive compensation of the 
error of closure of the satisfaction of (3.1). 

A. V. Khovanskii [31] put forward arguments for the regulariza- 
tion of the algorithm of the solution of Eq. (3.1), not the operator A 
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(which is the basis of the theory of [1]). The following quotation is 
of interest: "Tikhonov's regularization contains in an inseparable form 
two completely different notions, accuracy and stability, and there is a 
transformation of one into another. Nevertheless, there exists for a long 
time an idea of the predetermination of the operator [32], although only 
in the context of conjugate gradients and in a multiplicative form". 

In fact, this implies the use of Eq. (3.6) with a = ao, where ao is 
the minimal value allowing us to obtain information about the unknown 
function ip (x) of indirect character. With its use, one can determine the 
regularization parameter a matching the error of the free term / (x). 

However, the method of conjugate gradients is, in fact, Fridman's 
iterations of the type (3.14). Note that nonlinearities contained therein 
facilitate the smoothing of a well-known slow-down of the convergence 
of the procedure (3.10) with approaching the solution to Eq. (3.1). 
This effect was demonstrated by A. D. Myshkis [33] with the help of 
the representation of the components of (3.10) by scries in terms of the 
eigenfunctions of the kernel k{x,^). This leads to the relations 

Cn+l,m — ~ ^/ ^m) Cn,m + ^frm 171 — 1,2, ... , 

where c„ „ and are coefficients of the above-mentioned expansion of 
ipn {x) and / (x), respectively. 

When the number of the terms in the representation of the solution 
increases, which seemingly had to improve the accuracy, the coefficient 
of convergence (1 — A/A„i) approaches unity and, as a result of the 
accumulation of errors, the iterations become counterproductive. 

Note an effective method of the suppression of instability of the al- 
gorithm of the solution of the Predholm integral equation of the second 
kind 

b 

i;{x)^\jk{x,i)i;{i)di + f{x), xe[a,h], (3.19) 

a 

"positioned on the spectrum", i.e., in the case A = A„, developed by P. 
I. Pcrlin ([34], pp. 105-107). 

This problem is incorrect both with respect to the uniqueness of 
the solution and as a result of the degeneracy of the system of linear 
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algebraic equations obtained by discretization. Nevertheless, a pertur- 
bation of the right-hand side of / (x) by a zero (within the limits of the 
accuracy of calculations) component 



where i/j^ (x) is a normalized eigenfunction of the conjugate kernel 
k {^,x), allows one to improve radically the situation. 

The essence lies in the fact that, theoretically, the solution to Eq. 
(3.19) is expanded in a power series of A. Provided that computational 
procedures, matching this situation, are identical, one can compensate 
for the errors. 



3.6 A comparison between the main con- 
cepts of A. N. Tikhonov and V. M. 
Pridman 



A. N. Tikhonov's original suggestion (1943) admitting of the consider- 
ation of incorrectly formulated problems by an a priori restriction on 
the class of possible solutions is a kind of refraction of general method- 
ology of investigations of the issues of existence and uniqueness into 
the sphere of numerical analysis.^ Note that A. N. Tikhonov's proof of 
the well-known theorem on the uniqueness of the solution of the inverse 
problem of thermal conductivity in an infinite n-dimensional domain 
under an additional condition of the type \d^u\ < M dates back to 
1935. A clear illustration of these considerations is provided by the 
algorithm of the search for a quasisolution (3.7)-(3.9) that artificially 
subjects the data of Eq. (3.1) to conditions of the type of those that 
figure in Picard's theorem. 

^There is a translation of the first edition of [1]: A. N. Tikhonov and V. Y. 
Arsenin, Solutions of Ill-Posed Problems (Winston, Washington, 1977). See also V. 
A. Morozov, Solutions of Incorrectly Posed Problems (Springer, New York, 1984). 
An English translation of the article of V. M. Fridman [17] is given in Appendix. 



b 




a 
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Behind A. N. Tikhonov's method of rcgularization (1963), there is 
a global idea of a limiting transition to the exact solution with respect 
to a small parameter of the problem, which is unambiguously pointed 
out in ([1], p. 56): "Note that regularizing operators, dependent on a 
parameter, have been employed in mathematics since Newton's times. 
Thus, the classic problem of an approximate calculation of the deriva- 
tive u' (x) by means of approximate (in the metric C) values u (x) can 
be solved with the help of the operator 

^ ^x + a^-^a)„_ 

a 

Then, instead of the exact value of the function u (x) , an approx- 
imate one Us (x) = u{x) + Au{x) with |AM(a;)| < 5 is substituted. 
On the basis of these calculations, one makes the statement: "If a = 
S/rj {S), where r]{5)^0 for 5^0, then 25 /a^ 2r] {5) ^ for 5 ^ 0. 
Thus, for a = cxi {8) = 5/rj (S), G {us, ai {5)) u' (x)" . 

It should be noted that, using the methodology of a small parame- 
ter, A. N. Tikhonov obtained fundamental results in the field of inves- 
tigations of differential equations with a singular perturbation of the 
type 

eii^ f {u,v,t); v^g{u,v,t), 

where e is a small parameter; f {u,v,t) is a nonlinear function (1948)- 
(1952)3. 

The solution of the system of equations does not depend continu- 
ously on the parameter e. Proceeding to the limit e ^ creates a new 
object of investigations with completely different properties. In the first 
place, it imphes the issue of the so-called violation of the stability of 
the root of the equation f {u,v,t) = 0. Nevertheless, A. N. Tikhonov 
managed to develop a rather constructive theory that served as a basis 
for a number of productive approaches of both fundamental and ap- 
plied character. The importance of A. N. Tikhonov's achievements in 
the sphere of system analysis is analyzed in detail by N. N. Moiseev 
([36], section 5). 

However, properties of the integral equation (3.6) for a = also 
change radically. In this regard, generally speaking, a certain analogy 

^See the review by A. B. Vasileva [35]. 
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emerges. One can suggest that A. N. Tikhonov undertook an attempt 
to use the techniques of his theory of singular perturbations for the 
solution of incorrectly formulated problems. 

This suggestion is supported by the following quotation from the 
monograph by S. A. Lomov ([37], p. 12): "Now it is becoming clear how 
to isolate in singularly perturbed differential equations small terms that 
can be neglected. It turned out that one needed additional information 
about the solution to do this." 

Note J. Hadamard's remark that an extension of methods of the 
theory of ordinary differential equations to problems of mathematical 
physics should be done with great care ([38], p. 38). At the same 
time, at the turn of the 1950s, the theory of singular perturbations 
became an efficient tool in investigations of complicated problems of 
partial differential equations (publications by M. I. Vishik and L. A. 
Liusternik, O.A. Olejnik, K. O. Fiedrichs, and others). By the way, 
explaining the conceptual basis of their method of quasiinversion, R. 
Lattes and G.-L. Lions ([23], p. 11)^ refer to these authors and A. N. 
Tikhonov. 

Simultaneously, they pointed out that A. N. Tikhonov's priority 
publication on the method of regularization [8] (see also [39]) was pre- 
ceded by D. L. Phillips' article [40], whose results with respect to inte- 
gral equations were analogous. In the monograph by F. Natterer [41] 
this regularization figures as Tikhonov-Phillips' method. V. A. Moro- 
zov estimated the achievements of the latter author in a much more 
restrained manner ([6], p. 10): "Some recommendations on the use 
of this method are contained in the publications by L. V. Kantorovich 
[42] and D. L. Phillips [40]. There is no theoretical justification of this 
approach in the above-mentioned publications". 

The chronological reference to the most important results in the 
field of the construction of stable algorithms for the solution of integral 
equations of the first kind ([10], p. 234) gives the following informa- 
tion: "1962, Phillips's publication [40], where he suggested a variational 
method of conditional minimization of the functional (with the use of 
restrictions on the smoothness of the solution) and put forward the idea 



^Ideological closeness of quasiinversion and Tikhonov regularisation was pointed 
out by M. M. Lavrentiev [23, p. 5]. 
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... of a choice of the regularization parameter a". 

Turning to V. M. Fridman's achievements, note that it is rather 
difficult to evaluate the premises that form the basis of the iteration 
procedure (3.10). At the first sight, such a computational method has 
a lot of analogs. However, its adequacy, in a sense, to the object of 
investigation, the incorrect problem for the Frcdholm integral equation 
of the first kind, turned out to be rather unexpected. 

Later on, with the aim to improve convergence, V. M. Pridman 
also employed the nonlinear algorithm (3.14). In our opinion, different 
ways of the determination of the number of the final iteration and of the 
increase of the rate of global convergence (sec [14, 19, 20]), despite their 
actuality for practical application, should be interpreted as a technical 
complement to V. M. Pridman's methodology. 

Nowadays, the algorithm of conjugated gradients is considered to 
be nearly the most efficient one for the solution of large ill-defined 
sparse systems of linear algebraic equations, obtained by the reduction 
of, apparently, of most problems of numerical simulations [32, 43, 44]. 
As is pointed out by J. Ortega [32], this method was proposed by M. 
P. Hestens and E. L. Stiefel (1952). However, for certain reasons, it 
was not employed for a long time. It attracted considerable interest 
at the turn of the 1970s, when one realized the actual sphere of its 
applications, the potential of the above-mentioned predetermination 
and adaptivity with respect to paralleling of computational operations 
in combination with the architecture of modern computers. 

Thus, the priority of the method of conjugated gradients ensured 
its refraction to a class of problems of linear algebra, characterized 
by the instability of the numerical realization, that is, in fact, incor- 
rectly formulated. In this regard, we emphasize that V. M. Fridman's 
"methods of the saddle-point type" [21] can be interpreted as somewhat 
simplified representatives of the family of the methods of conjugated 
gradients ([20], section 2.1; [43], section 7.1). It seems that V. M. Prid- 
man, who was the first to use systematically iterations for the solution 
of incorrectly formulated problems, essentially foresaw the development 
of computational mathematics that followed. 

In light of above, the position of M. A. Krasnoselskii and the co- 
authors is worth noting [18]. The role of V. M. Pridman in the devel- 
opment of the iteration procedure (3.11), which is an analog of (3.10), is 
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described as follows: "A transition to the equation [ip = {I — vAi) ip + z//i] 
was pointed out for some cases by I. P. Natanson [45]. For Fredholm 
integral equations of the second kind, it was already employed by G. 
Wiarda [46] . For integral equations of the first kind, it was, essentially, 
employed in the publication by V. M. Fridman [17]" (p. 73). There 
is no comment on a qualitative difference between the objects of the 
investigation. 

The nontriviahty of V. M. Fridman's approach is noted in the re- 
mark of I. P. Natanson [45] : " Our method does not apply to the solution 
of the integral equation of the first kind. This could be expected, be- 
cause the use of the method implies complete arbitrariness of the free 
term of the equation Axp = /, whereas Eq. (3.1) is solvable not for 
all / (a;)". In what follows, the author gives an extended proof of the 
degeneracy of the corresponding discrete problem. 

The gradient algorithm of V. M. Fridman [21] is mentioned by the 
authors of [18] exclusively in the context of the equation Aip = f, where 
both the operators A and A~^ are bounded (p. 115). We quote the 
abstract to V. M. Fridman's paper [21]: "We present a new proof of the 
convergence of methods of the saddle-point type for a linear operator 
equation. We do not assume, unlike L. V. Kantorovich [47], M. A. 
Krasnoselskii and S. G. Krejn [48], that zero is an isolated point of the 
spectrum of the operator".^ 

3.7 Ill-defined finite-dimensional problems 
and issues of discretization 

In this subsection, Aip = f denotes a system of linear algebraic equa- 
tions. The conditionality number of the matrix A (see, e.g., [49]) 

cond(^)=ma.J!M/„i„ll-'''/'ll 



where ^0 is a manifold of vectors of the Euclidean space, represents a 
raising coefficient between a relative error of the data and the solution. 
At the same time, cond(y4) characterizes the measure of closeness of 

^This is equivalent to the boundedness of the operator A~^. 
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A to a degenerate matrix, for which the sohition of the corresponding 
system of algebraic equations does not exist or is nonunique. 

An algorithm of the solution of a degenerate system of linear alge- 
braic equations, based on the method of least squares, is presented in 
the book by A. N. Malyshev [50]. First, the matrix A is transformed 
to a two-diagonal one by means of a special transformation, and one 
finds its eigenvalues that arc subdivided into two groups, cxi, (T2,..., cr„ 
and (7„+i,..., such that an/ (cr„ — cr„+i) is not very large. Then, with 
the help of a rather laborious procedure of the exhaustion of the sec- 
ond group of the eigenvalues, one constructs a matrix A^ that is stably 
invertiblc beginning with a certain value n. The accuracy of the thus 
obtained generalized solution tjj is determined by the error of closure 
~ / / II^IL using heuristic considerations. 

It seems that in a methodological sense this scheme reminds of B. K. 
Ivanov's algorithm [15] that reflects computational relations (3.7)-(3.9). 

L. Hageman and D. Young [43] studied the approach of predetermi- 
nation, employed for the solution of systems of linear algebraic equa- 
tions, close to degenerate ones, to accelerate by the method of conju- 
gated gradients iterations of the type 

where P — I — Q~^A] g — Q^^f. It is assumed that this procedure 

can be symmetrized in the sense that there exists a non-degenerate ma- 
trix W such that the matrix W {I — P) is symmetric and positive 
definite. 

By use of W, the initial problem can be reduced to the solution of 
much better defined systems of algebraic equations B(f — q, where 

B = W {I - P)W-^; ip = Wilj; q = Wg. 

Formally, a choice of the predeterminer does not pose problems. 
However, in practice, one has to resolve a contradiction between the 
conditions imposed on the matrix W: "closeness'' to A^'^ to reduce the 
number of iterations; a "rapid" calculation of a product of the type 
Wip [51]. In the above-mentioned publication, I. E. Kaporin analyzes 
different approaches to the construction of predeterminers for systems 
of linear algebraic equations of a general type. An analogous issue, in 



3. 7. ILL-DEFINED FINITE-DIMENSIONAL PROBLEMS AND ISSUES OF DISCRETIZATIOmi 



the interpretation of J. Ortega [32], is oriented mainly towards sparse 
matrices. 

The complexity of problems of linear algebra that arise in the re- 
alization of modern methods of investigations in the field of the me- 
chanics of a continuous medium are characterized as follows [51]: "The 
matrices of corresponding systems are rather large (up to a hundred 
thousand nonzero elements), rather densely filled (up to hundreds or 
even thousands of nonzero elements in each line), have no diagonal pre- 
dominance, are not M-matrices and are rather ill-defined. In general, 
one can expect only symmetry and positive definiteness of the matrix 
of the system".^ 

Note that, for example, in seismic tomography [44], one has to be 
satisfied with a numerical realization of discrete analogs of integral 
equations of the first kind, because their kernels cannot be represented 
analytically and parameters of the considered models are determined 
with the help of natural experiments. 

In light of the above, the considerations of R. W. Hamming ([52], 
p. 360) may seem to be archaic: "A system of linear equations is said 
to be ill-defined, if, roughly speaking, the equations are almost linearly 
dependent. Many efforts were made to investigate the problem of the 
solution of ill-defined systems. However, one may pose the question: Is 
it necessary to solve such systems in practical situations? In what phys- 
ical situation may the solutions prove to be useful, if they depend in 
such a substantial manner on the coefficients of the systems? Usually, 
the following is true: Instead of the solution, one is looking for a system 
of almost linearly independent equations. In light of this information, 
the problem can be better understood and is usually reformulated again 
in a more satisfactory way. It is rather probable that ill-defined sys- 
tems of equations, provided that round-off and measurement errors are 
eliminated, are actually linearly dependent and thus do not reflect the 
physical situation". 

Note that, in contrast to the above-mentioned constructors of meth- 
ods of computational mathematics, the renowned practitioner adheres 
to the position of correctness according to Hadamard. Let us quote P. 



^The non-diagonal elements of an M-matrix are non-positive, and all the ele- 
ments of its inverse are non-negative. 
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S. Gutcr's preface to [52]: "The name of R. W. Hamming, a renowned 
American scientist, former President of the Computer Association, Head 
of the Mathematical Service of BeU Telephone Laboratories, and his 
works in the field of computational mathematics and the theory of 
information arc rather well-known and do not need special recommen- 
dations. ... The book 'Numerical Methods for Scientists and Engineers' 
is without any doubt an outstanding phenomenon in mathematical lit- 
erature" . 

Of special interest is R. W. Hamming's opinion about the priority 
of computational procedures ([52], p. 90): "It is often believed that the 
main problems of numerical analysis are concentrated on interpolation, 
but this is not the case. They are mostly related to such operations as 
integration, differentiation, finding zeros, maximization, etc., in those 
cases when all we have or can compute are some nodes of functions 
that are usually known not exactly, but approximately, because they 
are spoiled by the round-off error". 

Thus, the problem should be posed correctly under the conditions 
of an inevitable error in the data, which is equivalent to the preference 
of algorithmic efficiency to the quality of initial information. Interpo- 
lation, mentioned in the above quotation, implies approximate repre- 
sentation of the latter for the performance of computer operations by 
means of a finite-dimensional approximation. 

However, in computational mathematics, alternative concepts arc 
rather wide-spread, which is refiected in K. I. Babenko's remark [25]: 
"In some spheres of numerical analysis, the theory of approximation 
serves as the foundation for the building of the numerical algorithm" 
(p. 138). "Information, inputted into the algorithm, is characterized, 
in the first place, by its volume... All other characteristics, such as, 
e.g., accuracy, are its derivatives and do not present a true picture of 
the input" (p. 281). 

Here, information is understood in the sense of Kolmogorov's the- 
ory of e-cntropy that identifies it with the length of a given table or 
an alphabet, whose words are manipulated by the algorithm. Corre- 
spondingly, the issue of numerical analysis is interpreted in terms of, 
figuratively, the deficiency in the search for necessary words and of the 
deletion of tables in the course of operations. 

Nevertheless, R. W. Hamming's point of view on the relation be- 
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tween the method of investigations and the employed information is 
actively developed by a group of specialists with J. Traub and G. 
Wasilkovski at the head. The authors of [53] point out (pp. 9, 6): 
"In this book, we construct a general mathematical theory of optimal 
reduction of uncertainty. Wc interested in the two main questions: 1) 
Is it possible to reduce imcertainty to a given level? 2) What will it 
cost? The aim of the theory of informational complexity is to provide 
a unified approach to investigations of optimal algorithms and their 
complexity for the problems that involve incomplete, imprecise or paid 
information and to employ the general theory to concrete problems 
from different fields". 

Here, complexity implies the number of arithmetic operations, the 
time of their realization, computer memory resources, etc. It is interest- 
ing to note that the interpretation of the notion of probability [52, 53] 
correlates with the expressive statement of R. Bellman and S. Dreyfus 
([54], p. 342): "Fortunately, in some cases, there is a very simple way to 
overcome this difficulty. Instead of trying to study information as the 
"smile of Cheshire Cat", we consider the actual physical process, where 
information is used to work out solutions.^ The value of information 
can then be measured by the efficiency of the solutions. 

Thus, the usefulness of information depends on its application, which 
is the most reasonable concept!" 

It should be noted that the procedure of finite-dimensional approx- 
imation of problems of mathematical physics is, of course, also very 
important, which is accentuated by K. I. Babenko. Indeed, the ob- 
tained discrete model can turn out to be incorrect, and the employed 
algorithms of the numerical realization may prove to be divergent even 
in the solution of rather ordinary problems. An example of instabil- 
ity of a finite-difference scheme is given by S. K. Godunov and V. S. 
Ryabenkn ([55], section 4.9). 

K. I. Babenko also emphasized the absence of any general methods 
of the construction of finite-dimensional analogs ([25], p. 622): "...the 
provision of an approximation alone is insufficient"... one has to ensure 
that the discrete problem "retains the type of the original continuous 

'^The smile of Cheshire Cat, according to L. Carrols "Alice in th Miracle Land", 
existed separately from this cat (editor's note to [54]). 
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problem". In his opinion, to achieve the above goal, "a detailed inves- 
tigation in each concrete case is required, which is the most nontrivial 
part of work" . 

3.8 The crisis of the technology of numer- 
ical simulations 

Of considerable interest is, in fact, a program statement of O. M. Be- 
lotserkovskii and V. V. Stchennikov in the preface to [56]: 

"A rapid development of computers, especially during the last 10-15 
years, with a special acuteness posed the problem of the construction of 
a principally new technology of the solution of problems by computers. 
... Historically, the problems of numerical simulations (in this notion, 
we include the actual mathematical simulations related to a numerical 
experiment), being rather advanced already in the "precomputer" pe- 
riod and rapidly developing during the next periods, turned out to be 
the most conservative component of the modern technology of the so- 
lution of problems on the computer. Using, probably, redundant from 
the point of view of a mathematician expressiveness of the description, 
one can characterize the existent situation by two stable tendencies: 

- an increase of the complexity of mathematical models; 

- construction of rather sophisticated mathematical methods. 

Both the tendencies inevitably lead to a technological deadlock, be- 
cause they create complications in the solution of the problem of the 
construction of software-hardware means of the support of the whole 
technological chain. ... Without any pretension to profoundness and 
importance of the analogy, we dare say that the present situation in 
numerical simulations is similar to that in mechanics before the ap- 
pearance of main ideas and concepts of quantum mechanics" . 

In the introductory article [56] the same authors emphasize the phe- 
nomenon of the accumulation of the round-off error in the numerical 
realization of algorithms that include up to 10^^ operations and the 
absence of real means to estimate the error of solutions to, in partic- 
ular, evolution problems. In their opinion, "...the following conclusion 
is quite justified: a priori, any evolution problem for large times is nu- 
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merically (or computationally) incorrect in the sense of the absence of 
a practically important solution... 

In the case, when a priori or a posteriori information about the error 
of an approximate solution is absent, it is impossible to claim that the 
solution exists. This conclusion fairly agrees with A. N. Tikhonov's 
theorem that states that the problem with the data on the operator 
and the right-hand side has no solution in the manifold of approximate 
numbers" . 

O. M. Belotserkovskii and V. V. Stchennikov regard as construc- 
tive the idea that discrete models of the considered problems should 
be assembled with the aim of increasing the accuracy of information 
by means of special superposition. They also suggest to search for the 
solution in the class of function with a bounded variation, with would 
endow the difference operator of the problem with smoothing proper- 
ties. 

As is well-known, N. N. Yanenko paid considerable attention to the 
methodology of mathematical simulations (see [2]). His concept of over- 
coming the above-mentioned crisis is explained by O. M. Belotserkovskii 
( [57], p. 106): 

"An investigation of finite-difference schemes, approximating differ- 
ent classes of equations of mathematical physics, led N. N. Yanenko 
to an extension of the notion of the scheme. For the first time, he 
begins to consider the finite-difference scheme as an independent ob- 
ject of the investigation, as a mathematical model, adequate to this or 
that physical model. This fundamental concept is based on profound 
understanding of the fundamentals of differential and integral calculus. 

Indeed, physical and mathematical models, described by differen- 
tial, integral or integrodifferential equations, are obtained from discrete 
models by means of averaging and passing to the limit with respect to 
certain parameters. This is the case, for example, in the model of a con- 
tinuous medium, where for a sufficiently large number of elements in 
the unit volume one comes to the notion of the continuous medium by 
averaging and passing to the limit with respect to the volume. In this 
regard, one can interpret a finite-difference scheme as an independent 
mathematical model with certain properties" . 

Note the fundamental, as it seems, considerations of N. N. Yanenko 
[2]: "The objects of modern mathematics, whose theoretical "nucleus" 
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comprises topology, geometry, algebra and functional analysis, are ideal 
logical constructions forming a certain operational system. We will call 
them ideal objects, which underlines, on the one hand, their practi- 
cal inaccessibility and, on the other hand, their excellent operational 
properties that allow one to make operations without loss of informa- 
tion. Ideal objects of mathematics are essentially infinite and require 
an infinite number of operations" (p. 12). 

"The development of the experimental foundation and the tool of 
investigations, the computer, increased interest in such objects as com- 
puter numbers, programs, finite automata. In this regard, the definition 
of mathematics as studies of the infinite, accepted in the 20th century, 
should be replaced by another one, more correctly refiecting its essence, 
i.e., as studies of the relationship between the finite and the infinite" 
(p. 18). 

We think that on the basis of the above one can come to a very 
important conclusion: In their construction of the conceptual basis of 
mathematical simulation, the leading specialists were guided by the 
concept of inapplicability of Banach's theorem on the inverse operator. 
Note that N. Dunford and J. Schwartz considered this theorem as one 
of the three principles of linear functional analysis, characterized as 
being rather fruitful ([58], p. 61).^ 

A quotation from K. Maurin's manual ([59], p. 51) reads: "This 
theorem [on the closed graph], in the last years, has gained itself a 
reputation of being the most important theorem of functional analysis, 
if this one is considered from the point of view of applications" . 

An attempt to renew the above-mentioned fundamentals in the 
context of the accentuation of peculiarities of computational mathe- 
matics was made by A. V. Chechkin [60], who suggested a division of 
sections of mathematics into classical and non-classical ones, respec- 
tively: "arithmetics, mathematical analysis, algebra, geometry, proba- 
bility theory, etc.; mathematical logic, the theory of information and 
statistics, the theory of sets, the theory of algorithms and recursive 
functions, methods of computational mathematics, the theory of finite- 
difference schemes, the theory of cubic formulas, methods of the solu- 



^The other two are the principle of linear boundedness and Hahn-Banach's 
theorem. 
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tion of incorrect problems, etc." (p. 8). As a criterion, the authors 
choose the fact of availabihty of absolutely complete or partial infor- 
mation about the considered objects (points, functions etc.). 

Let us quote the abstract of section ([60], p. 78): " We define 
and study a new type of mappings that generalize classical notions. 
Classical mappings realize correspondence between the points of a set. 
This implies that the points are known with absolute precision. The 
new mappings, termed ultramappings, realize correspondence between 
pieces of information about points of sets. The main construction of 
the ultramappings, termed ultraoperators, allows one to obtain separate 
information about the image point from separate information about the 
inverse image point. 

Ultracontinuity of ultraoperators is defined, which is a broad gener- 
alization of the notion of the stability of methods. It is found that, for 
an arbitrary base operator, one can construct an ultracontinuous oper- 
ator over it. A class of ultracontinuous operators, termed Tikhonov's 
operators, is singled out. For these operators, the base operators are 
not continuous". Furthermore, "they are related to A. N. Tikhonov's 
ideas and methods of the solution of incorrect mathematical problems" . 

Returning to the question of adequate discretization, wc quote the 
abstract of the monograph by A. A. Dezin [61]: "It is devoted to the 
description of the basic structures of multidimensional analysis and 
to the consideration of internally defined discrete problems of analysis 
and mathematical physics. It implies not merely an approximation of a 
given continuous object, but the construction its analog, starting from 
the notion allowing for discrete interpretation". 

Arguments for contradiction to physical sense of differential models 
of certain classes of problems of the mechanics of a continuous medium 
are given by M. A. Zak [62]. In this regard, he developed a general 
approach, wholly based on the concepts of theoretical mechanics with 
a special interpretation of Gauss' principle of least action. 

The position of C. Truesdell is alternative. He thinks that contin- 
uum mechanics of a deformed body "is, in essence, not only subtler, 
more beautiful, majestic than a rather sparse particular case, called 
" analytical mechanics" , but it is much more suitable for the simulation 
of real bodies" ([63], page 10). 
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Chapter 4 



Short comments on the 
material of the above 
sections and some general 
considerations 



4.1 The correctness of the formulation of 
problems of mathematical physics 

The conditions of correctness, formulated by J. Hadamard at the turn of 
the 20th century (see [1]) and insistently advocated by him thereafter, 
amaze us by their ever- increasing importance for practical applications. 
It forms the conceptual basis for methods of the simulation of physi- 
cally meaningful problems, which, in fact, is disputed by nobody. At 
the same time, nowadays, the prevailing opinion is that Hadamard's 
concepts are principally invalid. 

This concerns the basic statement that the properties of existence 
and uniqueness, considered by Hadamard as inherent to mathemati- 
cal models of real processes, lead to the correctness of the formulation 
of adequate boundary-value (initial-boundary-value) problems, which 
implies the stability of the employed algorithms of a numerical realiza- 
tion. Although, of course, not all generally quite valid algorithms can 
ensure continuous dependence of the solution on the data of a concrete 
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problem: Implied is a possibility of their construction. 

The above-said concerns, in the first place, the key notion of this 
statement, namely, correct formulation of problems of mathematical 
physics. Indeed, only the existence of such formulation is claimed un- 
ambiguously, which implies, in particular, that the Fredholm integral 
equation of the first kind is simply unsuitable as a reflection of objective 
reality. 

Naturally, the course of investigations with the aim to confirm or 
disprove the hypothesis, or, in our opinion a prophecy, of J Hadamard, 
seemingly had to be conducted from the position of variability of formu- 
lations of the considered problems, which was not the case. The main 
reason is, in our opinion, an erroneous understanding of a special mis- 
sion of computational means of numerical simulations that lightheart- 
edly neglected even one of the main principles of functional analysis, 
i.e., Banach's theorem on inverse operator ([2, 3],section 9, and [4]). 

One can hardly explain the absence in special literature even of 
a thesis of the necessity to coordinate the formulation of problems of 
mathematical physics with algorithms of their numerical realization. 
The roots of this situation seem to be in systemic character of the 
giant computer-supply complex oriented at commercial efficiency at 
the expense of high costs of provided services. 

As a result, the alternative school of A. N. Tikhonov builds up the 
criticism of J. Hadamard ideas according to the following scheme: 

- the solution of the Fredholm integral equation of the first kind 

1 

{Ai;) (x) =Jk {x, iO d^ = fi^), ^ e [0, 1] (4.1) 



is, in general, an incorrect problem (which is undispTitablc); 

- integral equations of this type are adequate to a variety of real 
phenomena, which is actually supported by a rather transparent in- 
terpretation of corresponding direct problems (calculations of / from 
given k and 

However, what are the grounds for the formulation of the problem, 
inverse to the calculation of /, by means of trivial renaming the given 
and sought functions in (4.1)? The fact that the procedure of the 
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restoration of ip {x) for given f (x) and k{x,^) is incorrect does not 
imply any consequences. 

The Predholm integral equation of the first kind is, without any 
doubt, a rather graphic object for an illustration of the suggested 
arguments. As regards problems of mathematical physics for partial 
differential equations, their evaluation is not so unambiguous. Thus, 
Hadamard's example (2.3), (2.4), concerned with the restoration of the 
function u {x, t) from the initial data, is interpreted as a direct problem, 
whereas the inverse one, in contrast, has to be related to the determina- 
tion of the right-hand sides from the corresponding information about 
the solution. Hence, the formulation of the direct problem is incorrect. 

It is usually believed that the problem (2.3), (2.4) has physical appli- 
cations (see, e.g., [2]). In this case, following J. Hadamard"s method- 
ology, its formulation can be reduced to a correct one. As will be 
shown below, the difference between integral equations of the first kind 
and boundary-value (initial-boundary-value) problems of mathematical 
physics are, in a sense, conditional. Indeed, the problem (2.3), (2.4) can 
be reduced by elementary means to an equation of the above-mentioned 
type with respect to or d^u. 

The reproaches to J. Hadamard, whose typical elements are repro- 
duced in section 2.1, can be summarized as follows: The great scientist 
slowed down the progress of science by refusing to admit that incor- 
rectly formulated problems were adequate to a variety of real processes 
(see [3, 4, 5]). However, J. Hadamard had no doubts that such prob- 
lems existed. As a matter of fact, put forward the thesis of inadequacy 
of the employed means of numerical simulations. Indeed, is it justified 
to raise the issue of someone's monopoly of a mathematical model as 
regards its objective correspondence to the reaUties that we observe? 

Clearly, the answer can be only negative: Mathematical models 
are constructed with different degrees of detail, they have the form 
of integral equations, boundary- value problems, etc., inducing, by the 
way, rather non-identical difficulties of computational character. At the 
same time, the results of reliable versions of numerical simulations for 
the considered processes must, in general, be identical. 

One can raise an argument that a transition from an incorrect for- 
mulation of the problem to a correct one is a cardinal step and, ac- 
cordingly, should be accompanied by a substantial deformation of com- 
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putational relations, which in turn should influence the solution. This 
conclusion is seemingly supported by the well-known fact of functional 
analysis that the property of incorrect solvability is stable with respect 
to small perturbations [6]. However, this property is connected with 
a pair of spaces and, generally speaking, looses its validity when Eq. 
(4.1) is considered from the point of view of the mapping between the 
spaces L2 (0, 1) and I2, as it was mentioned above. 

It seems that we have pointed out an important issue. As a devel- 
opment of this issue, we state the absence of serious arguments against 
a possibility of the correct formulation of the problem, e.,g., of the 
restoration of the function ip {x) from the kernel k {x, ^) and the right- 
hand side / (x), obtained by integration according to Eq. (4.1). Purely 
heuristically, if small variations of the data can substantially influence 
the solution of incorrectly formulated problem, why not suggest that 
this effect can be overcome in a satisfactory way by a rather nonessen- 
tial correction of the traditional mathematical model? 

Among supporters of studies of problems of mathematical physics 
exclusively in the correct formulation are: A. Poincare, D. Hilbert, V. 
A. Steklov, I. G. Petrovsky, I. Prigogine [7, 8, 9, 10, 11]. On the other 
hand, the role of the three absolutely independent conditions of the 
correctness (existence, uniqueness and continuous dependence on the 
data of the problem), introduced by R. Courant and D. Hilbert [12], 
can hardly be called positive. We think that the constructive potential 
of the fact that the third condition is a corollary of the previous ones 
could facilitate the activation of studies of different aspects of the cor- 
rect formulation of problems of mathematical physics in the classes of 
functions inherent to them. 

When considering the Predholm integral equation of the flrst kind 
(4.1), one had to be more careful with respect to operations involving 
the function / (,x) from I2 as opposed, flguratively, to a surrogate of 
continuous inversion under the condition that it belongs to L2 (0, 1). 
Accordingly, an analysis of the premises of the correct solvability of this 
equation could lead to a conclusion that the problem of the restoration 
on its basis of the function ip (x) from the given / {x) and k {x, ^) cannot 
be interpreted as the inverse of the calculation of / (x) from Eq. (4.1). 
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4.2 A relationship to the theorem on the 
inverse operator 

The above-mentioned fact that the third condition of the correctness 
has the character of a corollary results from Banach's theorem on the 
inverse operator [13] whose optimistic meaning consists in the following: 
If the solution to Eq. (4.1), with D {A) = Bi and / e B2, where Si, B2 
are Banach spaces, exists and is unique, the inverse operator A"^ from 
B2 into Bi is bounded. As a consequence, the ordinary procedure of 
the solution of the Fredholm integral equation of the first kind (4.1) is 
by Picard's theorem stable under the condition that B2 coincides with 
k. 

However, in the general case both the verification and the fulfill- 
ment of the condition f E I2 are practically infeasible. Therefore, such 
spaces are called "inconvenient" (see [14, 15]). Thus, we are obviously 
in a principle dilemma as to the choice of the methodology of the in- 
vestigation: 

- an orientation at overcoming the difficulties resulting from the use 
of the space I2 related to the boundedness of the operator A^^; 

- the loss of this property in exchange for a possibility of studying 
mathematical models in "convenient" spaces. 

With the beginning of large-scale applications of computational 
methods to mathematical investigations, the second way became dom- 
inant. Correspondingly, the canonical formulation of discussed S. Ba- 
nach's theorem, given, e.g., by A. N. Kolmogorov and S. V. Fomin 
[16] (see section 2.3), has been accompanied by a specification that it 
is valid only in the case when the operator A realizes a mapping on 
the whole space B2 [17, 18]. It should be noted that for f E L2 (0, 1) 
this condition is rather important, because the space I2 inherent to Eq. 
(4.1) is a subspace or part of L2. 

The dynamics of the point of view of S. G. Mikhlin, reflected in 
his courses of mathematical physics and the theory of errors of 1968, 
1977 and 1988 [19, 20, 21]. At the beginning, the author considers Eq. 
(4.1) under the traditional assumption that the operator A is compact, 
that is, according to the type of the mapping inside the space L2 (0, 1). 
In this case, the inverse operator A~^ is unbounded. As a result, the 
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standard algorithms of the numerical realization are inapplicable, and 
one has to turn to the methodology of A. N. Tikhonov. 

Later, S. G. Mikhlin drew attention to the fact that if the Predholm 
integral equation of the first kind (4.1) is interpreted from the point of 
view of a mapping between the spaces L2 (0, 1) and I2, the corresponding 
operator A is no longer compact and, as a consequence, the inverse A~^ 
can be bounded, and the problem of the restoration of the function ip (x) 
becomes correct. In this way, the completeness of the conditions of 
correctness is restored, whereas the third condition was initially singled 
out by the author. 

Thus, the use of the pair of spaces L2 (0, 1) — I2 in a sense transfers 
the canonical incorrect problem to the mainstream of fundamentals of 
functional analysis. Note the fact that, being a mathematician, S. G. 
Mikhlin did not devalue the importance of the outlined step by banal 
reasoning in terms of "convenient/inconvenient" or "bad" and "good" 
spaces. 

Such a position apparently incurred criticism: In his concluding 
monograph, S. G. Mikhlin somewhat irritably readdresses actual for- 
mulation of problems of mathematical physics to specialists in applied 
sciences, including sociologists, who are interested in their solution. Si- 
multaneously, the author has found it reasonable not to consider infinite 
dimensional models with inherent aspects of incorrectness. 

The statement that a mathematical problem, being posed (in terms 
of premises), must be solved by rigorous methods belongs to V. A. 
Steklov. At the same time, why not consider the procedure of the for- 
mulation of problems of mathematical physics as an additional reserve 
of increasing the efficiency of employed techniques of numerical real- 
ization? Moreover, maybe rigidly predetermined formulations of prob- 
lems themselves prose artificial complications of computational charac- 
ter under the conditions when physical considerations admit a small, 
in a sense, variation? In our opinion, the formulation of problems of 
mathematical physics and the algorithm of its numerical realization are 
essentially interrelated categories. 

The material of section 2.4 develops the arguments for the advan- 
tages of studies of the Frcdhohn integral equations of the first kind in 
the pair of spaces L2 (0, 1) — 12- Thus, an obvious mismatch, in the 
general case, between the free term / (x) and the space I2 allowed us 
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to outline a constructive algorithm of numerical realization by means 
of an adequate representation of the arising error of closure. 

4.3 The methodology of the solution of in- 
correct problems 

Incorrectly formulated problems of mathematical physics arc decep- 
tively transparent from the point of view of the interpretation of con- 
sidered processes. This results from their adequacy to spaces that in the 
computational sense are practically infeasible. If the data of such prob- 
lems are specified in their natural classes of functions, the correspond- 
ing formulations loose a mathematical sense because of insolubility. In 
this situation, of crucial importance is a role of general methodological 
concepts, that is, one has to be guided by a system of global principles. 

From this point of view, J. Hadamard's insistence on the correct 
formulation of problems describing physical phenomena [1] still can be 
interpreted as a kind of hypothesis, whereas in fact related Banach's 
theorem on the inverse operator is a universally accepted element of the 
foundation of modern mathematics [22]. Nevertheless, A. N. Tikhonov 
chose an alternative orientation by a revision of the actual notion of the 
solution of the incorrect problem and the use with respect to it of special 
algorithms of numerical realization [2] . It seems that this choice was in 
full harmony with the attitude of the scientific community of the last 
decades towards the unprecedented revolutionary role of computational 
mathematics in natural science (see [3, 23, 24, 25]). 

A notion of correctness according to A. N. Tikhonov appeared that 
played up a rather simple version of a search for the solution in a re- 
duced class of functions meeting the criterion of the correct formulation 
[14]. By the way, any a priori premises for finding such a class on the 
basis of reasonable information arc missing. 

Shaky character of the conceptual basis led to the failure of the 
idea of a limiting transition with respect to a small parameter in the 
solution of a family of problems related to incorrectly formulated ones 
(the method of regularization [2]). The reason, in our opinion, lies 
in the same inadequacy of the use of functional spaces. Given that 
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I2 is characterized by an infinite number of features, whereas L2 is 
characterized by only one, is it possible, even on the heuristic basis, 
to expect to overcome this cardinal disagreement with the help of the 
regularization parameter al 

The situation in the sphere of activity of numerous followers of A. N. 
Tikhonov looks rather deplorable. The actual efforts are concentrated 
on a surrogate with a small factor a, formed on the basis of (4.1): 

1 

ai;{x) + Jk{x,Oi^{Od^ = fi^), xe[0,l]. (4.2) 


This is called the Fredholm integral equation of the second kind, with- 
out any mentioning of its insufficiency in this respect. Despite a large 
number of investigations devoted to the determination of the regulariza- 
tion parameter a, any more or less constructive algorithms are absent. 
The main reason seems to be the inconsistency of the thesis that implies 
a possibility of efficient matching between the solution and the data of 
incorrectly formulated problems (see, e.g., [2, 26]). 

As a matter of fact, one has to be satisfied only by a comparison of 
solutions to (4.2) obtained in the range of the decrease of a. Because 
of great labor input of numerical realization for small values of the 
regularization parameter, a large-scale application of A. N. Tikhonov 
methodology to the practice of scientific investigations incurred consid- 
erable economic damage. As regards attempts to investigate the Fred- 
holm integral equation if functional spaces of its correct solvability, they 
were isolated and were not accompanied by constructive implementa- 
tion [27]. 

V. M. Fridman, whose papers [28, 29] are considered in section 3.3, 
approached the solution of the Fredholm integral equation of the first 
kind as a mathematically formulated problem regardless of its applica- 
bility to modeling of concrete processes. From the point of view of our 
consideration, the iterative algorithms of V. M. Fridman may be of in- 
terest, because they allow one to achieve maximal possible efficiency in 
the framework of the chosen object of investigation, which is indirectly 
confirmed by their simplicity and brevity. 

In other words, it is hardly possible to obtain more from the tra- 
ditional interpretation of Eq. (4.1) (under the condition that formal 
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convergence exists). By approaching the solution, the employed cor- 
rections become negligibly small compared to the values of the sought 
function: 



As is well-known, in the absence of a timely halt of such a procedure, 
computational "noise" from operations with incommensurate numbers 
can radically distort the solution [5, 15]. It becomes obvious that the 
Predholm integral equation of the first kind, by virtue of its nature, 
contains an inherent defect that principally disagrees with pithy for- 
mulation of the problem of the restoration of the function ijj (x) from 
the right-hand side of (4.1). 

Therefore, it seems to be quite natural to complete the integral 
equation (4.1) by a component that, by adaptively compensating for 
the negative factor of computational "noise", would not substantially 
infiuence the initial mathematical model by virtue of its smallness (see 
the constructive premises in section 2.4). 

In section 3.5, we have given the argument of K. I. Babenko [23] 
for the necessity to take into account the fact of the loss of information 
when evaluating comparative efficiency of computational algorithms. 
This argument seems to be even more important at the stage of the 
formulation of the problem. Since calculations of / (x) from (4.1) ob- 
jectively delete the information on the function ip (x), its restoration in 
the framework of the traditional approach quite naturally reduces to 
an incorrect problem. 

If we hypothetically assume that the sought function ip (x) is ex- 
plicitly present in Eq. (4.1) [i.e., the parameter a of Eq. (4.2) is 
commensurate with the values of the kernel], all the outlined problems 
disappear. Nevertheless, any actions of the type of a coarse substitution 
are absolutely unacceptable. However, such an appearance of ip {x) can 
be viewed in the context of the modeling of computational "noise" with 
the participation of the integral term also. Below, we will specify this 
exclusively important point. 
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4.4 Concepts of numerical simulations at 
the modern stage 

The predetermined method of conjugate gradients is considered to be 
one of the most efficient methods for the solution of ill-posed systems 
of linear algebraic equations that appear as a result of discretization of 
different problems of mathematical physics [30]. The predeterminer, a 
non-degenerate matrix, allows one to reduce the procedure of numerical 
realization to a sequence of algebraic problems with desired favorable 
properties. On the other hand, however, the number of necessary iter- 
ations and their difficulty increase (section 3.7). 

It should be noted that the structure of the predeterminer is de- 
prived of adaptive basis with respect to the orientation at specific char- 
acteristics of a concrete computational procedure. In our opinion, an 
analogous situation is, in general, typical of discrete models. From this 
point of view, the potential of the techniques of continuous analysis is 
at a qualitatively higher level. 

One of the key problems of computational mathematics is the de- 
velopment of the conceptual basis for a relationship between a repre- 
sentation of the data and the efficiency of the employed algorithms. In 
this regard, the ideas of K. I. Babenko [23], completely based on a qual- 
itative interpretation of the notion of information put forward by A. N. 
Kolmogorov, can be estimated as rather pessimistic. Indeed, almost all 
computational operations of this guide are accompanied by a "colos- 
sal" loss of information, whereas rare exceptions correspond only to a 
special representation of initial tables, which, as a rule, is not realized 
in practice. 

The position of R. W. Hamming [31], who is a direct follower of 
the ideas of J. Hadamard in the field of computational mathematics, is 
alternative. In his opinion, methods of numerical realization must be 
adapted to the available information. As regards principal difficulties, 
such as the incorrectness of the formulation, the main attention should 
be concentrated on a modification of mathematical models. The argu- 
ments of P. Bellman and S. Dreyfus for the expediency of the evaluation 
of the quality of information on the basis of its efficiency indices [32] 
are also rather attractive. 
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O. M. Bclotscrkovsky and V. V. Shcnnikov [24] stated a crisis in 
the sphere of numerical simulations resulting from the complexity of 
both the formulations of practical problems and the techniques of their 
numerical realization, (section 3.8). As a reason, they have pointed 
out an inapplicability of methods of "pre-computer" mathematics to 
situations, when owing to the accumulation of round-off errors actually 
any algorithm becomes computationally incorrect. As a matter of fact, 
the authors proposed to develop more intensively approaches in the 
style of A. N. Tikhonov, without any mentioning of the alternative way, 
i.e., matching the formulations of considered problems with Banach's 
theorem on the inverse operator. 

Note that generations of specialists in different fields of mathemat- 
ical physics were brought up under slogans of the type "all real prob- 
lems of the mechanics of continuum medium are ill-posed" that were 
repeatedly reiterated without any explanations by "greats" at different 
conferences. As a result, we have an implementation at a folklore level 
of the thesis supported only by the practice of scientific research. 

N. N. Yanenko, who, in contrast to some colleagues, was well aware 
of the losses of numerical simulations from the breakup of ties of the 
techniques of numerical realization with the basics of functional anal- 
ysis, can be called a flagship of this ideology. However, he considered 
to be of crucial importance the principal difference between classical 
and computational mathematics consisting in the fact that the former 
dealt with abstract symbols without the loss of information, whereas 
the objects of the latter were numerical arrays whose transformation 
was inevitably accompanied by errors of different kinds (see [3, 33]). 

A methodological orientation of the arguments of the works of N. 
N. Yanenko allows us to suggest that a certain role in the formation of 
his ideas was played by ambitious motivations of being a co- participant 
of the emergence of "new" mathematics that, while partly employing 
the "old" one, was, in general, substantially superior. A grotesque 
manifestation of this position is represented by the monograph of A. V. 
Chechkin [25] completely based on A. N. Tikhonov concepts. 

It seems that we facing a distortion of the essence of the problem, be- 
cause Banach's theorem on the inverse operator is an entity of a higher 
level than numerical operations and, at the same time, is indispensable 
to them. Indeed, the boundedness of the inverse operator yields prac- 
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tically a unique possibility to prevent both inadequate dependence of 
the solution on the data and the accumulation of computational errors. 

The outstanding ideas of J. Hadamard and S. Banach who foresaw 
the development of the technique of numerical simulations should be 
regarded as its fundamental basis. 

4.5 Ideas of the development of a con- 
structive theory 

Let us suppose that the kernel k{x,^) of the Fredholm integral equation 
of the first kind (4.1) is closed in L2 (0, 1) and its solution exists.^ We 
also assume a possibility of exact evahiation of the right-hand side by 
means of integration, or the function is assumed to belong to the space 
I2, which implies the condition 

00 \ 

(^n^n < Oin= f (x) Ipn (x) dx, (4.3) 

n=l i 

where A„, ipn (x) are the characteristic numbers and the eigenfunctions 
of the kernel k {x, 

Because of the closure of k{x,^), the solution to (4.1) is unique 
as well; the operator A that maps from L2 (0, 1) into I2 is continuous: 
Hence all the conditions of Banach's theorem on the inverse operator 
are fulfilled. This theorem states that the inverse operator A~^ that 
maps from I2 into L2 (0, 1) is continuous as well. In other words, the pro- 
cedure of the evaluation of ip (x) is stable against small perturbations 
of the given k (x,^) and / (x). Therefore, it can be realized without an 
accumulation of round-off errors of significant digits. 

In this regard, the Inverse World of S. Banach is rather captivating 
and, at the same time, it is based on the absence of any differentia- 
tion of the employed spaces with respect to preference. The dominant 
ideas in the sphere of computational mathematics are purely alterna- 
tive. Therefore, both openly and mainly implicitly, the arguments for 

^In the interests of continuity of our consideration, we repeat some of the above- 
discussed points. 
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non-constructive character of this fundamental theorem are put for- 
ward. 

The arguments are seemingly reasonable. Indeed, the function / (x) 
is determined mostly with an error of measurement. Thus, spaces of the 
type L2 that restrict certain superpositions on an interval are natural 
for its estimate. As already mentioned, the space I2 is illusive because 
it deals with an infinite set of features of the data that, as a result, 
cannot be identified in practice. Even a comparatively small variation 
of / (x) may violate the criterion (4.3) and, as a consequence, distort 
the solution to (4.1) by Picard's theorem. 

Thus, the actually given function / (,x), as a rule, does not belong 
to the space I2. A trivial case of the representation of f {x) as series 
in terms of the elements ipn (x) for exactly determined A„ is an excep- 
tion. However, such premises cannot serve as a basis for the neglect 
of the space I2 in the studies of Eq. (4.1). It seems that constructive- 
ness is possible here only in the context of the agreement of, generally 
speaking, different orientations: 

- the function / (x) G L2 (0, 1); 

- the operator A maps from L2 (0, 1) into I2. 

The motivation is obvious: to preserve the potential of continuous 
inversion of the operator A for practical realization. At the same time, 
the outlined contradiction is clear, and it cannot be overcome exclu- 
sively in the framework of the Fredholm integral equation of the first 
kind (4.1). In this situation, it is quite natural to turn, figuratively 
speaking, to the origin of this equation, that is, to the issues related to 
the formulation of the problem. 

Consider a certain process described by the operator A. The direct 
problem consists in the evaluation of the integral according to (4.1) un- 
der the substitution of the given function ip (x) . If it can be expanded in 
a series of sin (mrx) , in particular, corresponding terms acquire a factor 
of n^^ and, as a result, the convergence strengthens. This procedure 
has a lot of physical and other interpretations and is mathematically 
correct. 

A key element is the formulation of the inverse problem for the 

same operator A, which is related to the restoration of the function 
i(j {x) from given realization of the procedure of integration, that is, 
/ {x). What is actually implied is the determination of the cause from 
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its consequence. Whereas the scheme of the sohition of the direct prob- 
lem is transparent, the status of the inverse problem is diametrically 
opposed. A priority of its solution is the actual algorithmic procedure 
on the basis of an adequate mathematical model that is not an analog 
of the process occurring in the regime of real time.^ 

In general, the traditional formulation of inverse problems by means 
of formal renaming of known and unknown components of mathemati- 
cal models describing objectively occurring processes has no grounds. 

In hght of the above, it is quite natural to recall the statement of 
J. Hadamard that all problems having practical interpretation admit a 
mathematically correct formulation. From this point of view, since the 
solution of the problem, inverse of the evaluation of the integral (4.1), 
objectively exists and is unique, it has to be only correctly posed. At the 
same time, J. Hadamard did mot give corresponding recommendations 
of constructive character, and at the present stage of the development 
of mathematical physics his methodology turned out to be, in essence, 
completely rejected."^ 

Let us try to follow the formulation of the problem, inverse of the 
evaluation of the integral (4.1), that is carried out, in general, with a 
certain error: 

A^^f + Sf, xe[0,l]. (4.4) 

In the direct formulation, taking into account this error has no 
principal importance. Nevertheless, solutions to the Predholm integral 
equations of the first kind (4.1) and (4.4) can be completely different. 
At this point, any quantitative interpretation of Sf is senseless. More- 
over, even under the condition of the existence of an analytical solution 
(Sf = 0), its restoration by means of a discretization reduces to the 
solution of ill-posed systems of linear algebraic equations. 

By general considerations, the presence of 5f' in (4.4) increases the 
potential of the formulation of the inverse problem, and the question 
of a model of the allowed error arises alongside. One must take into 
account that the mechanism of its generation is governed by the factor 

^Indeed, the cause as an outcome of the consequence has no physical sense. 

^The first pubhcation of J. Hadamard on this problem dates back to 1902. One 
may suggest that the idea of incorrectness initiated S. Bamach's investigations in 
1920s. 
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of smoothing of ijj (x) by the integration procedure. Therefore, the 
structure of 6f' must reflect, in a sense, informational incompatibihty 
between this function and its mapping Atp (i.e., between exphcit and 
imphcit representations) . 

By these arguments, we use an operator model of the error in the 
form 



where B is an integral operator, A is a parameter that serve for an ad- 
equate approximation of 6f' in the range from zero to finite distortions 
related to measurements. 

In this way, instead of the Fredholm integral equation of the first 
kind (4.1), we propose to solve the following problem: 



Here, the parameter like A, serves to prevent the Fredholm integral 
equation of the second kind, obtained as a result of transformations, 
from positioning itself on the spectrum, which is equivalent to existence 
and uniqueness of its solution (see, e.g., [34]); 6f = ^5f'. 

Note that the traditional formulation of the inverse problem in the 
form of the Fredholm integral equation of the first kind (4.1) is by 
no means violated, because we have only added a function representing 
zero in L2 (0, 1) to its free term. At the same time, the transformation of 
the incorrect problem (4.1) into the formulation (4.6) creates conditions 
for a radical change of the situation. Indeed, we can demand, generally 
speaking, that 5f adaptively compensate for the errors of numerical 
operations that take f {x) out of the space h- As a result, prospects 
for a constructive realization of the operator A~'^ emerge. For f + Sf ^ 
R{A), the negative factor of the incorrectness of Eq. (4.1) is fully 
neutralized. 



As will be shown below, it is reasonable to represent the operator 
B from (4.5), when 5/ = in L2 (0, 1), as 



6, = I-XB 



(4.5) 



fxA^^fxf + Sf; Sf^O, [0,1]. 



(4.6) 



1 




-1 



imposing the condition of closure on the kernel h {x, C,), transformed by 
means of a linear change of variables. As a result, the problem (4.6) 
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takes the form 

1 1 
k{x,0^{Od^ + >^ J h{x,0^{Od^-l^fi^), (4-7) 

-1 

1 

^l:{x) = \j h{x,0^p{^)di, a;e[0,l]. (4.8) 
-1 

It is easy to notice that the condition on 5/, equivalent to Eq. 
(4.8), is supposed to be satisfied by means of an extension of {x) on 
X e [—1,0), which is equivalent to the use of a new unknown function. 

There exists a well-known opinion that prospects of obtaining new 
substantial results by simple transformation of mathematical relations 
are not great. Indeed, by applying to Eqs. (4.7), (4.8) a subtrac- 
tion operation we again obtain the initial problem which is incorrect. 
However, first, we are not going to do this, and, second, the integral 
equation with the sought function in an explicit form will allow us to 
feel a constructive potential. 

Prom this point of view, a "refusal" of the well-known demonstra- 
tion of smoothing of peculiarities of the sought solution by means of 
integration of (4.1) seems to be very significant. Indeed, assuming that 
the function ijj = ijj^ {x) satisfying the system of equations (4.7), (4.8) is 
known, we give it a perturbation of the type esin {mix). A substitution 
into (4.6) shows that this perturbation influences the free term / {x) 
both via a reduction coefficient and without it, at the expense of 
corresponding integration and an explicit presence of ijj (,x) . 

What is said docs not apply to il-){x). x G [—1,0). However, the 
determination of this function is beyond the scope of the considered 
problem. We want to emphasize that the above arguments are exclu- 
sively heuristic, and, as will be shown, the fulfillment of the condition 
(5/ = is related to some rather subtle points. 

In conclusion of this section, we want to point out the inconsistency 
of the wide-spread opinion that the formulation of problems of numer- 
ical simulation should be left to specialists in applied sciences, whereas 
pure mathematicians should be concerned exclusively with rigorous an- 
alytical investigations, the development of computational methods and 
participation in their realization. 
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It seems that specialists in applied sciences should be concerned 
with the formulation of direct and, generally, correct problems. The 
factor of incorrectness is directly related to the method of the solution. 
Therefore, the main concern of pure mathematicians should be a reduc- 
tion of formulations of problems describing the considered processes to 
the conditions of efficient implementation of Banach's theorem on the 
inverse operator. 
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Chapter 5 



A method of the reduction of 
problems, traditionally 
associated with Fredholm 
integral equations of the first 
kind, to Fredholm integral 
equations of the second kind 

5.1 The structure of the representation of 
the error 

In light of the arguments of section 4.5, we proceed with the consider- 
ation of the same Fredholm integral equation of the first kind 

1 

{Ai;) (x) = Jk {x, iO d^-^f{x), xe [0, 1] , (5.1) 



with a closed kernel k{x,$^). In the determination of the right-hand 
sides by one of methods, including the integration of the initially given 
function ip {x), an error was made Sf/fi, where is a constant. In other 
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words, the property of correctness is inherent to the equation 

f,{A^){x)^f,f{x) + {5f){x), xe[0,l] (5.2) 

in its natural space l2- 

Following (4.5), we present the error as a difference between the 
sought function and the integral component 

(5/)(x)=^(x)-A(5^)(a;), x e [0, 1] , (5.3) 

where A is a constant, by setting^ 

1 

B.^ J h{x,0.d^. (5.4) 

-1 

The substitution of 6f from (5.3) into (5.2) leads to the equation 

(x) = /X (AV) (x) + A (B^l;) (x) -iif(x), xe [0, 1] (5.5) 

with two unknown functions: ip {x) on x & [0, 1] and on x G [—1,0). 
As the second equation of the system, equation (5.3) could be used, 
however, the function {Sf) (x) requires a specification. 

If exclusively the loss of information is implied, resulting from the 
evaluation of / (x) by means of the procedure of integration of (5.1), it 
is logical to assume that 

{Sf){x)^0, xe[0,l], (5.6) 

or 

II'^/IIl.(o,i) = 0. (5.7) 

In this modeling of the error occurs as a result of mutual inad- 

equacy of the components of the right-hand side of (5.3) of a qualitative 
character (in the first place, with respect to smoothness). It should be 
also noted that, except for the trivial case f & h, the substitution (5.1) 
has no mathematical sense, and the problem can be considered only 
from the point of view of a search for an approximate solution: 



min||A^-/||i^(o,i). 
^The use of an indefinite integral in Eq. (5.4) will be explained below. 



(5.8) 
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Taking into account (5.6), equation (5.3) can be given the form 

V' (x) = A (Bi/j) (x), xe [0, 1] , (5.9) 



or 



u 

j h{x,i)^l:{i)di = g{x), a;e[0,l], (5.10) 
-1 



where 



g {x) = h^{x)- jh (x, ^ (0 di. (5.11) 



Making in (5.10) a change of variables C, — 27rx — n; 9 — 27r^ + tt, 
we get 

TT 

\{C,e)^{e)de = g{C), Ce[-7r,7r]. (5.12) 

As is obvious, the satisfaction of (5.9) by ip{x) on x e [—1,0) is 
equivalent to the solvability of Eq. (5.12). Let the kernel h{(,6) be 
closed and g {() G L2 (— tt, vr).^ Then Eq. (5.12) is a Fredholm integral 
equation of the first kind, whose solution, if it exists, is unique (see [1]). 
To simplify the following transformations, it is reasonable to assume 
that its kernel is symmetric, depends on the difference of the arguments 
and satisfies the condition of periodicity on ^ e [— tt, tt] . For example, 

h (C, 0) = TT-r. ^ 77 — 7^^ ^, < Irl < 1. (5.13) 

^ 27r[l-2rcos(C-^)+r2]' ' ' ^ ^ 

In this case, equation (5.12) takes the form of the Poisson integral 
[2]. Correspondingly, in (5.4), 

l-r^ 

^ l-2rcos[27r(x-0]+r2 ^ ^ 

By Picard's theorem, the necessary condition of convergence to the 
solution of Eq. (5.12) in L2 (— tt, tt) of the series 

00 

Y^anXnl/jniC) (5-15) 
n=0 



^Under these conditions, by Eq. (5.11), tjj (x) G L2 (0, 1), which is, by the way, 
implied. 
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IS 



oo r _ 

a'X <oo, g (C) i>n (C) ^^C- (5-16) 

Here, the characteristic numbers and the eigenfunctions of the kernel 
(5.13) have the form,respectively [[1]]: 

A, = r-", n = 0,l,2,...; 

7/;o(C) = l/A Vin(C) = {cos(nC);sin(nC)}, n = l,2,.... 

However, because of (5.11), the free term g{C) and, hence, the co- 
efficients depend on the sought function il){x). Therefore, it is 
impossible to verify the fulfillment of the condition (5.16). From this 
point of view, E. Goursat's remark on Picard's theorem ([3], pages 
141-143) is rather appropriate. It consists in the following. 

The kernel h (C, 9) is assumed to be closed, whereas the function 
g{C,) is such that the condition (5.16) is not fulfilled. Nevertheless, 
there exists a function ip ((^), represented by the convergent series (5.15), 
whose substitution assures the fulfillment of the Fredholm integral equa- 
tion of the first kind (5.12) in L2 (— tt, tt). 

A proof is based on the fact that 

TT 

gn{0^ j h{c,e)iJr.{o)de, (5.17) 

— TT 

where 

n 

iPn{e) = Y,c^i\'4^i{0), (5.18) 

i=0 

coincides with the sum of the first n terms of the Fourier expansion of 
the function g [Q in elements i/^i {(). 

It should be noted that by Mercer's theorem [4] 

=f:M^. (5.19) 

Because of this fact, n can be chosen such that the integral 

TT 

/ b(C)-^n(C)]'^^C (5.20) 
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will be smaller than an infinitesimal e > 0. 

It seems that this result can be considered as a specific addition to 
the Riesz-Fischer theorem with preceding Weyl's lemma [4] . According 
to the classification [5], the series (5.18), by (5.17), converges to the 
solution of (5.12) in the sense of a distribution. 

If the condition (5.16) is not fulfilled, the scries (5.15) diverges in 
L2."^ In this regard, it is interesting to establish restrictions on the data 
of (5.1) to satisfy (5.16). First, one has to make Eq. (5.1) symmetric 
by multiplying by the conjugate operator 

1 



We get: 

1 

Jk'{x,0i^{0d^ = fi^), xe[0,l]. (5.21) 



Here, 

1 1 
k' {x, i) = jk (C, x) k (C, i) dC; f {x) = jk (e, x) f iO dC, 



and, correspondingly, k' (x,^) — k' (^, x). 

Suppose that the function /' (x) is such that the solution of the 
Fredholm integral equation of the first kind (5.21) is a correctly formu- 
lated problem in the space l2- Then 

oo J. 

Y: a'X < oo, < = / /' {x) < {x) dx, (5.22) 

n=0 i 

where and ■0^ (x) are the characteristic numbers and the eigenfunc- 
tions of the kernel k' {x, ^). 

In this case, by Picard's theorem, the solution to (5.21) is deter- 
mined by the series 

oo 

V' (^) = E <Ki^n (^) (5-23) 

n=0 

^This is not equivalent to the conclusion tl) (x) ^ L2 {—1,0). 
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that converges in L2 (0, 1). Correspondingly, in (5.12), 



00 _ " 

^ (C) = E <^'n T< iO- h (C, 9) (9) d9 



and, taking account of (5.19), 



1 -27r 
A 



E cx'm^L / C (C) dC; 





E / (C) (C) rfc, 




= 1,2,.... 



(5.24) 



Thus, the use of the additional condition (5.9) implies that not only 
the condition (5.22) but also (5.16) are fulfilled with the coefficients 



However, in the framework of the formulation (5.8), it is more cor- 
rect to assume that the stated requirements are not fulfilled in principle. 
Correspondingly, the series (5.18) turn out to be divergent for n — > 00. 

At the same time, if the value n, although arbitrarily large, is still 
finite, the situation changes diametrically. Indeed, the conditions of 
the solvability of the Predholm integral equations of the first kind (5.1) 
and (5.21) are fulfilled automatically; the series (5.18) is not merely 
convergent in the sense of (5.20), but is a L2-function. 

We emphasize that these are only heuristic arguments that do not 
allow us to claim that ip (x) e L2 (— 1, 0). For what follows, this condi- 
tion is necessary. Therefore, we turn below to the modeling of the error 
from a somewhat different position. 

In the case of an experimental determination, the right-hand side 
of (5.1), as a rule, differs from the expression for (x), obtained as a 
result of hypothetically exact integration of the given function i/j^ (x), by 
a certain quantity (A/) (x). We only have to require that this quantity 
be bounded, otherwise, on the one hand, the adequacy of the actual 
mathematical model is violated, and, on the other hand, a following 
reduction to the Predholm integral equation of the first kind loses a 
practical sense. 



(5.24). 
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In general, compared to 6f, the error A/ is not of principle impor- 
tance. Indeed, if 

||A/||« 11/11, L2(0,l), 

under the condition of the construction of a stable algorithm of the 
numerical realization, the error of the sought solution ip {x) turns out 
to be comparatively small. If this condition is not fulfilled, both the 
formulation of the concrete considered problem and its solution will be 
merely distorted. 

If we assume in (5.4) 

X 

B. = J h{x,0.d^, (5.25) 

-1 

this fact will influence the modeling of the error via a modification of 
the expression (5.11): 

g (x) = -^[x)- j h (X, i) ^ {i) di. 



The following method of the representation of the error can also be 
considered: 

{5f){x)=il^{x)-\{Bil^){x)-a{x), a:e[0,l], (5.26) 
where a {x) is another unknown function; the operator is given by 

X 

B.^l h(x,0,d^. 



It is obvious that a (x) ^ 0, because otherwise the condition (5.6) 
for the L2-kernel h (x, ^ reduces to the homogeneous Volterra integral 
equation of the second kind 

X 

i,{x) = Jh{x,0i^{0di, xe[o,i], 



whose solution is trivial. 
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5.2 A transformed formulation of the prob- 
lem 

Let us turn to the model of the error given by (5.3), (5.4) and (5.6). 
In order to use more completely the possibilities of the techniques of 
integral equations, we introduce new unknown functions (p (x) and x {^) 

by 



Here, 



* (x) 



^{x), xe [0, 1] ; 

(p{x), xe[-i,o), 



where ip {x) is the solution of the problem in the sense of (5.8); ip (x) is 
determined via ip (x) from the Predholm integral equation of the second 
kind 



ip{x) = X I h{x,0^{Od^ + P{x): a;e [-1,0), 
-1 



with the free term 



J. 

p{x)^xj h{x,0i^{0dC 



and the characteristic numbers A„ = r~", n = 0, 1, 2, . . . (see a transi- 
tion to the Poisson integral in section 5.1). 

In other words, x (•'^) is a combination of the hypothetically known 
solution to the considered problem and of the function ip (x), by means 
of which ^ {x) satisfies an equation analogous to (5.9), however not in 
the mean but uniformly on x e [—1, 0). In this case, equation (5.27) is 
a Fredholm integral eqiiation of the second kind with respect to ^ (x). 

Taking into account (5.1), on the basis of (5.27), we can obtain an 
equation 



* (x) = A {B^) {x) + 



At (A^) (x) - iif {x) + X (x) , xe [0, 1] ; 
0, xe[-l,0), 

(5.28) 
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which, by (5.8), strictly speaking, is neither integral nor of the Fredholm 
second kind type. At the same time, this objection can be easily over- 
come, if, analogously to (5.2), we subtract jj, {5f) (x) in the right-hand 
side of (5.28). 

By the well-known stability of the properties of Fredholm integral 
equations of the second kind against small perturbations [6], the fulfill- 
ment of the condition (5.7) is fully sufficient. Thus, the actual absence 
of Sf in (5.28) will not substantially influence the results of subsequent 
transformations with the use of this equation. 

Using the solution of the Fredholm integral equation of the second 
kind 

*W = A(i5*)W + { "(^^JW^^WlJ,^ P-^l^ (5.29) 

on a; G [—1,0) that will be denoted as (po{x), we can eliminate the 
function x i^) in (5.27). 

As a result, we get an equation of the type 



Here, 




with 



if' {x) = Lp{x) - ipo (x) 



and K (x) being a new unknown L2-function.^ 

In this case, k, ^ 0: otherwise, the solution to Eq. (5.30) would 
be trivial. One can note that a relationship between Eqs. (5.30) and 

(5.27) , (5.28) is reahzed by means of a procedure of a specific "leakage" 
of the function x i^) into k (x) via their solutions on x e [— 1, 0). 

Thus, starting from the traditional formulation (5.1), we have pro- 
ceeded to the Fredholm integral equations of the second kind (5.27), 

(5.28) and (5.30). In light of the general methodological considerations 



^This is determined, in particular, by the type of Eqs. (5.27)-(5.29). 
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of sections 2-4, these equations reveal powerful potential of a construc- 
tive orientation. 

Indeed, the structure of the free terms is such that any solution on 
a part of the domain [x e [0, 1] or x e [—1, 0)] contains explicitly new 
unknown variables x and k, (x) . Simultaneously, in the expressions 
for the solutions on the remaining part, only integral terms are present.^ 
This makes a premise for an analogous representation of 6f, which is, 
in fact, a strategic objective of the proposed method. 

Returning to the questions of section 5.1, note that, on the basis of 
(5.30), we can make an immediate extension of Eqs. (5.5), (5.9): 



11 {A%Ij) (x) - fiq{x), X e [0, 1] ; 
k{x) , a; e [—1, 0) ; 




(5.32) 



In this case, in fact. 



ip' (x)=^(a:)eL2(-l,0). 



In what follows, we will need the inverse of the operator I — XB, or, 
taking into account the change of the variables ( — ttx, 9 — ttx, 



TT 




(5.33) 



— TT 



where 




^We must take into account the form of the solution of the Fredholm integral 
equation of the second kind by means of the resolvent (see, e.g., Ref. [7]). 
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cos (nC) dC 



(r) = [ h (C) cos (nC) dC = I - . o 

— TT 

cos (inC) o?C _ 1 - ^ (^"0 



— TT 

27r 



27r J 1 - 2rcos(C) +r2 

— 27r 



2r cos (C) + ' 



n = 1,2 



It is known ([2], p. 161) that 
l-r2 



2[l-2rcos(C) + r2] 2 



+ ^ r™ cos {m() , 



hence 



2n 



= ^ / 
TT J 



m=l 



COS c^C, n = 1, 2, . . . . 



As a result, A„ = 0.5r n = 0, 1, 2, . . ., and the resolvent of the 
considered operator is [4] 



^ (x, ^, A) - }^ A -A " ^ 1 - 2Ar- 

n=0 n=Q ^ 

where 



ijn{x)^Pn{0, (5.34) 



'j/'o (a;) = 1/v 2, (x) = {cos (2n7ra;) ; sin (2n7ra;)} , n=l,2, ... 



5.3 A constructive algorithm of practical 
realization 

By inverting the operator I — XB in (5.27), (5.29), we get, respectively. 







1 

i^{x)^x{x) + xj H{x,i,\)x{i)di, xe[0,l] (5.36) 
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and 

1 1 
ifo {x) =^iXjl {x, ^0 iOdC + xjH {x, e, A) X iO d^, xe[-l,0), 



(5.37) 

where 

1 

lix,0^ I Hix,C,X)k{C,OdC. (5.38) 



Prom Eq. (5.30), 

1 

«(x)=V'' (^)-A / h(x,0^'(Od^-^J h{x,0^{Od^. 

-1 

Upon substitution of the expressions (5.31) and (5.36), taking into 
account (5.35) and the functional relation [7] 

1 

xjh{x,C)H (C, X)dC^H (x, e, A) - /i (x, , (5.39) 
-1 

we get 



K{x)^\j h (x, (0 di - ipo (x) . (5.40) 

-1 

As the basis for the construction of the system, i.e., of one more 
Predholm integral equation of the second kind in addition to (5.37), we 
will use (5.28): 

1 

x(x)=V'(^)-A / h{x,0'p{Od^-X J /i(x,OV'(0^e 

-1 

1 

-l,lkix,0^iOd^ + l^fix), xe[0,l]. (5.41) 



It is obvious, that the explicitly appearing function ipi^x) cannot be 
represented with the help of (5.36), because x (^) in the right-hand side 
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cancels out.^ Therefore, we employ the solution of Eq. (5.30) 



^Ij(x)^xJ H(x,^,X)K(Od^, xe[0,l]. (5.42) 
-1 

Taking into account (5.40), (5.39), the property of orthogonahty 

j {x) i)m {x) dx = j Ipn {x) V'^ {x) dx = - j ^„ {x) tjjm (x) 



dx 



_ j 0, n ^ m; 
■~ [ i, n = m, 



(5.43) 



that follows from the fact that the kernel (5.14) is periodic not only on 
X e [—1, 1] but also on a; G [—1, 0), x G [0, 1], and, besides, using 

1 1 
J h{x,C)H{C,^,X)dC = J H{x,^,X)h{C,OdC 



that follows from (5.19), (5.34), we find 
1 ? 

^ (^) = J ^) + ^ ^' ^ ^ [°' 1] ■ (^-^^^ 

-1 

To determine the other components of the right-hand side of (5.41), 

we can employ the expression (5.36). Substitution of (5.35), (5.36) and 
(5.44) into (5.41), taking into account (5.39), leads to the equation 

X (.) ^ - / { « , A) + , U, + / * « (c. « . A) dc] U (0 rfe 

OK L } ) 

1 ? 

--A y [/i {x, O + H {x, e. A)] ipo (0 + fif{x), xe [0, 1] . (5.45) 



^In this way, we would trivially return to the initial object (5.1). 
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Note that on the basis of (5.34), in (5.38) and (5.45), 



H {x, c, A) k (c, = 2 E T3|r^^- (^) 

n=o ^ "^"^^ 



J k {x, c) H (c, i,\)dc = 2± 0^i^n (0 , 

n n=0 



where 



C'n iO-Jk (C, V'n (C) dC; < {X)^jk {X, C) V'n (C) d^ 



As a result of symmetry and periodicity H {x, ^, A), x e [0, 1], the 
operator 

1 

7 + A I H{x,C,X),dC 



can be inverted analytically (by analogy with the above-described pro- 
cedure, one can construct the corresponding resolvent). At the same 
time, the realization of this procedure would not yield a considerable 
simplification of Eq. (5.45). 

The inversion in (5.37) of the operator 

1 

I-fiX J I{x,0.dC 



for fi ^ /In, where are its characteristic numbers, leads to the equa- 
tion 

1 r 1 

Vo{x) = xJ H{x,^,X) + fijL{x,C,fi)H{C,^,X)dC xiOdt 

L 

(5.46) 

Here, L (x, ^,/x) is the resolvent of the kernel (5.38). 
In practice, for its construction, the solution can be represented as 
a series 

oo 

V'O (x) = E '^"V'n (x) , 
n=0 
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where a„ are constants to be determined from a system of linear al- 
gebraic equations, obtained by a reduction of factors of the elements 
ipn (x) . In this regard, we mention also a method of numerical realiza- 
tion of the resolvent developed by S. G. Mikhhn ([8], section 12). 

Substitution of the expression (5.46) into (5.45) allows us to re- 
duce the considered problem to the solution of the Fredholm integral 
equation of the second kind with respect to x (^)- 

In particular, in order to give the system of the Predholm integral 
equations of the second kind (5.37) and (5.45) a canonical form, we set 

A = //^ {0.5r-'';r-"}, < |r| < 1; n = 0,l,2,.... (5.47) 

We also introduce the notation 

X (x) = xi (x) , (po (x) = Xi (x) ; (5.48) 



Ku (x,0 = - 



± 

H{x,C,l,) + k{x,0+l^ I k{x,C)H{C,C,l^)dC 



{x, = -^A [h (x, C) + H {x, e, //)] ; (x) = /// (x) ; (5.49) 

1 

K21 {x, = H {x, /i) ; K22 {x,^) = iij H (x, C, A) k (C, i) dC,. 



By virtue of (5.43) and the obvious property 

ipn {x) =iJn{x + 1) , n = 0, 1, 2, . . . , 

equations (5.45) and (5.37) take the form, respectively, 
1 

Xi {x)=fij [Ku {x, Xi iO + Ku {x, X2 m + ^1 (x) ; 


1 

X2 ix)^i^l [i^2i (x, XI (0 + K22 (x, X2 (01 d^; X e [0, 1] . 

" (5.50) 
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As shown in ([4], pages 195-196), this system is equivalent to the 
Fredholm integral equation of the second kind 

2 

x{x)=fijK{x,OxiOd^ + F{x), xe[0,2], (5.51) 



where 



Fix) 



Xi (x), xe [0, 1] ; 
X2{x-1), xe(l,2]; 

f Fi (x) , xe [0, 1] ; 
1 0, xe(l,2]; 



Kn{x,0: [0,1]; 

Ku(x,^-l), xe(0,l), ee(l,2]; 
K2i{x-i,0, xe(i,2), ee[o,i); 

X22(x-l,e-l), a;,ee[l,2]. 



In the case 



M<^, M'=Y. J J \KiAx,Ofdxd^, 
*>i=io 



1 1 



(5.52) 
(5.53) 

(5.54) 



(5.55) 



the solution of Eq. (5.51) can be obtained by means of successive 
approximations [9] 

2 

Xn+i{x)^l^l K{x,OXniOdC + F{x), xe[0,2], (5.56) 



where xo {x) G L2 (0, 2) is arbitrary.^ 
If, along with (5.55), 



± 

J \ Kij {x^^)f < aij — const ^ 



(5.57) 



this convergence is regular. Of course, the parameter /j, cannot coincide 
with characteristic numbers of the kernel (5.54). 

''The so-called method of simple iterations. 
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It should be noted that the procedure (5.56) is, in fact, equivalent 
to the construction of the resolvent [7] 

oo 

r(a;,e,/x) = E/^""'^n(a;,0, 

n=l 

where 

2 

K, {x, i)=K (x, ; Km+1 (x, = 1 Km {x, C) K, (C, dC 



are iterated kernels K {x,C,). 

Iterations of the type (5.56), under the conditions (5.55), (5.57), can 
be directly apphed to the system of the Predholm integral equations of 
the second kind (5.50). For the determination of the appropriate value 
of /i, the means of numerical simulations can be used. The identification 
(5.47) somewhat restricts its variability. 

Note that there are a number of efficient algorithms for the solu- 
tion of Predholm integral equations of the second kind, including the 
determination of spectral characteristics and the evaluation of quadra- 
tures. (The literature on this issue is immense [10, 11, 12, 13].) The 
fundamentals of general theory of approximate methods of numerical 
realization of this type of equations are comprehensively presented in 
[14]. 

It should be emphasized that the structure of the kernel K{x,0 
does not enforce on it any additional features or oscillations beyond 
the framework of the initial formulation (5.1).® 

Given xi {x), from Eq. (5.36) we get 

1 

V'l {x) = xi (x) + H (x, A) xi (e) di. (5.58) 



The function ipi {x) is the solution of the considered problem (5.1) in 
the sense of (5.8), that is, exactly for -0 (x) = V'l (x), min ||/ — /i||i2(o i) 

®In other words, those that result from internal difficulties of the employed 
method. 
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is achieved, where 

1 



As a matter of fact, we determine in this way a maximal possible 
approximation in L2 of the result of integration of (5.1) to the free 
term /(x), containing some errors (of measurement, approximation, 
etc.), that objectively belongs to the space 12- Correspondingly, as 
an extension of (5.8), the solution of the problem (5.1) acquires the 
following meaning: 

PV^-/-5/|L,(o,i) = 0. (5.59) 
This issue seems to be of principal importance. 

5.4 The issue of justification and the mech- 
anism of achieved efficiency 

As follows from (5.50), by virtue of (5.1) and (5.49), the function x (^) 
satisfying (5.51) is also a solution of the system of homogeneous Pred- 
holm integral equations of the second kind 

1 

Xi {x)=iij [K[, (x, Xi (0 + ^12 (x, X2 m di-, 


1 

li J [K21 {x, xi (0 + K22 {x, X2 m d^, X e [0, 1] , 

° (5.60) 

K[, (x, 0=^11 {x,0 + kix,0- 

The solution (5.60) is trivial, if we do not take into account all the 
subtleties related to the character of the dependence of the kernels K[i 



X2 {X) = 

where 
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and K22 on ji whose value does not coincide with the characteristic 
numbers of the homogeneous equation 

2 



where the kernel is determined by (5.54), (5.49), with Kn [x, ^) replaced 
by K[i (j'-C)- This implies a contradiction. 
Let us turn to Eq. (5.27) 

1 

-1 

(5.61) 

assuming that the function ip (x) is explicitly determined via (5.36), 
while the rest of the components are determined as in (5.41). Tak- 
ing into account (5.36) and (5.44), as a result of subtraction from Eq. 
(5.61), we get 

1 . 

X{x)^-pi I H{x,i,tx)x{i)di--ii I [/i(x,0 



-1 
+H {x, e, //)] cpo (0 d^ + f^q{x), xe [0, 1] . 



Here, 



q{x)= n 



fix) -I, J k{x,Oi^i iOd^ 



(5.62) 



(5.63) 



with t/ji (x) being the solution of the problem (5.1), calculated as de- 
scribed in section 5.3. 

Substitution of (5.46) into (5.62), taking into account (5.43) and 
(5.47), leads to the Fredholm integral equation of the second kind 
1 

x{x)=fijQ{x,Ox{Od^ + q{x), xe[0,l]. (5.64) 


Here, 

1 } 

Q (x, = -H {x, e, //) - -/i J [h {x, O+H {x, C, //)] [H (C, e, 1^) 
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1 








with L {x, ^, fi) being the resolvent of the kernel (5.38). 

Suppose that the free term of (5.1) is determined without the error. 
In other words, we are dealing with the situation when Aip^ = /^,. Then 
the algorithm of section 5.3 must lead to ipi = ip^, and as a result of 
(5.63), equation (5.64) becomes homogeneous. 

Contextually, this situation is rather natural, because it reflects the 
equivalence of different representations, namely, of Eqs. (5.36) and 
(5.44), and the function ip (x) on a; G [0, 1]. However, from a formal 
point of view, the solution of the homogeneous equation 



is either trivial or, for /x = with Hn being the characteristic numbers 
of the kernel Q {x,$^), coincides with its eigenfunctions. 

It should be pointed out that, using the assumption of exact ful- 
fillment of (5.1), we cannot consider Eq. (5.65) in an abstract way: 
Its solution should be sought in a narrower class of functions that are 
adequate to the convergent series (5.23). Substitution of the latter 
into (5.65) under the condition of hypothctically exact realization of 
corresponding transformations should lead to an identity. It is logical 
to interpret the cancellation of the parameter /i in the context of the 
relation (5.39). 

As mentioned above, in the general case of approximate determina- 
tion of f (x), equations (5.28) and (5.29) should be complemented by 
the error {5f) (x). Therefore, in light of the discussion of (5.2), we add 
to the free term of Eq. (5.28) and subtract from the same term /iSf. 
By virtue of this fact, correction components appear in (5.37), (5.40) 
and (5.42). Using the technique of transformations of section 5.3 and 
the relation [7] 



1 




X e [0, 1] 



(5.65) 







1 



I H{x,C,f^)H{C,C,fx)dC^d,H{x,^,f,), 



-1 
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we get, respectively, 

1 

(Sipo) {x)^i,Jh{x,^,i,) (Sf) iOd^; 



1 } 

{Sk) {x) = J [h {x, + H {x, e, /.)] {Sf) io d^; 



1 

(^V') {x)^i,lt{x,0 (Sf) iOdC, 



where 

t = [h {x, 0-H {x, e, 1^) - dxH {x, e, 1^)] . 

Simultaneously, the term /iSf appears in (5.63), and, taking into 
account the definition (5.2), instead of (5.65), we get 

1 

x{x)^l^lQ(x,Ox(Od^ + (S^lj)(x), xe[0,l]. (5.66) 



By the use of the notation 

(Sf) (x) = a; (x) (Sx) (x) 
and of the representation 

u (x) (5x) (x) = (5a;) (x) x (x) 

(i.e., the arising correction is achieved as a result of the variation of the 
connecting function cu), equation (5.66) takes the form 

1 

X {x) = /i / [Q {x, + 1 {x, {Sij) m X iO d^, X e [0, 1] . 



In this regard, the algorithm of section 5.3 can be interpreted as 
follows. There exists a function xi {x) that by means of (5.58) satisfies 
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Eq. (5.1) cither exactly or with a certain error. The existence of a 
homogeneous Fredholm integral equation of the second kind with the 
eigenfunction xi (x) is also quite obvious. These considerations allow 
us to draw a conclusion that (5.65) plays the role of such an equation. 

Quite naturally, this is true within accuracy depending on the error 
of the determination of the free term of (5.1).^ In this case, an adap- 
tive connection to the value of the parameter is realized by means of 
variations of Sip, which is equivalent to a shift of spectral characteris- 
tics. In light of the above, the outlined contradiction can be considered 
resolved. 

It is clear that, as a rule, the error Sf and, respectively, the correct- 
ing component 6ip in (5.66), arc unknown, and wc arc actually dealing 
with the Fredholm integral equation of the second kind 

1 

x{x)=i^J[Q{x,0+k{x,0]x{Od^-i^fix), xe[0,l]. (5.67) 



The function x (x), obtained from this equation, allows us by (5.36) 
to calculate the solution of the problem under consideration in the sense 
of (5.59). 

Note that the structure of Eq. (5.65) is determined by the separation 
of the representation of ip (x) on a; G [0, 1]: with the function x (x) in 
the explicit form, Eq. (5.36), and by means of its integration, Eq. 
(5.44). As a matter of fact, this is a key clement of the algorithm. The 
possibility of its realization has been ensured, in the first place, by the 
transformation of the formulation (5.1), namely, the fact that the free 
term of (5.27) is equal to zero on a part of the domain. The latter 
follows from the model of the error (5.3), which is the main concept of 
our investigation. 

Thus, the problem (5.1) is reduced to the solution of the Fredholm 
integral equation of the second kind (5.51) [or the system of analo- 
gous equations (5.50)], whose components are determined according to 
(5.52)-(5.54) and (5.49), as well as (5.14), (5.19), (5.34), under the con- 
ditions that the parameter /i satisfies (5.47) and does not coincide with 
characteristic numbers of the kernel (5.54). To simplify the procedure 



That is, of its misfit to the space h- 
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of the numerical realization, it is reasonable to assume the restrictions 
(5.55), or even (5.57). In Eq. (5.67), the value of ii should not coincide 
with characteristic numbers of the kernel (5.38). 

Summarizing, our arguments for the representation of the error in 
the form (5.3), using the transformed formulation (5.27), (5.28) and 
(5.30) instead of (5.1), have been justified. It should be emphasized that 
before the actual reduction of the considered problem to the solution of 
Predholm integral equations of the second kind the arguments for the 
usefulness of our approach were merely heuristic. 

Prom a practical point of view, of main importance is the modeling 
of the error of the calculation of the function / {x) according to (5.1) by 
a difference between the sought function ip [x) in the explicit form and 
the integral of this function, Eq. (5.4), that is, a qualitatively adequate 
interpretation of the mechanism of smoothing of information. 

The representation of the operator S by a Poisson integral is, of 
course, not the only possible one. However, this point is not crucial. 
Therefore, the expression for the kernel (5.14) can be considered uni- 
versal. Should the importance of the above-mentioned structure of the 
representation of (5.3) be characterized in an ideological and method- 
ological sense, the operator B with the kernel h (x, ^) facilitates the 
implementation of constructive and technical function, which is also 
necessary for the realization of the algorithm. 

In the construction of the Predholm integral equation of the sec- 
ond kind (5.51), any additional information on the data of the problem 
(5.1), except for belonging to the space L2 (0, 1), was not used. There- 
fore, it can be logically considered as being objectively equivalent to 
the formulation that is erroneously supposed to be unambiguously in- 
correct. As far as errors of the determination of / {x) and k (x, ^) are 
concerned, their influence on the solution of both (5.51) and related to 
the formulation (5.1) problems is transformed into the sphere of stable 
dependencies of the theory of operator equations of the second kind 
[14]. 
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5.5 On other versions of the transforma- 
tion 

Applying to Eqs. (5.27), (5.30) the subtraction procedure, we get the 
Predholm integral equation of the second kind 



^ (x) -^'{x)^\jh {x, b (0 - (01 di + K{x), xe [-1, 0) , 
-1 

whose solution has the form 



^ (x) - ^' (x) = k{x) + xJ H' (x, e, A) K (e) d^, X e [-1, 0) . (5.68) 

-1 

Here, the resolvent 

oo j.n _ 

H' (x, e, A) = E J n (x) (0 (5.69) 

n=0 

is obtained by analogy with section 5.1. 
From Eq. (5.30), we get 



ip'{x) = K{x) + xjH'{x,^,X)K{Od^, a;e[-l,0). (5.70) 
-1 

The substitution of Eqs. (5.35), (5.70) into (5.68) leads to the Pred- 
holm integral equation of the first kind 



I [H(x,^,X)-H'(x,^,X)]K(Od^ 
-1 

1 

= Jh{x, e. A) X (0 d^, X e [-1, 0) , (5.71) 



that, nevertheless, allows us to establish a relation between the Fourier 
coefficients of the entering functions: 


j K{x)'lpn{x)dx; Xn ^ j X (x) 'Ipn (x) dx; 71 = 0,1,2,..., 

-1 -1 
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which may seem rather captivating. 

However, the substitution of (5.34), (5.69) into (5.71) yields 

1 - Ar" 

f^n = Xn, n = 0,1,2,..., 

and the expression (5.42) reduces to (5.36). Correspondingly, the func- 
tion X (x) in the left-hand side of (5.41) cancels out and the uselessness 
of the transformations is obvious. 

The reason is that any attempts to represent 5f by use of an explicit 
function for which the Predholm integral equation of the second kind 
is constructed, without any relation to the data of the problem (5.1), 
as it is done above, are futile. 

In (5.3), we replace the operator (5.4) by (5.25) and turn to the 
system of equations (5.32), where the kernel h{x,^) can be different 
from (5.14), with the parameter A = 1. 

It is obvious that 

X 

V'(x)=«(x) + I H{x,OK(Od^, [-1,0), 
-1 

where H (x,^) is the resolvent of the kernel h{x,^), and the problem 
under consideration is reduced to the formulation (5.26) with a new 
unknown function 



a{x)^J H{x,^)^{Od^. 
-1 

The corresponding system of integral equations takes the form 

1 X 

i;{x)^lijk (x, ^(Od^ + Jh (x, V' (0 d^ + (^ (x) - lif {x) ■ 



X 

^Ij{x)^ lh{x,Oi^{OdC + a{x), xe[0,l]. (5.72) 



10 



The use of Eq. (5.29) is implied. 



122CHAPTER 5. A METHOD OF THE REDUCTION OF PROBLEMS, TRADITION. 



Now we use the procedure of successive approximations 



with 



1 X 



X 

= jh{x,i)'il;n{i)di + a{x), 71 = 0,1,2,..., 



V'O {x) = (7 (x) - -Hf {x) , 



to get 



ii){x) = Y, 

n=0 



X 1 

j h{x,i) + ]^ii J k{x,0 



^ (0 - 2/^/ (0 



de, (5.73) 



where the value of is supposed to be small enough to ensure the 
convergence of this series. 

The substitution of (5.73) into (5.72) leads to the cancellation of 
both the function a (x) and the indefinite integrals of this function. If 
we retain only n — N terms of the series, a term of this type, namely, 

X 

j hN+1 {x, o- (0 d^, 


will survive. Here, the iterated kernel has the form 

X 

hN+1 {X, 0^ j hm {X, C) hN-m (C, ^^C, 771 = 1, 2, ... , N. 


Taking h{x,^) such that 

hN+1 {x, x) 7^ 0, X e [0, 1] , 



5.5. ON OTHER VERSIONS OF THE TRANSFORMATION 123 



by differentiating (5.73), under tlie condition that dxk (x, E,) and d^f (x) 
exist, we can get a Fredliolm integral equation of tlie second kind for 
the function a {x). Its structure is determined by evaluating the power 
of the operator sum by the Newton binomial: 

(a + 6)" = Cy + Cla''-^b+.. . + C>"-'"6"^ + . . . + Cr'«^""' + 

(5.74) 

where 



" m!(n — m)! 

To increase the efficiency of the procedure of the numerical realiza- 
tion, it is reasonable to maximize explicitly the multiplier of cr (x), i.e., 
k{x,x), which is achieved at the expense of the strengthening of the 
oscillations of the kernel h {x, ^). This situation manifests itself also in 
the growth of the integral terms, which less important for the terms 
with prevailing coefficients C^. On the other hand, the convergence of 
the series (5.73) objectively worsens. 

In general, for sufficiently large N, the obtained integral equation 
can be represented in the form that does not contain a small multi- 
plier of the sought function a (x) that enters explicitly. At the same 
time, both its actual structure and corresponding transformations are 
rather cumbersome. Prom a constructive point of view, this approach 
is disadvantageous compared to the algorithm of section 5.3. 

Note that, as follows from (5.74), the integral terms of the above- 
mentioned equation have the form 

1 



Here, 

X 

tnm {X, 0^ J K {X, C) km (C, ^^C, (5-76) 


with hn (xjC), km (x,^) being the iterated kernels of the corresponding 
order. However, it is impossible to extract any advantages by means of 
a priori adaptation of /i„ {x,^). 
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The reason lies not only in the difficulties of the restoration of 
hn{x,^) from the integral equation (5.76): In practice, it is impossi- 
ble to find a non-singular kernel tnm {x, ^) for which the integral (5.75) 
could be effectively inverted (see [15]). 

Indeed, a specific exception is provided by Schlomilch's integral 

TT 

2 f 

TT J 


which, however, does not fit into the framework of the suggested scheme. 
We briefiy note one more version resulting from the representation 

X 

a(x)^ lh(x,0i;(0d^ 



in (5.26), where the kernel h {x, ^) is to be determined. 

The problem reduces to the solution of the system of equations 

1 X 

i^{x)^,ijk (x, ^(Od^ + fh {x, V (0 - P^f {x) ; 



X 

i:{x)^ jh {x, V' (0 di, X e [0, 1] , (5.77) 



where, aside from ip (x), the kernel h (x,^) is also unknown. 

The procedure of successive approximations, analogous to that out- 
lined above, for a (x) = 0, yields a solution in the form of the series 
(5.73). Substitution of this series into (5.77) leads to an equation for 
the sought function h {x, ^) outside the integration sign, i.e., of the sec- 
ond kind, however, essentially non-linear. Its higher dimension is yet 
another factor further complicating its solution. 

The model of the error (5.26), analogous to the use of the oper- 
ator (5.25) in (5.3), is more general in comparison with (5.4). This 
implies intuitive attractiveness of the interpretation of Sf in a uniform 
metric. However, the above transformations witness that less rigid con- 
ditions discarding (5.3), (5.4) and (5.6), and also (5.7) turn out to be 
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fully sufficient for their implementation. Therefore, the representation 
(5.26) does not reflect the qualitative peculiarity of the mechanism of 
smoothing of information by integration. 

Nevertheless, the prospects of the reduction of the considered prob- 
lem to the Predholm integral equation of the second kind on the basis 
of (5.26), in principle, can be clearly seen. From the point of view of 
general considerations, this situation is rather remarkable. 
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Chapter 6 



A reduction of linear 
boundary- value and 
initial-boundary- value 
problems to Fredholm 
integral equations of the first 
kind 

6.1 Problems for ordinary differential equa- 
tions 

Consider, for example, 

u" -a{x)u = f{x), xe [0, 1] ; (6.1) 

li' (0) = M (1) = 0, (6.2) 

where a (x) and / {x) are given L2-functions. 
From the notation 

u" {x) =^{x), (6.3) 
129 
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we get 



X 



X 

u{x) = j (x - ^ (0 + cix + Co, 



(6.4) 
(6.5) 



where cq, ci arc the constants of integration. 

The substitution of expressions (6.3) and (6.5) into (6.1) leads to 
a Volterra integral equations of the second kind for a new unknown 
function: 

X 

i^ixj^a {x) j {x-i)'il) (0 di + (cix + Co) a{x) + f{x), xe [0, 1] , 



(6.6) 

whose solution is 

X 

{x) = (cix + Co) a{x) + f{x) + lQ {x, [(ci^ + co) a (0 + / (01 d^, 



(6.7) 

where Q (x,^) is the resolvent of the kernel a{x) {x — 

Taking into account (6.4), (6.5) and (6.7), we find from the bound- 
ary conditions (6.2): ci — 0; 



Co = 






f{0 + IQ{^X)f{C)dC 




di 


i + /(i-0 




«(0 + /Q(e,C)«(C)^^C 




di 



(6.8) 



One can act in a different way: Namely, upon the substitution of 
expressions (6.4), (6.5) into (6.2), we get ci = 0; 

1 



Co = 



{l-i)'^{i)di, 



and, as a result. 



u {x) 



X L 

/(^-0-/(i-0 



^ (0 di. 



(6.9) 
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Thus, the problem under consideration reduces to the Fredholm integral 
equation of the second kind 



ip {x) = a {x) 



1 



/ (^-0-/(1-0 



LO 



(6.10) 



Its solution is 

1 

^(x)^f(x)+ lQ(x,0f(0d^, (6.11) 



where Q (x,^) is the resolvent of the kernel 
—a (x) 



1-e, x<^<l; 
1-x, 0<^<x 



for the parameter A = 1 in the right-hand side of (6.10).^ 

The substitution of (6.7) into (6.5), taking into account (6.8), or the 
substitution of (6.11) into (6.9), allows us to calculate the solution of 
the problem (6.1), (6.2). 

Note that the outlined approach is substantially indifferent to the 
order of the differential equations, in view of initial or boundary con- 
ditions, and to the data of the problem under consideration, such as 
the functions a{x) and f{x). The reduction of the problems to inte- 
gral equations of the Fredholm type, in general, requires an analysis 
of their solvability. In other words, it is necessary to verify whether 
A = 1 belongs to the set of characteristic numbers of the corresponding 
homogeneous equation, which is implemented in the framework of the 
general theory. 

Analogously transformations are traditionally discussed in courses 
of the theory of integral equations (see, e.g., [1, 2]). At the same time, 
as far as the solution of applied problems is concerned, the methodol- 
ogy of the reduction to integral equations of the second kind did not 
gain sufficient popularity, which can be characterized as a kind of a 
paradox. It is rather surprising in light of rather active attempts of its 



^It is clear that this resolvent is not identical to the resolvent of Eq. (6.6). 
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popularization: see, e.g., publications of S. E. Mikeladze, I. A. Birger, 
and A. N. Golubentsev [3, 4, 5]. 

In our opinion, the reasons for this are, on the one hand, purely 
technical difficulties of numerical realization of integral equations be- 
fore wide-spread implementation of computers, and, on the other hand, 
a systematic orientation on the use of the techniques of fundamental 
solutions, which was formed under the influence of general methodolog- 
ical concepts of studies of problems of mathematical physics. Indeed, 
integral equations, as a rule, were constructed on the basis of the the- 
ory of the potential, when the role of the kernel k{x,^) was played 
by the solution for a localized source and the sought function was the 
distribution of the intensity of the corresponding action. 

In this context, it is rather interesting to point out the development 
of the ideas of mutual relationships between the notions of the deriva- 
tive, the indefinite and definite integrals, as it was followed in detail 
by F. A. Medvedev [6]. As a matter of fact, these notions, although 
initially invertible, diverged under the influence of increasing sophis- 
tication of the structure of employed kernels and became practically 
independent. It should be noted that the representation of the solu- 
tion by means of different kinds of local actions was in agreement with 
practical methods in a number of scientific disciplines and was widely 
employed before the appearance of technical methods of efficient real- 
ization of computational operations. 

In general, the idea of a possibility of a comparatively simple reduc- 
tion, with respect to the highest-order derivative, of boundary-value 
problems for ordinary differential equations to Volterra integral equa- 
tions or Fredholm integral equations of the second kind did not receive 
adequate understanding in the field of applications. 

6.2 Illustrations of the procedure of the 
reduction 

Let us turn to a model of the bending of a membrane stretched along 
a contour by a uniform load or, alternatively, of a bar under torsion: 



dlu + d^yU = -l, 



(6.12) 
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u{0,y)=u{l,y) = 0, 
u (x, — u (x, 1) = 0. 
Introducing the notation 

we get 



(6.13) 
(6.14) 

(6.15) 



u {x, y)^ j {x-0^ it y) + xgii (y) + gu (y) , (6.16) 



where gij (y) are functions of integration. 

In view of (6.15), equation (6.12) takes the form 



9yU (x,y) = -l-'ijj(x,y), 



and, respectively. 



u 



1 ^ 

{x, y) = --y"^ - j (y-v)^ {x, v) dv + y92i (x) + ^22 (x) , (6.17) 



where g2j (y) are also functions of integration. 

The substitution of expressions (6.16), (6.17) into the boundary 
conditions (6.13) and (6.14), respectively, allows us to determine gi2 — 
922 = 0; 

1 ^ 1 

5ii(y) = -/ {i-0^{ty)d^; 92i{x) = 1^ + j {l-v)^{x:V)dv- 



As a result. 



u (x, y) 



X J. 

j{x-i)-xj{l-i) 



i^{C,y)dC; (6.18) 



u(x,y) = -y(l -y) - 



j {y -ri) - y j {l-rf) 



ip{x,r))dr). (6.19) 
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Eliminating u {x, y) from these expressions, we get a Fredholm in- 
tegral equation of the first kind with respect anew unknown function: 



X i 



+ 



y 1 

/ iy-v)-y j i^-v) 



ij{x,r))dr) = -y{l -y) . (6.20) 



Thus, a principal difference from the one- dimensional case consists 
in the reduction of the problem (6.12)-(6.14) to an incorrectly formu- 
lated one. However, here we will be interested not in a numerical real- 
ization of Eq. (6.20) (note that the algorithm of the previous section 
applies to it as well) but in the universality of the procedure of trans- 
formation. 

Indeed, let the domain of the problem be different from the canonical 
one, and let, for example the second condition (6.13) have the form 
u (7i y) = Oi where x — ^ (y) a certain single- valued function. Instead 
of (6.18), we have 



u {x, y) = 



Jix-O-x J biy)-C] 



^{ty)dt 



and any specific changes are absent, because a transition to an ordi- 
nary procedure of the evaluation of the integral requires only a non- 
orthogonal mapping of the type x = '-fx, y = y. 

It is not difficult to notice that each of expressions (6.18) and (6.19) 
satisfy identically the pair of boundary conditions (6.13) and (6.14), 
respectively. The rest of the conditions are fulfilled approximately, 
depending on the accuracy of the determination of ip {x, y) . At the 
same time, the solution can be represented in the form that satisfies 
identically both the conditions (6.13) and (6.14): 



Ui {x, y) = ui {x, y)-{l- y) ui {x, 0) - yui {x, 1) ; 



U2 (x, y) = U2 [x, y)-{l- x) U2 (0, y) - xu2 (1, y) 
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Here, the functions Ui {x,y), U2 {x,y) are determined by (6.18) and 
(6.19), respectively. 

The norm of the error of closure of the values of Ui {x, y) or Ui {x, y) 
allows us to estimate the error of the approximate solution: 

^^ 2\\U^{x,y)-U2{x,y)\\ 
\\Ui{x,y) + U2{x,y)\\ ' 

However, if instead of (6.13) the conditions 

d:cU (0,2/) = d:cU (1,|/) = 

would be imposed, they could not be satisfied by the expression for the 
derivative 

X 

dxu {x, y)^ j (C, y) + 9ii (y) 



that follows from (6.15). 

Nevertheless, this complication can be easily overcome by the use, 
in particular, of the relation 

dlu + u — 'iIj 

that allows as to retain both the functions of integration gij (y). 

In the general case, it is reasonable to turn to an equivalent formu- 
lation of the problem (6.12)-(6.14): 

d^ui + dyU2 = -I; ui{x,y) = U2{x,y) , (6.21) 

ui (0, y) = ui (1, y) = U2 (x, 0) = U2 (x, 1) = 0, (6.22) 
using a representation of the solution of the type 

(x, y)^ I ki (x, y, V'l (C, y)d^ + Yl l^ij (x) Qij (y) ; 
J=i 

y 2 

U2 {x, y)^ k2 {x, y, rj) ip2 {x, rj)drj + J2 /^2j (y) 92j (x) . 
n J=i 
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We assume that the kernels are given and satisfy the conditions 
h {x,y,x) = ki {x,y,y) = 0; 

d^ki{x,y,x) ^0; dyk2{x,y,y) ^ 0, x,ye[0,l]; (6.23) 

fiij {x), fi2j {y) are also given; gij (y), g2j (x) are to be determined from 
the boundary conditions as discussed above. 

Thus, in agreement with the previous scheme, we can set /in = x, 
/^2i = y, A*i2 = A*22 = 1- The conditions (6.22) lead to the expressions 



ui {x,y) = 

U2 (x, y) = 
and, respectively. 



X 1 

j ki {x,y,i)-x j /ci(l,|/,0 



1 

k2 {x,y,rj) - y j k2 {x,l,ri) 



i^i{C,y)dC; (6.24) 



ip2 {x, rf) drj, (6.25) 



a^/ci (x, y, x) V'l (x, y) + j dlki (x, y, V'l (C, y) = dl^i (x, y) ; 

° (6.26) 



dyk2 {x, y, y) V^2 {x, y) + J dyk2 {x, y, Vj) V^2 {x, rj) drj = dyU2 (x, y) . 

° (6.27) 
Let, in addition to the conditions (6.23), diki {x, y,^) and dyk2 {x, y, rj) 
arc L2-kcrnels. Taking into account the a priori information about the 
solution of the considered problem, (6.26), (6.27) are Volterra integral 
equations of the second kind with respect to the functions ipi {x,y), 
ip2ix,y), whose solutions, according to the general theory, exist and 
are unique. Therefore, expressions (6.24) and (6.25) are adequate to 
the physical content of the problem (6.21), (6.22). 

Note that, for a different choice of the elements Hij (x), ii2j {y), 
integral equations constructed in this way could belong to the Predholm 
type of the second kind. In this case, ki {x,y,^), k2 {x,y,ri) must also 
satisfy a solubility condition with the parameter A = — 1, which, in 
fact, does not pose any problem. 
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Taking into account (6.23), we can set 
^1 {x, y, = - K (x, y, ; ^2 (x, y, r])^{y- 7]) (x, y, rj) , 
where 

k[{x,y,x) 0; k'^^x.y^y)^^, x,yG[0, 1], 

using expressions k'^ {x, y, x), k'2 {x, y, y) to refract the a priori informa- 
tion about the solution in order to smooth the sought functions tpi (x, y) 
and, in general, to simplify the procedure of calculations. It is clear that 
this point is important for more complicated problems with different 
singularities of the solution, and we mention it here only for the sake 
of completeness. 

The substitution of expressions (6.24), (6.25) into (6.21) produces 
a system of integral equations 



^2 {x,y) 



where 



dyk2 {x, y, y) 



d^ki {x, y, Vj) i)2 (x, r})dr) + F {x, y, ipi) , 

(6.28) 



F{x,y,'ijji) = 



dyk2 {x, y, y) 

X 

+ 1 8% (x, y, e) ^1 (e, y) di 



[1 + d^ki {x, y, x) ipi {x, y) 



X 1 

j h {x,y,^)-x j A;i(l,y,0 



y i- 

j h {x,y,r])-y J /cs (x, 1,77) 



Prom Eq. (6.28), we find 



ip2 {x, T}) dfj ^0, x,y e [0, 1] . 

(6.29) 



y 

fp2 {x, y)^F (x, y, -01) + y <5 {x, y, ??) F (x, 77, -^i) d-q, (6.30) 
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where Q (x, y, rj) is the resolvent of the kernel — <9^^2 {x, y, rj) /dyk2 {x, y, y). 

Substitution of the expression (6.30) into (6.29) leads to a Fredholm 
integral equation of the first kind with respect to the function ipi {x, y). 

The described procedure can be directly extended to differential 
equations of other types. As an illustration we consider the simplest 
problem of thermal conductivity: 

dtu - dlu2 = 0, (6.31) 

u{x,Q) =Uo{x)] 'u(0,t) = M(l,t) = 0. (6.32) 
From tfj = d^u, equation (6.31) and conditions (6.32), we get 



u {x, t) 



u 



X ± 

J {x-O-x J (1-0 i^{U)d^; 

.0 



As a consequence, 

- X in i 

j {x-C>-x j {l-i) il^i^i, t) di-j ijj (x, r]) df] = Uo (x) ; x,y e [0, 1] 

.0 J 

In order to make an analogous reduction of the problem of the 
bending of a rectangular plate of variable stiffness D, fixed along a 
contour, we write [7]: 

DAAu + 2d^Dd^Au + 2dyDdyAu + ADAu 

- (1 - I/) {dlDd^yU - 2d,yDd,yU + SFyDdlu) = g, (6.33) 

a> (0, y) = 9> (a, y) = d'^u {x, 0) = d'^u {x, 6) = 0, n = 0, 1, 

(6.34) 

where A = d1 + d"^; v is the Poisson coefficient; q (x, y) is the intensity 
of the transverse load. To calculate the derivatives with respect to x, 
we can set 

X 4 

u{x,y)^ I k (x, y, V' (C, y)di + Y^ x^'^gij (y) . (6.35) 
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Here, 

d^h{x,y,x) =0, n = 0,l; dlk2 {x,y,y) ^ 0, a; G [0, a] ; ye[0,b], 

and the functions gij {y) are intended to satisfy the conditions (6.34) 

for X = 0, X = a. 

The second representation of the solution via ip {x, y) is determined 
by means of the substitution of (6.35) into Eq. (6.33) and four-fold 
integration over the variable y. The appearing functions g2j {x) allow 
us to satisfy the conditions (6.34) for y = 0, y = h. After that, u (x, y) 
is eliminated from the representation of the solution.. 

Note that with the help of a special structure of the kernel k {x, y, rj) 
one can easily satisfy conditions of the type u{xi,yi) = at isolated 
points inside the considered domain. The procedure of the reduction 
also applies to mixed boundary conditions ( a change of the type along 
a side) and to the case of a connection of plates. Analogously, three- 
dimensional problems of mathematical physics can also be reduced to 
Predholm integral equations of the first kind. 

6.3 Universality and analogous approaches 

Thus, a rather elementary method of the reduction of linear boundary- 
value and initial-boundary-value problems to Predholm integral equa- 
tions of the first kind is practically universal from the point of view of 
its realizations as far as the following aspects arc concerned: 

- the order and structure of differential equations; 

- the form of boundary conditions; 

- the availabihty of variable coefficients; 

- the form of the domain; 

- the dimensionality of the problem. 

In this situation, all the information about a concrete problem is 
transferred into a functional equation, whose solution does not require 
any conditions on the contour of the domain, which poses a substantial 
advantage. Thus, its solution can be sought in the form of a series 
of a system of coordinate elements intended exclusively to ensure the 
efficiency of the procedure of the numerical realization. 
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However, the problem obtained as a result of transformations is 
incorrect, hence its numerical realization requires adequate methods. 
At the same time, in applications, the solution of such a problem can 
be acceptably approximated by a series with the number of terms that 
does not affect the stability of the numerical algorithms. Therefore, 
one can hardly explain the absence of interest to a systematic use of 
this procedure, especially in the period before the general orientation 
at the discretization of problems of mathematical modeling. 

One may state that special literature did not point out the existence 
of a formalized method of the reduction of practically arbitrary initial- 
boundary-value problems to Frcdholm integral equations of the first 
kind. At the same time, there a number of examples of applications 
of analogous transformations in rather particular situations. As a rule, 
they were a given physical interpretation that considerably disguised 
the generality of this approach. 

Thus, Yu. V. Repman used as new unknown variables boundary 
forces of a plate of a canonical configuration that allowed one to satisfy 
conditions on an internal contour of complex configuration [8]. L. A. 
Rozin has developed a method of separation that admits a reduction of 
the problems of calculations of membranes under the forces of interac- 
tion of isolated bars to systems of Fredholm integral equations of the 
first kind ([9], section 9). Some publications point out the advantages 
of the approximation of higher-order derivative of differential equations 
with respect to one of the variables that, compared to numerical dif- 
ferentiation, are much more accurate (see [10]). It should be noted, 
however, that an actual transition to an incorrect formulation, as a 
rule, passed unnoticed. 

In general, we think that there occurred a kind of assimilation of 
the discussed procedure of the reduction by the methods of the the- 
ory of the potential, based on the use of integral relations along the 
boundaries of the domains and by the techniques of fundamental solu- 
tions [11]. The reduction of the dimensionality of the sought functions, 
achieved in this way, seems to have overweighed by its importance the 
above-mentioned universality. Moreover, the construction of integral 
equations with strong singularities in the kernels that partly smooth 
over the factor of incorrectness attracted certain attention [12]. 

Some problems for differential equations, and, in particular, the 
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following one: 

dxyU = adxU + hdyU + cu + f, 

where a, b, c and / are given functions of the variables x and y, can be 
reduced directly to Volterra and Fredholm integral equations of the sec- 
ond kind with respect to the higher-order derivative {ip — d^yu). These 
issues are studied in detail by G. Miintz [13]. Of considerable interest 
is the fact, established by this author, that analogous transformations 
cannot be extended to the case of the simplest equation of the elliptic 
type. 

6.4 A connection to the algorithm of the 
previous section 

Fredholm integral equations of the first kind, appearing in the above- 
discussed case of the reduction of two-dimensional boundary- value prob- 
lems, can be represented in the form 

1 1 

{All)) (x, y) = Jti (x, y, V' it v) j ^2 {x, y, rf) i/j (x, r])dr]^ f (x, y) , 



(6.36) 

x,ye [0, 1]. 

A typical discontinuity of the kernels ti {x,y,^), T2 {x,y,r)) on the 
diagonal ^ — x, r) — y, as well as other analogous complications, are of 
no principal importance. 

The algorithm of section 5 docs not require in this case any substan- 
tial changes. The differences are related only to the structure of the 
operator A. The analog of the system of Fredholm integral equations 
of the second kind (5.37) and (5.45) takes the form 

1 1 
cpo {x, y)^fxXjH {x, e. A) (Acpo) (C, y)di + \j H (x, A) x {i. y) dC; 


(6.37) 

1 

X y)^-\ j H (x, A) X (C, y)d^-i^ (M) {x, y) 
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J + H (x, A)] ipo (e, y) + l^f {x, y), x,y e [0, 1] , 



where 



^ (x, y) = x{x,y) + xj H {x, e, A) X (C, v) (6-38) 



Under the condition (5.47), taking into account (5.43) and the no- 
tation (5.48), we can reduce the problem to the solution of a system of 
two dimensional Predholm integral equations of the second kind: 



Xi {x,y) 



1 1 

j Kn (x, y, Xi (C, y)d^ + j {x, y, 77) Xi {x, v) dr] 



1 1 



+ y j {x,y,^,rj)xi d^drj + j Ku (2;, 0X2 (tv) d^ 







+Fi{x,y); 



X2 (x, y)= n 



1 

K21 {x, Xi (0 + / i^22 {x, y, X2 (e, I/) 



1 1 



+ j J K23 {x, y, 77) X2 (C, ?7) c^Cc^?7 







, X, ye [0,1]. 



Here, 



Ku {x, y, = - 



1 

-f^ {x: t /^) + Ti (x, y, + / (2;, y, C) -f^ (C, f^) dC 



K12 {x, y, rj) = -T2 (x, y, rj) ; X13 (x, y, v) = -7-2 (2;, y, r^) /7 (x, //) ; 

-fCu (a;, = (^> + (^> /^)] ; -^1 {x, y) = l^f {x, y) ; 

1 

K21 (x, i)^H (x, ^, ; is:22 (a;, y,i)^ n j H (x, (, )u) ti (C, y, c^C; 



-f^23 (a;, y, ^) = /^-f^ {x, ^, IJ) T2 (C, y, v) ■ 
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In addition to (5.47), the parameter /i should not coincide with 
the characteristic numbers of the above system of the equations. To 
simplify the procedure of the numerical realization, it is reasonable to 
choose this parameter in agreement with conditions of the type (5.55), 
(5.57). 

The function ip {x,y) satisfying Eq. (6.36) in the sense of (5.59) is 
determined from (6.38) with x(a;, y) = xi {^iV)- It should be noted 
that the techniques of the numerical realization of Predholm integral 
equations of the second kind, pointed out in section 5.3, can be directly 
applied to this case. One can also use special methods of the evaluation 
of two-dimensional integrals ([14], section 3). 
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Chapter 7 



Other classes of problems 



7.1 The initial-boundary- value problem for 
the Korteweg-de Vries equation 

Consider the formulation 

dtu - Gud^u + dlu = 0, (7.1) 

u{x,Qi) ^uq{x)] u{0,t) ^ ui{t) ; dxu{0,t) ^U2{t); u{l,t) = U3{t) , 

(7.2) 

where the given functions Uq{x); Ui(t), i = 1,2,3 are supposed to 
impose no additional restrictions on the model of the process under 
consideration. 

It should be noted that up to now there is no any general theory 
that allows one to investigate a priori the solvability of the problems of 
the type (7.1), (7.2) in adequate classes of functions. In this situation, 
the main means of the specification of the mathematical model are re- 
sults of numerical simulations as well as solutions of specially simplified 
equations near the boundary (see [1], section 10). 

Using the procedure of the previous section, it is not difficult to 
reduce the problem (7.1), (7.2) to an integral equation of the first kind 
with respect to 

ip {x, t) = dlu {x, t) , 
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which yields 

{x, t) = ^J{x-Cfip (e, t) + -x^gs (t) + xg2 (t) + gi (t) , (7.3) 



u 



with the functions of integration determined from the boundary condi- 
tions: 

gi (t) = ui (t) ; g2 (t) = U2 (t) ; 

1 

gs (t) = 2 [m (t) - U2 (t) - m (t)] - I (1 - 0' V' (e, t) dC. 



Substitution into (7.3) leads to the expression 

1 



u {x, t) 



1 



.0 



^{^,t)d^ 



Wuz (t) +X{1-X)U2 (t) + (l - x'^) Ui (t) . 

Now we rewrite Eq. (7.1) in the form 

dtu = GudxU — d^u. 



(7.4) 



(7.5) 



The substitution of (7.4) into the right-hand side of (7.5) and in- 
tegration from to t under the initial condition (7.2) allows us to 
determine 



t V X 



u{x,t)=6j J {x-O' + x'J (1-0' 



LO 



X 

_0 



1 

-^/(i-O 



V' (C, V) di + 2x1^3 (77) + (1 - 2x) U2 (77) > dr] + [2x1^3 (r?) 



+ (1 - 2x) U2 iv) + 2xui 



X 1 



^ iC, V) dC > dr) 



- y ip {x,r}) drj + 6 j x^us (77) + x{l — x)u2 (rj) + (l - x^^ ui (77) 



7.1. THE INITIAL-BOUNDARY-VALUE PROBLEM FOR THE KORTEWEG-DE VRIES EQUATIC 



X [2x^3 (7;) + (1 — 2x) U2 {rj) — 2xui (77)] dr] + Uq {x) . 



(7.6) 



Eliminating u{x,t) from (7.4), (7.6), we get the equation 



{Aip){x,t) = f {x,t) , x,t efl:0 <x,t <1, 



(7.7) 



where the expressions for the operator A and the free term / are rather 
obvious. 

Let the data uq (x), Ui (t) be such that / (x, t) G L2 (fl) and, besides, 
the solution of the problem (7.1), (7.2) be a L2-function. In this case, 
following section 5.1, we can use the relation 



where the operator B and the parameter A are the same. 

The algorithm of section 5.1 allows us to reduce the considered 
problem to the solution of a system of two nonlinear integral equations 
of the second kind (6.37), where the variable y should be replaced by t. 

By virtue of continuous dependence of the integrand of (7.6) on 
■0 {x, t) and the principle of contracted mappings for sufficiently small 
values of the parameter /i satisfying the condition (5.47), the solution 
of this system of equations is achieved by means of simple iterations 
with an arbitrary initial approximation from L2. An estimate of the 
influence of the error of calculations, the issue of the acceleration of the 
convergence, as well as other aspects of the numerical realization are 
discussed in [2]. 

It should be emphasized that the realization of the conditions of 
convergence with the help of the parameter n, which, in fact, allows 
one to determine particular solutions by comparatively simple means, 
is of special importance in the nonlinear case. It is implied that other 
solutions can be determined in the framework of numerical simulations 
for large /i, related to studies of implicit functions, branching and bi- 
furcations. 

Using expressions (7.4), (7.6), the solution can be given a form that 
satisfies identically both the initial and boundary conditions (7.2). We 
can also estimate the error by analogy with the discussion in section 



7/j{x,t) ^ X{B7p){x,t) , x,teQ, 



6.2. 
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7.2 A boundary- value problem for a sub- 
stantially nonlinear differential equa- 



Hcrc, we discuss nonlinearity related to higher-order derivatives. As an 
example, consider Monge- Ampere's equation that plays an important 
role in a number of problems of geometry: 



where Si, i — 1, 2, 3 and q depend on the variables x, y, the sought 
function u {x,y) and its first derivatives dxU, dyU [3]. 

For definitiveness, we set Sj — Si{x,y), q = q {x, y) and 



Without going into analysis of the conditions imposed on the data 
of this problem to ensure the existence and uniqueness of the solution, 
we outline a scheme of a numerical simulation. It should be noted that 
corresponding estimates are rather non-trivial and, as a rule, cover only 
certain particular cases [4]. 

We use the notation 



tion 



dludyU - {dxyuf = sidlu + S2dyU + s^dr^yU + q, 



(7.8) 



u (0, y)^u (1, y)^u {x, 0) = w {x, 1) = 0. 



(7.9) 



d^u {x, y) = Vi {x, y) ; dlu (x, y) = -02 {x, y) ; 

d^yU {x,y) = 'tp {x,y) , x,y E Cl : < x,y < 1. 
As a result, taking into account (7.9), we get 



X 



1 




Lo 







y 



1 





Lo 







X 



y 













7.2. A BOUNDARY-VALUE PROBLEM FOR A SUBSTANTIALLY NONLINEAR DIFFERENTIAL 



Upon the substitution of these expressions into (7.8) and the ehm- 
ination of the function u {x, y), we reduce the problem to the following 
system of equations: 

-01 y) -02 {.X, y) - -0^ (x, y) = Si (x, y) -01 (x, y) 

+S2 {x, y) iIj2 (x, y) + S3 (a;, y) i) (x, y) + q (x, y) ; (7.10) 

' X 1 "I X y 

j{x-i)-xj{l-i) MU)di- jdiji:{i,r))dr^^Q- (7.11) 





X y 



iy-v)-y j i^-v) ^2{x,r])di]-Jd^Jtl){^,i])dri = 0, x,yen. 







(7.12) 



The existence of d^u, dy U implies a possibility of two-fold differenti- 



ation of Eqs. (7.11), (7.12) with respect to x and y, respectively, which 
yields 



Equation (7.10) takes the form 



d^ip {x, T}) dr] 



dyip (C, y) 



y 

-ip"^ (x, y) = si {x, y) j d^ip (x, r]) dr] 



X 

+S2 (x, y) j dyip (^, y) d^ + S3 (x, y) ip (x, y) + q (x, y) , x,y eil, 



and after integration in the limits 0, x and 0,y: 

{M){x,y) = f {x,y) , x,yen, 

where 



(7.13) 



X y X y 

(Alp) (x, y) = J tp (C, y)d^ (x, rj)dri- Jtp (^, rj)d^ j ^p (x, 77) ^77 
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y V 



- j ^ {i^y) di j ^ (C, r]) dri-j drj j [si (x, vi) ijj ({, Vi) - si (0, rj) ijj (0, m) 



X 1 X ^ 

- J ^^s^ (e, V) ^ (e, ^i) di dm- Jd^J [S2 (e, y) ^ (Ci, - S2 (e, 0) (6, 0) 







y X y 



f {x,y) = J d^ J q r;) dr?. 







The above implies the existence of the derivatives Si, dyS2. 
As far as a numerical reahzation of Eq. (7.13) is concerned, the 
considerations of the previous section can be applied. 



7.3 Nonlinearity of the boundary condi- 
tion 

Consider a typical problem of the irradiation of an infinite plate with a 
thermally insulated surface into a medium whose absolute temperature 
is equal to zero [5]: 

ardlu - dtu = 0, (7.14) 
u{x,0) ^uo{x); pd^u{0,t)+u"'{0,t) ^0; d^u{l,t)^0. (7.15) 

Here, u {x, t) is the temperature gradient; Uq (x) is a given function; 
qt is the temperature conductivity; p = A/a, with A, a being the 
thermal-conductivity and the heat-transfer coefficients, respectively. 

Introduce the notation 

dluix,y)^,ljix,t), (7.16) 

which leads to 

X 

u{x,t) = J {x-0^{C,t)dC + xgi {t)+g2{t), 





7.4. A SMALL PARAMETER BY THE HIGHEST-ORDER DERIVATIVE OF THE DIFFERENTIAl 



where Qi [t) are functions of integration. 
The boundary conditions (7.15) yield 

1 

gi{t) = -Jij{U)d^; 9iit)=P9Tit), 



and, respectively, 



u {x, t) 



{x-^) -X 



Using (7.14), (7.15) and taking into account the initial condition 
(7.15), we get 

t 

U (x, t) = ttT {x, Tj) drj + Mo (x) , 



and the problem under consideration is reduced to the solution of the 
nonlinear integral equation of the first kind (7.7), where 



A, 



11 



j {x-i)-x j di + [p j ,di\ -ttT j ,dri; 

.0 J , V 

/ {x,t) = Uo (x) . 



7.4 A small parameter by the highest-order 
derivative of the differential equation 
of the problem 

As an illustration of general considerations, wc consider the problem of 
heat transport induced by the processes of thermal conduction and con- 
vection (the first and the second terms, respectively, of the differential 
equation) [6]: 



dtu = ed^u + pdxU. 



(7.17) 
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Here, /3 > is a constant; e is a small parameter, 

u{0,x)^0; u{t,0)^0; u{t,l) ^ Ui{t) , 



(7.18) 



with Ml (t) being a given function. 

The notation (7.16) under the boundary conditions (7.18) leads to 



u {x, t) 



/ (^-0-/(1-0 



i;{^,t)d^ + m{t). (7.19) 



The integration of (7.17) in the limits 0, t with the use of (7.19) and 
of the initial condition (7.18) yields 



u {x, t) 



L 



eip {x,ri)+i3 ip ({, ri) 



df]. 



The problem is reduced to the linear integral equation of the first 
kind (7.7), where 

A^Ai + A2, 

t X t V X 1 " 

Au^eJ.drj; A2. ^ (3 J d^ J . d?] - J (x - ^ - J - 

Lo . 

f{x,t)^ui {t). 

The algorithm of section 5 allows us to go over to a Predholm integral 
equation of the second kind 

X {x, t)^ix [(ei?i + R2) x] {x, t) + F {x, t) , (7.20) 

where Ri and R2 are corresponding operators. 
As a result of the expansion [7] 

00 

X i^: t) = e™Xm {X, t) , 

m=0 

we get a sequence of recursion relations 

Xo {x, t)= n {R2X0) {x, t)+F {x, t) ; 
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Xi {x, t)= fJ' (-R2X1) {x, t)+iJ, (-R2X0) {x, t) ; 

Xm+l (a;, ^) = (i?2Xm+l) (a;, ^) + (-R2Xm) (2^, t) 

that are integral equation of the same type of the canonical structure. 

The above allows us to state that the proposed approach is rather 
efficient in the problems of mathematical physics with a singular pertur- 
bation, whose numerical realization, as a rule, meets with considerable 
difficulties (see, in particular, [8]). Indeed, we managed to transform 
the singular perturbation (7.17) into the regular one (7.20), which fa- 
cilitated a radical simplification of the problem.^ 

7.5 Equations of a mixed type 

Boundary-value problems for equations of this type are characterized by 
non-triviality of the investigation of the issues of existence and unique- 
ness (see [9]). As a consequence, one has to consider such equations 
on rather special domains, which restricts the field of practical appli- 
cations. 

As an illustration, consider well-known Tricomi's equation 

ydlu + d^u = 0. (7.21) 

this equation belongs both to the hyperbolic and elliptical types for 
y < and y > 0, respectively. We restrict ourselves to the framework 
of a numerical simulation, that is, we do not consider a priori the solv- 
ability of the problem and employ, for example, the following boundary 
conditions: 

u{0,y)=u{l,y) = u{x,-l)^0; u{x,l) = u{x), (7.22) 

where the function u (x) and satisfies the conditions v (0) = p (1) = 0. 
Introduce the notation 

^ (x, y)^^ (x, y) , (7.23) 

^According to a standard classification, singular and regular perturbations af- 
fect, respectively, main and dependent terms of the operators. 
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which, taking into account (7.22), leads to 



u {x, y) 



X 1 

/ (x-O-xJil-O 



Two-fold integration of Eq. (7.21) in the limits —1, y under the 
conditions (7.23) and (7.33) yields the expression 



u{x,y) = 



J {y-v)- / (1 - ^) 

.-1 -1 



rjijj {x,r)) dr)+- {l + y)u{x) 



The problem is reduced to the solution of the Predholm integral 
equation of the second kind (7.13) on the domain Q : < x < 1,-1 < 
y < 1 with the operator 



A. 



.0 J , L-1 -1 



ri,dr] 



and the free term 



f{x,y) ^ ^{^ + y)j^{x) ■ 



Note that the so-called condition of matching on the line of parabolic 
degeneracy y = 0, imposed on the solution of Eq. (7.21) ([9], p. 27), is 
fulfilled in a natural way: 



hni u {x, y) = hm u{x,y), x e [0, 1] ; 



lim^ dyU (x, y) = Ymi^ dyU {x, y) , xe [0, 1] . 

As in the previous subsection, this situation results from the fact 
that the singularity of the problem is transferred from the main term 
of its operator to the dependent one. 



7. 6. THE INVERSE PROBLEM OF THE RESTORATION OF THE COEFFICIENT OF THE DIFFL 



7.6 The inverse problem of the restora- 
tion of the coefficient of the differen- 
tial equation 

Small oscillations in the transverse direction of a stretched string of 
variable density are described by the equation 

d^u = a {x) dlu. (7.24) 

Here, t are dimensionless coordinates; 

a{x) =NT'^/p{x) /^ 

with N being the tension, p [x) the density of the material, 21 the length 
of the string, T the time interval. 

We assume that the ends of the string are fixed, whereas its density 
and the oscillations are symmetric with respect to the coordinate x — Q. 
The corresponding boundary conditions have the form 

d^u{Q,t) = 1^(1,^) = 

We also employ the following initial conditions: 

u (x, 0) = uq {x) ■ dtU (x, 0) = 0. (7.26) 

The coefficient a {x) is to be determined from (7.24)- (7.28) for given 
Uq{x), N, I, T and additional information on the oscillations of the 
middle cross-section of the string: 



7.25 



u (0, t)^u (t) . (7.27) 

It is known [10] that the solution of this problem can exist and 
be unique in the classes of functions that are fully adequate to the 
method our consideration. We assume that the necessary, in this sense, 
requirements are fulfilled. 

By analogy with what was done many time before, using the nota- 
tion (7.16) and (7.24)-(7.26), we find the representations 



u {x, y) 



X J. 

j[x-i)-xj{l-i) 



(7.28) 
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u [x, y) = a [x] j [t — ifjil) (x, 7]) dr] + uq (x) 





By eliminating u {x, t), we obtain an equation of the type (7.7). The 
substitution of (7.28) into (7.27) leads to the integral equation 

(Ai^) (t) = /i (t) , te[o,i], 

where 

1 

^. = /(i-0.rfe; fi(t)^-u(t). 



One can estimate the convergence of successive approximations to 
the solution of the thus formulated problem by the smoothing of the 
dependence of the coefficient a on the values of t in the process of 
calculations. 



7.7 The problem of the Stefan type 

Consider the classical model [11]: 

dtu = dlu, < x < 7 (t) ; < t < 1, (7.29) 

u {x, 0) = uo {x) ; u{0,t)^u (7 {t) ,t) ^0, Uq (0) = 0. (7.30) 

On the moving boundary that separates the phases an additional 
condition is imposed: 

ad,u{j{t),t)^^'{t), 7(0) =70, (7.31) 

where 70 > 0; the constant a can be both positive and negative; 7' (t) = 
0^7 (t) /dt. 

Thus, the data of the problem are uq (x), a and 70; the functions 
u {x, t) and 7 {t) are to be determined. 

In Eqs. (7.29)- (7.31), it is reasonable to make a non-orthogonal 
mapping 

X = x/7 (t) , t = t (7.32) 
on the canonical domain Q : < x,t < 1. 
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We get: 

dtu - [x-f' (t) /7 (t)] dsu = dlu, 
m(x,0) = 0; M(0,f) = = 0; 

ad^u(l,t)=-i' {t), 7(0) = 7o- 
By analogy with the above, the notation 

^x'^ y) = V' y) , 
conditions (7.34) and equation (7.34) lead to 



(7.33) 
(7.34) 
(7.35) 



u (x, y) 



{x-i)-x / (1-0 



LO 



V^(e,t)c^^; (7.36) 



u{x,y) 



^(x,77) + [x7' {t)h{t)] 



.0 



+1*0 {x) . 

The substitution of (7.36) into (7.35) yields 



(7.37) 



7' {t) = aJ^iP {^,t)d^, 



from which we get 



{i) = a j ^d^ j ({, r/) cir/ + 70, 



(7.38) 





and, respectively, in the expression (7.37), we have 

xi it) 



lit) 



ax J Ci^{C,t}dC/ a J CdC lij{C,v)dv + lo 



Loo 

The elimination of u{x,i) from (7.36), (7.37) leads to the integral 
equation of the first kind (7.7), whose solution was discussed above. 
The function ■0 determined from this equation should be approximated 
by an analytical dependence on x in order to make an inverse change 
of variables. The sought boundary is determined from the nonlinear 
integral equation (7.38). Finally, using Eqs. (7.36), (7.37) and (7.32), 
we calculate the function u {x, t). 
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Chapter 8 
Conclusions 



Let us briefly summarize the above consideration. Thus, the solution 
of the Predholm integral equation of the first kind 



in the "convenient" for the numerical realization space L2 is, without 
any doubt, an incorrectly formulated problem.In the case of the space I2, 
the situation is diametrically different: the data of equation. (8.1) may 
satisfy the conditions of its correctness, but, nevertheless, the solution 
determined by Picard's theorem will constitute a series that diverges 
as a result of the accumulation of round-off errors. 

It should be noted that even an actual verification of whether the 
data of (8.1) belongs to the space I2 is, as a rule, infeasible. More- 
over, as an objective factor of incorrectness, there appear the error of 
experimental determination of / (x) as well as different errors of the 
identification of the system under consideration, i.e., k (x,^). 

As ma matter of fact, the basis of our consideration is formed by the 
suggestion to connect adaptively Eq. (8.1) to the space I2 by means 
of a functional model of the allowed error, including the smoothing 
of information by the procedure of integration characterized by the 
condition 



1 



(Ai;) {x) = jk {x, ^ (0 di^f{x), xe [0, 1] (8.1) 







II^/IIl2(o,i) 



0. 
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Starting by heuristic considerations that were later supported by 
more firm arguments, we showed the usefulness of the representation of 
the error as a difference between the sought function and of the integral 
component: 

1 

(Sf) {x) = (x) - A (Bij) (x), B.^Jh {x, 0. dC, 

-1 

in order to satisfy (8.2) by (x) , a; G [1, 0). 

We developed the proposed concept by transforming the formulation 
of the problem under consideration: 

*W = A(B*)W + {^W-^^-,lM; (8.3) 

(8.4) 

*'(x) = A(B*')W + {,(°)/,"i°jil;0). (8.5) 

Here, 

^^""^-{vix), XG[-1,0); ^ ^""^ \^'{x), xe[-l,0). 

Note that the solution of the Fredholm integral equation of the 
second kind (8.3), expressed via the resolvent of the kernel h{x,^), 
ensures, in fact, the achievement of the set goal, because it contains 
on X e [0, 1] and [—1, 0) the function ip (x), respectively, explicitly and 
under the integral sign: 

1 

i^ix)^xix) + xjHix,C,X)xiOdC, ^e[0,l]; (8.6) 



1 

ip{x)^xlH{x,C,X)xiOdC, xe[-l,0). (8.7) 
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Later, using the data of the problem (8.1), which is important, 
we expressed by means of integration of x (^) the function ipQ (x) that 
constitutes the solution of the Predholm integral equation of the second 
kind (8.4) for / (x) = 0, as well as 



k{x)^X j h {x, ^) cpo (0 - (x) . 

-1 

From the Fredholm integral equation of the second kind (8.5), a new 
representation of the sought solution via x (^) emerged: 



i:ix) = xj H{x,^,X)K{Od^, xe[0,l]. (8.8) 
-1 

The measure of non-identity of this representation to (8.6) is re- 
flected exactly by the error {Sf) (x). Therefore, the solution of the 
problem (8.1), understood in the sense 

\\A^ - / - 5/11^^(0,1) = 0' (8-9) 
could be represented by expressions of principally different structure: 
with the function x (x) that enters explicitly, equation (8.6), and with- 
out this function, equation (8.8). 

This factor predetermined the reduction of the considered problem 
to the solution of the Predholm integral equation of the second kind 
1 

x{x)^iJij[Q (x, - k' (x, 0] X (0 di + /// (x) , xe [0, 1] , (8.10) 



where 

1 

k' (x, i)^k{x,i) + xjk (x, i) H (C, A) dC (8.11) 



by the substitution of expressions (8.6)-(8.8) into (8.4). 

The structure of Eq. (8.10) allows us to give the following inter- 
pretation. As a result of (8.6), (8.11) and (8.9), the obtained solution 
must also satisfy the homogeneous equation 

1 

x{x)^l^lQ{x,Ox{OdC, xe[0,l], (8.12) 
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with accuracy that depends on 6f, which is quite natural This means 
that in the process of transformation we build up, in a way, the kernel 
Q {x, ^) for which the solution of the problem is an eigenf unction, with 
equation (8.12) positioned on the spectrum. 

Thus, the traditional formulation (8.1) is transformed into equation 
(8.10), whose advantage lies in the fact of continuous dependence of 
the solution on the data of the problem, an additional advantage con- 
sists in a possibility to choose the parameter fi from the condition of 
convergence to the solution by simple iteration 

1 

Xn+l {x)^pLj [Q (X, - k' (X, 0] Xn (0 + l^f{x), XE [0, 1] . 


Note that it proved to be sufficient to allot the kernel h{x,^) the 
property of completeness, whereas its symmetry, dependence on the 
difference of arguments and periodicity were also very important and 
facilitated analytical transformations. 

We have shown that wide classes of problems of numerical simula- 
tions can be rather elementarily reduced to Predholm integral equations 
of the first kind. After that, the discussed procedure of the correct for- 
mulation and numerical realization can be applied to them without 
any substantial changes. By the way, a traditional difference between 
direct and inverse formulations of problems of mathematical physics 
practically vanishes. In the situation of nonlinearity, the realization of 
iteration algorithms of the solution to the obtained integral equation 
by means of the parameter fi is of special importance. 

In light of the above, we can draw a conclusion that, if the phe- 
nomenon (process) admits an adequate description by methods of nu- 
merical simulations, the restoration of its underlying cause or of differ- 
ent parameters from an objectively sufficient volume of additional infor- 
mation does not pose principal difficulties, because the corresponding 
problems can be correctly formulated. From this point of view, an anal- 
ysis of actually observed events, including multi-factor social-economic 
and ecological processes, can be done with much larger efficiency. 

Maybe, it would be reasonable to suggest that, in general, the pro- 
cess of the understanding of the World is much simpler than a wide 
audience usually supposes it to be under the influence of the sphere of 
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applied science that, at present, armed with means of electronic pro- 
cessing of information, constitutes, in fact, a natural monopoly with an 
almost dominant role of commercial component and, correspondingly, 
a systematic drive for investment? 

Thus, colossal means are invested in problems of the restoration of 
dependencies from empirical data and, in particular, in remote probing 
of the surface of the Earth by spacecraft. What is actually realized is a 
search for minimally and maximally acceptable values of the parameter 
a in the integral equation of the type 

1 

a'ilj{x) + J k{x,Oi^{OdC-fix), [0,1]. 



The essence lies in the necessity to establish a balance between 
computational and, respectively, financial abilities of the solution of 
an ill-posed system of linear algebraic equations (which implies a real- 
ization of one of discretization methods) and an approximation of the 
formulation of the problem under consideration to an exact one that is 
erroneously associated with the factor of incorrectness for a = 0. 

In this regard, we note that, of course, it would be incorrect to 
suppose that problems in science are altogether absent or that one can 
develop, irrespective of the circumstances, efficient means to overcome 
these problems. However, in our opinion, complications of principal 
character are inherent, in the first place, to direct formulations of some 
problems, that is, to the construction of mathematical models of insuf- 
ficiently studied processes and phenomena. 

It is clear that the solution of some classes of inverse problems may 
also pose substantial difficulties, but, nevertheless, the wide-spread the- 
sis that the procedure of the restoration of the cause from the conse- 
quence is incorrect, in general, seems to be manifestly erroneous. 

J. Hadamard's statement that the problems that adequately de- 
scribe real processes are correct is an ingenious idea, whose construc- 
tive development allows one to attain a qualitatively higher level of the 
potential of methods of numerical simulations. 

March 21, 2001. 
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Appendix A 

A method of successive 
approximations for Fredholm 
integral equations of the first 
kind 



Theorem. Let K (x, s) be a symmetric square integrable positive 
definite kernel, and let the equation 



be solvable. Then the sequence {ipn (2^)}, determined by the recursion 
relation 



V. M. Fridman^ 




(A.l) 



a 



(:r) + A[/(x)-/„_i(x)], 



(A.2) 



where 



(/?o {x) e L2 (a, h) , 




(A.3) 



a 



^Uspekhi, Mat. Nauk 11, No 1, 233-234 (1956). 
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< A < 2Ai, (A.4) 

and Ai is a minimal characteristic number of the kernel K [x,s), con- 
verges in the mean to the solution of Eq. (A.l). 
Proof. Let us set in Eq. (A. 2) 

multiply both its parts by an eigenfunction of the kernel Vi (x) and 
integrate over x from a to &: 

b b 

a" = Q;"~^ ~ J ^'^ ^"^-^ dx J K {x, s) Wn-i (s) ds, 

a a 

where 

b 

a" = Un (x) Vi (x) dx. 



As a result of the fact that K {x, s) is symmetric and Vi (x) satisfies the 
equation 



Vi (x) — Xi j K (x, s) Vi (s) ds — 0, 



we have: 

b 



j Vi (x) dx j K {x, s) Un-i (s) ds — j (x) dx j K {x, s) Vi (s) ds 



, n-1 



^ [ Un-1 {x) Vi (x) dx 
Ai J Ai 



Thus, 

a 



By virtue of the completeness of the system of functions Vi (x) , 

J ul (x) dx = Y: , (A.6) 



i=l 
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which is to be estimated. 

The scries EZk «)^ = T,Zk (l - j-f'' i^if is majorized by the 
series J2iLk i'^i )^^ because, by the inequahty (A. 4), 

Therefore, for any given small positive e, there exists k > K (e) that 
does not depend on n, such that 

oo 
i=k 

At the same time, by choosing n > N (e), we can make 

k k / \ \ 2" 



,2 e 
<2' 



because for finite i we have the inequality sign in the formula (A. 7). 
Thus, we arrive at the inequality 



j ul (x) dx < e, (A.8) 



which proves the Theorem. 

Submitted on March 29, 1954. 



